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ABSTRACT 


A simplified pol: energetic cell theory is formulated 
to determine spatially averaged energy dependent thermal 
fluxes in the moderator, cladding, and fuel regions within 
the unit cell of a reactor lattice. The derived spectra 
are then utilized in the calculations of the thermal integral 
parameters and average cross sections required for reactor 
computations. 


The cell theory, as formulated, postulates an infinite 
moderator region with the absorption cross section of this 
region appropriately modified to account for the neutron 
leakage into and absorption by the fuel element. The mod- 
ifications to the moderator absorption cross section are 
formulated both in terms of the net current at the fuel 
element-moderator interface and in terms of energy dependent 
moderator and fuel element escape probabilities, the latter 
approach offering physical transparency and ease of cal- 
culation. 


Analytic expressions for the escape probabilities are 
presented, integral transport theory being applied to the 
fuel element region, while diffusion theory is utilized in 
the moderator region. Using these analytic expressions, 
the theory is applied to actual lattices in the form of the 
light water moderated and uranium dioxide fueled cores of 
the Cornell University Zero Power Reactor. Room temperature 
parameters and their temperature coefficients are determined 
using both the monatomic gas model and the Nelkin water 
kernel to describe the energy transfer process in the moder- 
ator. Calculations are made with PROGRAM COUTH, a Fortran-63 
program written for use with the Control Data Corporation 
1604 digital computer. A typical lattice calculation 
including the computation of the spatially averaged fuel, 
cladding, and moderator spectra and the thermal integral 
properties and average cross sections takes approximately 
thirty-five seconds of computer time. This figure is exclusive 
of the compiling time and the time required to calculate the 
moderator scattering kernel. 


In an attempt to estimate the accuracy of the calculational 
results, the method is applied to the Brookhaven National 
Laboratory uranium dioxide cores and the results are then com- 
pared with those predicted by Honeck's THERMOS code. Dis- 
advantage factors agree to within 1.0% while the thermal 
utilizations agree to within 0.5%. A study of the sensitivity 
of the calculated integral parameters to variations in the 
input data leads to the assignment of rather small uncertainties 
in the results calculated with the simplified cell theory. 
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CHAPTER lL 
INTRODUCTION 


The analysis and design of any thermal reactor requires 
an accurate knowledge of the various and competing neutron 
reaction rates within the reactor assembly. This, in turn, 
implies a knowledge of the distribution of thermal neutrons 
in space and energy. In the calculation of the thermal neu- 
tron spectra one must consider the two processes which affect 
these distributions in a heterogeneous reactor lattice. The 
first of these is the moderation process wherein a neutron 
is scattered by a moderator atom or molecule with a subsequent 
change in energy. The second process is the transport or 
diffusion mechanism whereby the neutrons are carried into 
regions of relatively different absorption characteristics. 
Although the time independent Boltzmann's equation or the 


1/ 


integral transport equation— correctly describes both the 
moderation and transport phenomena, the difficulty in 
obtaining solutions to these equations have led previous 
investigators to develop their work along one of two lines. 


Initial attempts to characterize the distribution of thermal 


ter apae Tey ae. 
neutrons were done for infinite homogeneous pee eee 


L/ A. M. Weinberg and E. P. Wigner, The Physical Theory of 


Neutron Chain Reactors, pp. 221-234, The University 


of Chicago Press, Chicago, 1958. 


2/ M. J. Poole, M. S. Nelkin, and R. S. Stone, "The Measure- 
ment and Theory of Reactor Spectra,'' Progress in Nuclear 
Energy, Series I, Vol. 2, pp. 91-164, Pergamon Press, 
New York, London, Paris, and Los Angeles, 1958. 
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Chap. 1 2 


there is no net transport or diffusion and the spatial 
variable is of no consequence. This approach allows con- 
centration on the moderation or energy transfer process 
and sophisticated mathematical models may be incorporated 
into the calculation. 


The other approach, as exemplified by multigroup 


4/ 


diffusion calculations,— attempts to consider spatial 


effects at the expense of the details of the moderation 
process. In these calculations, the scattering kernel 
is often characterized by a simple model which allows 
ease of mathematical description but which may lack 
accuracy in its approximation of the energy transfer 
process. 
More recently, interest has been focused on the 
use of transport eon for dealing with the heterogeneities 


present in a reactor lattice. This approach is embodied 


in the work of signal where both the spatial and 


3/ Neutron Spectra in Lattices and Infinite Media, 

Proceedings of the Brookhaven Conference on Neutron 
Thermalization, BNL 719, Vol. II, Brookhaven National 
Laboratory (1962). 

4/ A. Hassitt, "Methods of Calculation for Heterogeneous 
Reactors,'' Progress in Nuclear Energy, Series I, Vol. 2, 
pp. 2/1-313, Pergamon Press, New York, London, Paris, 
and Los Angeles, 1958. 

S/ Neutron Spectra in Lattices and Infinite Media, see 
Footnote 3. 


6/ H. C. Honeck, ''THERMOS -- A Thermelization Transport 
Theory Code for Reactor Lattice Calculations" BNL 5826, 
Brookhaven National Laboratory (1961). 
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moderation effects on the spectra are considered through a 
space and energy calculation utilizing integral transport 
theory with only slight restrictions on the complexity of 
the scattering model used. 

Still another approach is that found in the work of 


8,9/ In this formulation the effect of moderation 


Leslie. 
and heterogeneity are both considered, but the spatial 
effects are mathematically subordinated by consideration 

of only the spatially averaged fluxes in the moderator and 
the fuel. This method of attack leads to physical trans- 
parency while allowing one to maintain and concentrate on 
the importance of the moderator scattering model. 

This work is an extension of Leslie's ideas, culmina- 
ting in a computer code which provides for the rapid cal- 
culation of the spatially averaged thermal neutron spectra 
in the moderator, fuel, and cladding regions of a reactor 
lattice. The simplified cell theory utilized in these cal- 
culations is formulated in Chap. 2. A general formulation 
is presented as well as a formulation in terms of escape 
probabilities defined for the fuel element and the moder- 
ator regions. The escape probabilities utilized are a 
polyenergetic generalization of those defined by Amouyal, 
iI H. C. Honeck, "The Calculation of the Thermal Utilization 

and Disadvantage Factor in Uranium/Water Lattices," 

Nucl. Sci. Eng. 18, 49-68 (1964). 

8/ D. C. Leslie, "The Spectrox Method for Thermal Spectra 


in Lattice Cells,'' AEEW-M211, United Kingdom Atomic 
Energy Authority (1962). 








Chap. l 4, 


10/ 


Benoist, and Horowitz— in their one velocity thermal 
utilization theory. Chapter 3 deals with the escape prob- 
abilities defined for the moderator region, while Chap. 4 
details the calculational methods available for the deter- 
mination of the fuel element escape probability. Mathe- 
matical models for the energy transfer or scattering kernel 
are discussed in Chap. 5. Emphasis is placed on the models 
actually used for calculational purposes; the monatomic gas 


ILL ab) 13/ 


kernel and the Nelkin water kernel. Chapter 6 


defines the calculations necessary in the determination of 
spatially averaged, energy dependent thermal fluxes for the 
moderator, cladding, and fuel regions in a lattice utilizing 
the simplified cell theory. Integral thermal properties are 
defined as are the averaged cross sections necessary for the 
accurate calculation of reaction rates. In Chap. /7 the 
formulated theory is applied to actual lattices. A study 

is made of the low enrichment, high density, uranium dioxide 


fueled and light water moderated Cornell University Zero 


2/ D. C. Leslie, 'The Calculation of Thermal Spectra in 
Lattices,'' Proceedings of the Brookhaven Conference 
on Neutron Thermalization, BNL 719, Vol. II, pp. 592- 
609, Brookhaven National Laboratory (1962). 

10/ A. Amouyal, P. Benoist, and J. Horowitz, "Nouvelle 
Méthode de Determination du Facteur d'Utilization 
Thermique d'une Cellule,'' J. Nucl. Energy 6, 79 (1957). 

it/ E. P. Wigner and J. E. Wilkins, Jr., "Effect of the 
Temperature of the Moderator on the Velocity Distri- 
bution of Neutrons with Numerical Calculations for H 
as the Moderator,'' AECD-2275 (1944). 








Chap. l ») 


14/ 


Power Reactor cores.—' Temperature effects are predicted 
through the use of temperature dependent models for the 
scattering kernel. The accuracy of these results are dis- 
cussed in Chap. 8. This chapter also includes a comparison 
of the results predicted with the simplified cell theory to 
those calculated by Honeck with THERMOS=2/ for the Brookhaven 
light water moderated, and uranium dioxide fueled lattices. 


Chapter 9 contains a summary of the work undertaken as well 


as the conclusions which may be drawn from it. 


i2/ J. E. Wilkins, Jr., "Effect of Temperature of the Moder- 
ator on the Velocity Distribution of Neutrons for a 
Heavy Moderator,'' CP-2481 (1944). 

13/ M. S. Nelkin, "Scattering of Slow Neutrons by Water," 
Phys. Rev. 119, 741-746 (1960). 

14/ D. D. Clark, "'The Cornell University Nuclear Reactor 
Laboratory,'' CURL-1, Cornell University (1961). 

15/ 


H. C. Honeck, see Footnote /. 








CHAPTER 2 
A SIMPLIFIED CELL THEORY 


2.1 General Formulation 


Recognizing that the neglect or oversimplification of 
the spatial effects in a heterogeneous lattice leads to 
answers which have little relationship to the physical 
processes taking place, one is prompted to overemphasize 


the spatial variations by resorting to multigroup die etedone od 


ef 


or transport—' codes. However, these codes consume a con- 
Siderable amount of computer time, particularly if the more 
complex models of the scattering kernel are used. 

Thus, there is a need for a simple polyenergetic 
formulation for determining the thermal distribution of 
neutrons in the various regions of a heterogeneous lattice. 
Hopefully one could retain the calculational ease of the 
homogeneous, infinite medium thermal spectrum computation, 
thereby offering only slight limitations on the complexity 
of the scattering model used, while still accounting for 
the heterogeneities of the lattice in a meaningful way. 

This section formulates a simplified cell theory for 


use in the calculation of thermal spectra and parameters 


in a lattice composed of uniformly arrayed unit cells. The 


16/ R. S. Varga, ''Numerical Methods for Solving Multi- 


Dimensional Multigroup Diffusion Equations," 


Proceedings of S osia in Applied Mathematics, 
Vole: Se . T6Z~ 180" Nuclear Reactor Theory, 


American Mathematical Society, Providence (1961). 


ay 
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heterogeneous effects due to neutron transport within the 
cell are formulated in a physically transparent manner, 
thereby allowing sufficient emphasis to be placed on the 
moderation process in any calculation. 

Consider first the neutron balance equation to be 


solved in the case of an infinite, homogeneous, moderating 
Pedium: 22! 


E 
e 


| 2g) | BE) | o'(E) = S(E) +[ S(E'-~E)o "(E')dE'. 
O 
(2219) 

E. is the thermal cutoff energy above which the moderator 
flux per unit energy, o (E), is assumed to be 1/E, while 
=(E'-~E) is the energy transfer or scattering kernel 
expressed in Gan. 5, (E) and 5° (E) are the macroscopic 
absorption and scattering cross sections respectively of 


the moderator, 5, (E) being defined by 
E 


c 
5. (E) -| (Edis (22) 


O 


iz! H. C. Honeck, "THERMOS -- A Thermalization Transport 
Theory Code for Reactor Lattice Calculations," 
BNL-5826, Brookhaven National Laboratory (1961). 

18/ M. J. Poole, M. S. Nelkin, and R. S. Stone, "The Measure- 
ment and Theory of Reactor Spectra,'' Progress in 
Nuclear Energy, Series I, Vol. 2, pp. - , Pergamon 
Press, New York, London, Paris, and Los Angeles, 1958. 








Sec. 2.1 g 


The source of thermal neutrons, S(E), due to down scattering 


from energies above Eo is given by 
OO 


S(E) | 3(E'>E)o0(E')dE' ; (2.3) 


Be 


where 9, (E) is pPrOporeional co i) Ee 


Neutron conservation requires that 
E E 


C Cc 
[ S(E)dE | 5, (E) o(E)dE , (2.4) 
O O 


which is obtained by integrating Eq. (2.1) over energy from 
0 to Ea: 

The problem as stated places no restrictions on the form 
of the scattering kernel, =(E*—»E), however, we require that 


it satisfy detailed balance: 
=(E'—~E)M(E') = =(E—>~E')M(E) , (225) 
where the Maxwellian distribution, M(E), is given by 


M(E) = E/(kT)* exp(-E/kT) . (2.6) 


The moderator temperature in electron volts is equal to kT, 
T being the moderator physical temperature while k is 
Boltzmann's constant. 

There are of course practical limitations which restrict 
the complexity of the kernel and the number of energy points 


to be utilized in the computer solution of Eq. (2.1). 
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The infinite medium calculation is now extended to the 
heterogeneous lattice utilizing the ideas embodied in Leslie's 


SPECTROX method .22/ 


Since the only effect which the fuel 
element has upon the spectrum in the moderator of a lattice 
cell is to modify the transport process in the vicinity of 
the fuel element, while contributing very little to the modera- 
tion process itself, one may postulate an infinite moderator 
region with a contrived energy dependent absorption cross 
section which includes the heterogeneous properties of the 
lattice. In other words, if the absorption cross section of 
the moderator can be appropriately modified to account for the 
neutron leakage into and absorption by the fuel element, one 
can then turn his attention to the choice of a proper 
scattering kernel. The calculational ease of Eq. (2.1) is 
still available along with the fact that the non moderating 
part of the reactor lattice has been accounted for. 
Considering the energy dependent absorption rates, one 
defines 


absorption rate in the fuel element 


7 at ener E er eV 
£(E) = absorption rate in the lattice cell ° (2.7) 


at energy E, per eV 
The "fuel element'' includes the fuel rod and its cladding, 
and in the case where there is no neutron absorption in the 
cladding, f(E) can be considered to be the energy dependent 


thermal utilization. In terms of f(E) the neutron 


19/ D. C. Leslie, "The SPECTROX Method for Thermal Spectra 
in Lattice Cells,'' AEEW-M211, United Kingdom Atomic 
Energy Authority (1962). 








Sec. 2.1 10 


loss rate per unit volume of the moderator due to leakage 
into and absorption by the fuel element may be expressed as 
SA(E)9 (E)£(E)/[1-£(E) J, where © (E) is the spatially 
averaged flux in the moderator region having a volume of 


yo: 


m 


< 


OE) = + | eee : (2.8) 


Equation (2.1) may now be rewritten for the cell in a 


heterogeneous lattice as 


| Be (E) + 5. (E) E + 2 os | ® (E) = S(E) + 
Ee 
if (EE) (E')dE' 


O 
(239) 


S(E) is now the spatially averaged slowing source of neutrons 
which appear in dE at E due to scattering collisions at 


energies greater than Ea: 


3(E) “2 S(Z,E)dv . (2.10) 
V 


As it stands, Eq. (2.9) is correct for an infinite 
reactor lattice and it could be modified to include the 
effects of thermal leakage out of the reactor assembly. 
However, since one is usually interested in the intensive 


properties of the lattice, conventional approaches neglect 








Sec. 2.1 imi 


this thermal leakage calculation. In addition, Eq. (2.9) 
becomes approximate as soon as one postulates appropriate 
mathematical models for the scattering kernel, =(E'—E), 
and for the relative fuel element absorption, £(E)/[1-£(E) ], 
in the quest for a practical and economic computer solution. 
The most general formulation for the relative fuel 
element absorption may be obtained by considering a unit 
cell in an infinite reactor lattice. Assume a source of 
neutrons in dE at energy E and position r, N(r,E), so that 
the spatially average source of neutrons, N(E), per unit 


volume, V, is expressed by: 


N(E) -¢ | sena (2.11) 


Given N(r,E), one can in principle calculate the net current 
into the fuel element, JE); at each energy. The total 
number of neutrons absorbed by the fuel element at energy 
E is then AJ (E) ; where A is the surface area of the fuel 
element. Since the total number of neutrons produced at E 
must equal the total number absorbed, the number of neutrons 
absorbed by the moderator is VN(E)-AJ, (E). The relative 
fuel element absorption is then expressed as: 

1s ia = AJ, (E) (2.12) 
-f(E) ~ UN(E)-AJ_(E) ° | 
In the assumption of a spatially dependent neutron 


source, N(¥,E), and the subsequent calculation of the net 








Secs. 2.1 - 2.2 17 


current into the fuel element, any of several means may be 
utilized including the use of multigroup transport codes, 
integral transport theory, or energy dependent diffusion 
theory. The accuracy of the results is directly related 

to the effort expended on these calculations, and this fact 
must be reconciled with the desire to produce economy in 
the use of the computer along with the knowledge that 
certain errors will be introduced with the use of mathe- 
matical models for the scattering kernel. 

It is to be noted that by incorporating and concen- 
trating all of the heterogeneous effects into the function 
f(E), the calculational simplicity of the infinite, 
homogeneous case has been retained as is seen in Eq. (2.9). 
A discussion of practical approximations for f(E) is found 


mmoees. (22.2) and (2.3). 


2.2 A Practical Approximation of Heterogeneous Effects 


In formulating an analytic approximation for f(E), the 
function which takes into account all of the heterogeneous 
effects in the lattice cell, one is forced to compromise 
between accuracy and computational ease. As an initial 
approach, the following simple model is postulated in an 
attempt to investigate the feasibility of this method in 
determining the thermal spectra within a lattice cell. In 
this work we have restricted ourselves to systems with 
cylindrical symmetry, recognizing that analogous results 


for spheres and slabs are obtainable. 
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Consider a cylindrical unit cell situated in the 
regular, rectangular or hexagonal, array of a lattice. 
It is characterized by the fuel element radius, ae and 
the moderator radius, r,- In determining ry from the actual 
lattice parameters, the Wigner-Seitz eriLemionee is applied. 
A spatially flat source of neutrons in the moderator, N(E), 
is assumed. The validity of this assumption for the 
calculation of integral properties has been discussed by 


21/ Les lie22/ has indicated that multigroup 


Pomraning. 
transport calculations yield moderator flux shapes which 
compare favorably with those predicted by one group 
diffusion theory. In light of this it is proposed to 
calculate the net current into the fuel element, J_(E); 
through the extension of monoenergetic diffusion theory to 
the polyenergetic case. The flux shape in the cell 
moderator at any energy E is assumed to be spatially dis- 
tributed in a manner given by solution of the constant 


cross section diffusion theory approximation, using cross 


sections evaluated at the energy E: 


D_(E)V“9" (rx, E) -22(E)o™ (x, E) ida) = (2.13) 


20/ S. Glasstone and M. C. Edlund, The Elements of Nuclear 


Reactor Theory, p. 265, D. Van Nostrand Company, Inc., 
Princeton, Toronto, London, and New York, 1952. 


21/ G. C. Pomraning, "On the Validity of the Constant Source 
Assumption for the Cell Problem," Nucl. Sci. Eng. 18, 
400-403 (1964). 


22/ D. C. Leslie, see Footnote 19. 
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The boundary conditions applied are 


Vo" (x7 ,E) = 0, (2.14) 
and 


m 
p (r,>E) 


Vo (r,,E) 


D_(E) = 1/3>,_.(E) is the moderator diffusion coefficient, 

5, (E) is the moderator macroscopic absorption cross section, 
and d(E) is the energy dependent linear extrapolation length. 
The linear extrapolation length may be determined from d.(E), 
the linear extrapolation length into a black cylindrical 
hole, 23/ by characterizing the fuel element with its escape 
probability, P(E). This has the added advantage of 
eliminating the need to state anything explicit about the 
shape of the flux within the fuel element. The fuel element 
Spatial flux is implied in any calculation of the escape 
probability, but by dealing strictly with P(E) one can 
extract and utilize in a simple manner the accuracy of 
transport theory calculations .24/ 
Considering the escape probability of the fuel element 


25/ 


as its albedo, and including the term in V6 which can 


23/ A. M. Weinberg and E. P. Wigner, The Physical Theory of 


Neutron Chain Reactors , p. 765, The University of 
Chicago Press, Chicago, 1958. 
24) K. B. Cady and M. Clark, Jr., "Neutron Transport in 
Cylindrical Rods," Nucl. Sci. Eng., 18, 491-507 (1964). 


25/ S. Glasstone and M. C. Edlund, p. 129, see Footnote 20. 
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AS 


be obtained in the derivation of Fick's Law,—’ one can write 


o™(r,,E)-2D,(E)Vo™(r,,E) + G(E)V*9"(r, ,E) 


P(E) = a 2.16 
‘ p) (r,,E)+2D_(E)Vo"(x,,E) + G(E)V“9"(r,,E) ——. 


For a black cylindrical hole, MS) is by definition zero so 
that G(E) is given by 


Vo (r,,E) 


aga) = | 20,(8)-a,(E) | (2.17) 


2m 
V'o (%,;E) 
In the actual cell, where P(E) is known from transport theory 


calculations, substitution of Eq. (2.17) into Eq. (2.16) 


yields with rearrangement 


o'(r),E) —  4D,(E)P,(E) 


d(E) = =< : 
@ (r);E) saa) 


+ di (E) . (2.18) 
An approximate expression for d,(E) is given in Appendix C 
while the calculation of the fuel element escape probability 
is discussed in Chap. 4. 

The relative fuel absorption, f£(E)/[1-f£(E)],may now be 
found from the solution of Eq. (2.13) utilizing the boundary 
conditions in Eqs. (2.14), (2.15), and (2.18). Fora 


cylindrical cell, Eq. (2.12) may be expressed as 


2 
AE) rr eG (2.19) 
2D_(E)Vo" (x, ,E) 


where y = r4/r,- Denoting 5, (E)/D_(E) by Ke, it follows 


that 


26/ A. M. Weinberg and E. P. Wigner, pp. 192-194, see 


Footnote 23. 
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- N(E)/,(E) 
Vo (r,,E) = Sy C2-20)) 
V(E)/« + d,(E) + TPE) 
where 
bez) K, (xr, ) 1, (er,) + I, (xr )K (xr,) ete 


Ky (xr, ) I, («r,) - I, («ry )K,(«r,) 


I I,; Ko? and Ky are the modified Bessel functions of the 


‘are 
First and second kind, respectively. Substitution into 
Eq. (2.19) yields 

L-£(E) 2% ol ¥*-L) SR (E) P(E) 


Es 
2 
Fe) ~— _I-F,ce). ‘+ 7 (y"-1) ex (E) 


2 m 
oy 1) 2 (aa) 
+ DCE = - 1. G27 


Insertion of Eq. (2.22) into Eq. (2.9) thus attempts to 
account for the heterogeneous effects in the lattice through 
use of diffusion theory in the moderator and transport theory 
in the fuel element. The formulation is not restricted to 
the use of diffusion theory in the moderator region but it 
is utilized because of its calculational ease and good approxi- 


mation to one velocity transport theory. 





=, 
= 
io 
7 2 a eae 
= > ee 
a —[—— ee ll 
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2.3 Formulation In Terms of Escape Probabilities 


The expression, Eq. (2.22), which approximates the 
relative fuel absorption, and hence, the heterogeneous effects 
in the cell may also be arrived at by the formulation of the 
energy dependent absorption rates in terms of moderator and 
fuel element escape probabilities. This approach is an 
extension of Amouyal, Benoist, and Horowitz's one velocity 
thermal utilization theory2_! to the polyenergetic ESA. 

One begins by defining the following energy dependent 
escape probabilities, all of which are calculated assuming 
the adjacent region is black: 

8) is the escape probability of neutrons incident 
on the fuel element in the actual lattice directional dis- 
eripDution, 

P.(E) is the escape probability of neutrons in the 
moderator at energy E due to collisions at other energies, 
and, 

P(E) is the escape probability of neutrons incident 
on the moderator from the fuel element. 

These escape probabilities are schematically depicted 
in Fig. 2.1 for a Wigner-Seitz cylindrical cell. If one 
27/ A. Amouyal, P. Benoist, and J. Horowitz, ‘''Nouvelle 

Méthode de Determination du Facteur d'Utilization 


Thermique d'une Cellule,'' J. Nucl. Energy 6, 79 
(1957). 
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Figure 2.1 The Fuel Element and Moderator Escape Prob- 
abilities 
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considers a spatial source of neutrons from scattering which 


has been normalized to one: 


[seme al (2.23) 


then reference to Fig. 2.1 shows that the current into the 


fuel element, J_(E), may be expressed as 
AJ_(E) = P,(E) + PE) Bey ee : C2224) 


oie. 
P(E) 


AJ_(E) = -——————__ 
at oe OD Pa OD) 


(2.25) 


where A is the surface area of the fuel element. The net 


current into the fuel element is given by 
J (CE) = J (E) - P(E) J_(E) ‘ (2.26) 


so that 


[1-P,(E) ]P, (E) 


AJ, (E) = 
1-P(E)P_(E) 


2 2) 


Since the source has been normalized to one neutron, Eq. (2.27) 
is the ratio of the absorption rate in the fuel element to the 
absorption rate in the cell. Therefore, 


[1-P,(E) ]P, (E) 
£(E) a “I-P (E)P_(E)_ ° (2.28) 
e m 





Be = eae a 


o a 
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Rearranging Eq. (2.28) or subsStitwtine  bqym(2-2)) eoneomtaE 
(2.12) yields: 


‘LAI P(E) | 1-P_(E) 1-P_(E) 
= + (2.29) 
E T-P,(E) | P, CE) “PL CE) 


Writing the heterogeneous effects in terms of these 
escape probabilities has several advantages. In the first 
place a certain amount of physical transparency is gained by 
separating the contributions of the fuel and moderator into 
easily defined quantities. Secondly, P(E) is a functional 
of the fuel element properties and its calculation depends 
only slightly on the moderator through the directional dis- 
tribution of the current entering the fuel elements. Bae 
and P(E) are functionals of the moderator and are very 
nearly independent of the fuel element properties. In 
addition, there are no restrictions on this simple model 
and the escape probabilities may be calculated with as much 
sophistication as is desirable. Calculational techniques 
for determining the moderator escape probabilities, P,(E) 
and P(E), are discussed in Chap. 3, while the fuel element 


escape probability, P(E), is discussed in Chap. 4. 








ia = > | a 





CHAPTER 3 
MODERATOR ESCAPE PROBABILITIES 


Sep ee) 

In this chapter the moderator escape probabilities, 
P(E) and PSE); as defined in Sec. 2.3 are discussed. 

It will be shown that by making appropriate assumptions, 
the formulation of the heterogeneous effects in terms of 
escape probabilities, Eq. (2.29), is identical to the 
expression derived utilizing diffusion theory in the 
moderator and characterizing the fuel element with its 
escape probability, P(E); as given in Eq. (2.22). 

The simplest model for 59 Bi) in cylindrical geometry 
is obtained by assuming, as in Sec. 2.2, a spatially flat 
source of neutrons in the moderator. Again as an expediency 
it is assumed that the flux shape in the moderator is 
given by one velocity diffusion theory applied at each 
energy. One finds P,(E) by calculating the net current 
into a black cylindrical hole of radius r,. For a spatial 


source normalized to one neutron, P(E) is given by 


P.(E) = AD_(E)Vq" (r,,E) (3.1) 


Solving the diffusion equation, Eq. (2.13), with the 


boundary conditions given by 


Von(r,,E) = 0 , (3.2) 


Se 
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and, 
9” (x,,E) 
= aE) 
Vo (r_,E) 
O 
yields 


1/(veze(E) | 


Von (r5,—E) = ————_—— 
° X(E)/« + do (E) 


(323) 


v" is the moderator volume, Ke 


- sl : 
= = (E) /D_,(E) and y(E) is 
given by Eq. (2.21). In accordance with Eq. (3.1), Bay 
is then given by 


2 1 
EAE) = — >... 1°) == (3.4) 
; yg) X(E) + Kd. (E) 


Or rearranging, 
1-P_(E) v6 1g 
Ey (y*-1)K7d (E) + 2 (y?-1)ex(E) - 1 
(3.5) 


As before, y = r,/r,- Figure 3.1 shows P,(E) as calculated 
for the Cornell University Zero Power Reactor cores. The 
effect of various water to uranium dioxide volume ratios 

is shown as well as the energy dependence. The parameters 
for the cores utilized in these calculations are given in 
Appendix A, while the cross section data incorporated in 


these computations are discussed in Chap. 7. 
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Br2 Pky 


The moderator escape probability for neutrons incident 
on the moderator from the fuel element, ge Oa) may be cal- 
culated from 5) through a reciprocity relationship under 
the following assumptions. If the flux in the moderator is 
isotropic, and if ray is calculated using a spatially 


uniform source, as it has been, then P.(E) and P_(E) are 


related by28:29/ 
1-P_(E) m 
Ey = 424 (E) a (3.6) 
S 


where for cylindrical cells 


m r 
Y= 2(y?-1) . (3.7) 
Ee CE ) is shown in Fig. 3.2 as calculated from Eq. (3.6) 
utilizing Way) as defined by Eq. (3.6) and depicted in 
mo, 3.1. 
Insertion of Eqs. (3.5), (3.6), (3.7) intonhaqenc2-2» 


yields an expression for the relative fuel absorption: 


1e£(E) 2% o(¥*-1) a (E)P, (E) 


iG 
je ie) Hy*-1) ex (E) 
e 


ro (y2-1)2"(E)d, (E) 


+) === at (Ep ee a (3.8) 
™m 


28/ K. M. Case, F. de Hoffmann, and G. Placzek, Introduction 


to the Theory of Neutron Diffusion, Vol. I, United 
States Government Printing Office, 1953. 


29/ K. B. Cady, "Neutron Transport in Cylindrical Rods," 
Massachusetts Institute of Technology Thesis, 1962. 








Zz 


ws 


Sec. 


O° T 


W 
(a) d ‘AIT[Tqeqorg edeosy A0ReAEepoW SUL z‘¢ 9INST A 
(A®) q ‘AB20Uq 
ay TO" 100° 


OTIeY sUNTOA 
OPTxXQ/1799eM 


Jo0¢ 
S910) [Teuro9 z 


G*0 


2 0 


L°O 


8° 0 


6°0 


O°, 


(a) d 





Sec. 3.2 26 


which is identical with Eq. (2.22) as derived in Sec. 2.2. 
This is the expression which will be incorporated in the 
calculations discussed in Chap. 7 in an attempt to account 
for the heterogeneous effects of the lattice. Due to the 
physical transparency and separability of effects achieved 
through the consideration of the various escape probabilities, 


all future discussion will be in terms of them. 








CHAPTER 4 
THE FUEL ELEMENT ESCAPE PROBABILITY 


4.1 The Fuel Rod Escape Probability 

As was noted in Sec. 2.2, it is in many respects desir- 
able to characterize the fuel element in terms of its 
escape probability, P(E), thereby eliminating the need 
to explicitly define the flux in the various regions of the 
element. In addition, computation of Ee (Ey allows one to 
complete the calculation of the relative fuel absorption as 
defined in Eq. (3.8). 

In the most general case, the fuel element is composed 
of a cladding region, a void region, and a fuel rod or fuel 
meat region. In the determination of P(E); the fuel rod 
may be characterized by its escape probability, P(E). The 
effect of the void may be accounted for by strictly geometric 
factors, while the cladding effects may be expressed in 
terms of several transmission factors. The results derived 
from the calculation of any of the above factors are of 
course dependent upon the assumption of an angular distri- 
bution of the current entering each respective region. 

The fuel rod escape probability depends on the fuel rod 
radius, a = St (E) rg, measured in mean free paths (mfp's) and 
the ratio, c = 5< (E) /2; (E) of the fuel scattering cross 
section to total cross section. The simplest expression 
for P(E) is derived from diffusion theory and it is just 
the albedo of the rod: 


E972 
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_ 1-2D¢(E)B(E) 


ae ooo (4.1) 
1+2D,(E)B(E) 
where 
ee Cera. 
Le Sate te 
BCE) =k (42) 
£ - Ket ¢ 


Ly and I, are the modified Bessel functions of the first kind, 
while x,* = SI (E)/D-(E), Dg(E) being the diffusion coefficient 
omeene fuel. The quantityxp may be obtained either from Pa! 


theorys2/ 


k¢/54(E) = Walesa) . (4.3) 


or from Wigner's "elementary theory," 


¢ tanh" (q/EE(E)) = /EE(E) - (4.4) 


Transport theory expressions for the fuel rod escape 
probability have also been developed from Boltzmann's one 
velocity transport equation under the assumption of 
isotropic scattering in the laboratory system. This results 


in the integral transport or Peierl's equation: 
aoa 


=F 


p(r') 


: | (4.5) 


(x) = 9, Cr) + on, | 


4n|r-r' 


30/ A. M. Weinberg and E. P. Wigner, The Physical Theory of 


Neutron Chain Reactors , The University of Chicago 
Press, Chicago, 1958. 
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where 9, (¥) is the uncollided flux due to incident surface 


sources. Escape probabilities in terms of the solution of 


Eq. (4.5) have been considered by several authors .24232233/ 


This work utilizes the treatment by Cady for an isotropic 


flux incident upon the fuel rod. The calculational techni- 


ques utilized in the computation of if (C20) are quoted eysikevas al! 


For thin rods of less than 0.5 mfp's a very accurate 


expression results from the uniform flux approximation. In 


ce 


nee this yields 


terms of Placzek's collision probability, 


the following expression for P(E): 


1-M 
P(E) = 1-2a(1-c) rates | (4.6) 
OO 


ae is calculated in terms of the modified Bessel functions 


of the first and second kind with an argument of a: 


_ 2 : 
Moo = L > F | 1K, +1-2al 1K +1-a(1,K,+1,k,)] ) 





(4.7) 
For rods which have radii, a, greater than 0.5 mfp's, 
the parabolic flux approximation yields more accurate results. 


In this case 


2 
M__-c(M__M,,-3M_“) 
P(E) = 1-2a(1-c) =) et 


(4.8) 


where the following quantities are calculated in order: 


f= 
_ 


/ K. B. Cady, "Neutron Transport in Cylindrical Rods," 

Massachusetts Institute of Technology Thesis, 1962. 
32/ K. B. Cady and M. Clark, Jr., "Neutron Transport in Cylin- 
drical Rods,"' Nucl. Sci. Eng. 18, 491-507 (1964). 
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A = a[ I, (a)K, (a) ae (a) K (a) J 5 (4.9) 

B = 2a(14K, le) (4.10) 

C = 3(1,K, + 1-38), (4.11) 

D= (I,K, -1+B+C), (4.12) 

B= e(—2t+c-p), (4.13) 

a rz 2 E) , (4.14) 

Whee ES Ges (4.15) 

i) = 1-2C+E, (4.16) 
and, 

Mon7 = 1- 30+ 3E -F . 

Consideration of the asymptotic expansions for the moments, 

ee? Moo and My9> yields the following approximations for the 


fuel rod escape probability. For the uniform flux approximation, 


P(E) ¥———+~—__ =» —___+—__ , (4.17) 
1 + 2a(1-c) 1+ 2f(E)r, 


and for the parabolic flux approximation, 
lt4a(1-c) + = aW(1- -c)? + 32 a(1- ya) 
(4.18) 


33/ G. W. Stuart and R. W. Woodruff, ''Method of Successive 


Generations," Nucl. Sci. Eng. 10, 388 (1961). 
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Figure 4.1 shows the fuel rod escape probability as cal- 
culated utilizing the expressions stated above. The effect 


of the 0.27 2? 


resonance is plainly evident. The parameters 
used in the computations are those of the low enrichment, 

high density uranium dioxide fuel of the Cornell University 
Zero Power Reactor. These parameters are listed in Appendix 
A, while cross section data are discussed in Chap. 7. The 


numerical methods needed for calculation of the modified 


Bessel functions are outlined in Appendix C. 


4.2 Void Considerations 


Any void space which is present between the fuel rod 
and the cladding is accounted for by considering the fraction 
of neutrons entering the void space from the cladding which 
are heading toward the fuel rod. The escape probability for 


the fuel rod and the void, 1B gE): may then be expressed as 


P(E) = (1-6) + ¢P.(£), (4.19) 


where ¢ is the fraction of the neutrons headed toward the 


fuel rod. For an assumed isotropic flux emergent from the 


cladding, 36/ 

e = re/t. (4.20) 

34/ K. B. Cady, see Footnote 3l. 

33/ K. M. Case, F. de Hoffman, and G. Placzek, Introduction 
to the Theory of Neutron Diffusion, Vol. I, United 
States Government Printing Office, 1953. 

Sey, 


K. B. Cady, p. 161, see Footnote 31. 
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Figure 4.1 The Fuel Rod Escape Probability, P_(E) 








Bees. 4.2 - 4.3 58 


where Ie is the fuel rod radius and r is the inner radius 
of the cladding. Examination of Eqs. (4.19) and (4.20) shows 
that the correct limits exist in that Eaeaeey) goes to Pe 


as r, goes to rg, and eee (OE), goes to one as r¢ goes to zero. 


4.3 Cladding Effects 


The fuel element cladding is treated in accordance 
with the cylindrical thin region theory developed by Cady 
The neutron reflection and transmission characteristics of 
the cladding can be divided into several fractions. For 
neutrons incident on the outer surface of the cladding, 

l. g is the fraction of neutrons, which return 

uncollided to the moderator, 

2. T is the fraction transmitted uncollided through 

the cladding, 

3, ores is the fraction scattered in and transmitted 

through the cladding, 

4. eae is:the fraction scattered in and reflected 

out of the cladding back to the moderator, and 


a Der URE echt oe) is the remaining fraction which 


is absorbed in the cladding. 


Similarly, for neutrons incident on the cladding from the 


fuel rod region, 


37/ K. B. Cady, p. 63, see Footnote 3l. 





” 
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1. t is the fraction transmitted uncollided through 


the cladding into the moderator, 


Ze P;,, is the fraction scattered in the cladding and 
reflected back into the fuel rod and void region, 
3. is the fraction scattered in and transmitted 


Pout 
through the cladding, and, 


4. 1-(t+p;,,+9,,,,) is the remaining fraction which is 


absorbed in the cladding. 


In terms of these energy dependent fractions, the 
escape probability for the entire fuel element may be 


expressed as 
a) cee ey aa (T+) op) I, (ED , (4.21) 


J, (E) is the net current entering the inner surface of the 
cladding from the fuel rod and is a function of lea (CO) 
the escape probability of the void and fuel rod region: 
J, (E) = a (4.22) 
hoe an 

The computational receipes used for the calculation 
of the various fractions mentioned above are stated below. 
These schemes assume isotropic fluxes on both sides of 
the cladding and they are performed at each energy where 
applicable. The outer radius of the cladding is denoted 
by E53 the inner cladding radius is ro: Macroscopic cross 


sections refer to the cladding material. 
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t= 1-EC(E)r, Sy + F(SC(E) 4)? Be (4.40) 
y = 0.577215 ..., Euler's Constant, (4.41) 
as ye (4.42) 
x 3 x2 
P,=% (y+ In(x) -5] +1-% 
x ae eo 
* 28 - 360 + 2880 » 08) 
and, 
u (4.44) 
ee ———EE ia 
unif =¢(E) 
- — (1-P,) 
ze (E) = 


Hoe x < .05 , 


g = (1-F,) - 22C(E)r (1-Fy) + 3(2E(E)r,)* + AP (1-F3) 





(4.45) 
mime for .05<x < 0.2, 
g = -.0212 - .0266 In(x)/@2°? . (4.45) 
C 
S 
A as Cre (US ee i) beats ; (4.46) 
2, (E) 
z. (E) 
Pin T Pout = Te7ny (1-7) Punif ? (4.47) 


E,(E) 


f= + sin™* 25 | (4.48) 





-_ 
-_ 
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Tin > S eaeeien vy tt 2 (4.49) 

Bear = (hgeece eae eae ee (4.50) 

Pin = Fin (in + Pout? » (4.51) 
and, 

O) A epee Uelrievens, [OP ery Miata (oir. (4.52) 


Calculations have been made using the above techniques 
for the fuel elements of the Cornell University Zero Power 
Reactor, the specifications of which are given in Appendix A. 
The results of these computations in the form of the fuel 
element escape probability are shown in Fig. 4.2. In these 
calculations, the parabolic flux approximation, Eq. (4.8), 
was used for the calculation of P(E). Again the effect of 


DE 
. resonances can be seen. 


the U 
In terms of previously defined quantities, the ratio 
of the number of neutrons absorbed in the fuel rod to the 


number absorbed in the fuel element is given by 


rod abs. rate _ J, (E) (1-P) (CE) J 


a (4.53) 
element abs. rate ee, Jey) | ke (08 | 


Similarly the ratio of the number of neutrons absorbed in 
the cladding to the number absorbed in the fuel element 


can be defined by 
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eladeabsnm Gaten se J, (E){1-P_(E)) 


; (4254) 
element abs. rate Bey 2) eee 
B, (E) 
= ae, Ces ey) ee 
S(E)tEZS(E)P,, 
(4.55) 


The results of performing the calculation expressed by 
Eq. (4.53) are shown in Fig. 4.3 for the Cornell fuel 
elements. Inspection of Fig. 4.3 shows that only about 1% 
of the neutrons absorbed by the fuel element are absorbed 
in the cladding. This is because of the low absorption 
cross section of the aluminum cladding. On the other hand, 
in the Brookhaven uranium dioxide fuel elementenee (see 
Chap. 8), where stainless steel is used as the cladding 


material, about 12% of the neutrons absorbed in the fuel 


element are lost to the cladding. 


38/ H. C. Honeck, "The Calculation of Thermal Utilization 


and Disadvantage Factor in Uranium Water Lattices," 
Nucl. Sci. Eng. 18, 49-68 (1964). 
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CHAPTER 5 
MODELS FOR THE SCATTERING KERNEL 


5.1 The Monatomic Gas Model 

By formulating the lattice heterogeneous effects in 
terms of contrived absorption rates as in Secs. 2.1 and 2.3, 
one is left free to incorporate as much sophistication as 
is desired into the calculation of the energy transfer or 
scattering kernel to be used in Eq. (2.9). With this type 
of an approach one can determine directly the effect of 
different mathematical models on the neutron spectra and 
lattice integral properties. One can then decide on the 
basis of computer time used and the accuracy of the results 
whether or not use of the more complicated models for the 
kernel is desirable for a particular computation. 

Numerous mathematical models for the scattering 
kernel have been proposed 222 tO eee In an attempt 
to bracket the degrees of computational difficulty and the 
accuracy of the results predicted with these many models, 
we have limited ourselves to calculations with only two 


of them; the Wigner-Wilkins monatomic gas modete and the 


Nelkin water es 


39/ E. P. Wigner and J. E. Wilkins, Jr., "Effect of the 
Temperature of the Moderator on the Velocity Distri- 
bution of Neutrons with Numerical Calculations for H 
as the Moderator," AECD-2275 (1944). 

40/ 


J. E. Wilkins, Jr., "Effect of Temperature of the 
Moderator on the Velocity Distribution of Neutrons 
for a Heavy Moderator,' CP-2481 (1944). 


sy Ae 
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The choice of the monatomic gas model is mainly 
dictated by its ease of calculation. This model treats 
the moderator as a free gas in thermal equilibrium at 
the temperature of the moderator. Although crystalline 
or molecular binding effects in the thermalization prop- 
erties of the moderator are neglected, infinite homogeneous 
media peeee dene aren this model have yielded spectra 
which are essentially consistent with those measured 
experimentally. 

The cross section for scattering from an energy E' 
to an energy E through an angle with a cosine of wu may be 
derived from the short collision time approx aera 
applied to the single oscillational mode of the trans- 
lational motion of the scattering particle. For the 


situation where E' is greater than E, this results int?/ 





O Lif 7 2 
o(E'—E,u) = 2 (= | exp | ty (E-E' + 5) | ’ 


(5,1) 


ib 


L/ M. S. Nelkin, "Scattering of Slow Neutrons by Water," 
Phys. Rev. 119, 741-746 (1960). 


IS 


/ A. Radkowsky, "Temperature Dependence of the Thermal 
Transport Mean Free Path,'' Physics Quarterly Report, 
(April, May, June, 1950), ANL-4476, p. 89, Argonne 
National Laboratory (July, 1950). 

eS! H. D. Brown and D. S. St. John, "Neutron Energy 

Spectrum in D,0," DP-33 (1954). 
gay The Scattering Law, Proceedings of the Brookhaven Con- 
ference on Neutron Thermalization, BNL 719, Vol. I, 
Brookhaven National Laboratory (1962). 
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where o(E'—~E,u) is in barns/(eV-unit cosine). The bound 
atom cross section in barns is Op» T is the moderator 
temperature in electron volts, and A is the ratio of the scat- 
tering nucleus mass to the neutron mass. The quantity X 


is related to the momentum transfer in that 


X = E+E’ - 2uVvEE' . (522) 


A similar result can be obtained for upscattering in energy 
by employing the principle of detailed balance, Eq. (2.5) 
as is discussed below. 

The scattering kernel is found by the angular inte- 


gration of Eq. (5.1), 
1 


o(E'—>E) -| o(E'-E,u)du . (5.3) 
-1 


This operation can be carried out analytically and it results 


in the formulation of the Wilkins' heavy gas node for 


the scattering kernel. For down scattering in energy, 


45/ E. P. Wigner and J. E. Wilkins, Jr., see Footnotes 39 


and 4O. 


46 / M. S. Nelkin, see Footnote 41. 


47/ L. G. de Sobrino and M. Clark, Jr., "A Study of the 
Wilkins Equation,” Nucl. Sez. Engs 10.5385 (on 
"Comparison of the Wilkins Equation with Experiments 
on Water Systems," Nucl. Sci. Eng. 10, 377 (1961). 


M. S. Nelkin, see Footnote 41. 


49/ F. D. Federighi and D. T. Goldman, "KERNEL and PAM - 
Programs for Use in the Calculation of the Thermal 
Scattering Matrix for Chemically Bound Systems," 
KAPL-2225, Knolls Atomic Power Laboratory (1962). 


50/ J. E. Wilkins, Jr., see Footnote 40. 
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E' > E, the monatomic gas scattering kernel is expressed in 
barns /eV asot/ 
, a gn 
o(E'—E) = rs ERF(nz-pz') + ERF( z+pz' ) 
2E 


+ exp(z!*-z*)(ERF(nz'-pz)-ERF(pztnz')] |, (5.4) 
where 

zo = E/T , (5.5) 

git = E'/T , Geo 

n = Hatta ll? (5.7) 
and, 

p = H(A-L)A 1/2 : (G8) 


ERF(z) is the error function defined by 
Z 


ERF(z) = 2 Re dy (5.9) 
Jr 3 


O 
while Of is the free atom scattering cross section for the 


moderating medium. 


(5.10) 





A 
A+1 





=e ap 
The scattering kernel for upscattering, E' < E, is cal- 
culated utilizing detailed balance as expressed in Eq. (2.5). 


In the case where one is considering moderation by 


hydrogen, Eq. (5.4) reduces to the Wigner-Wilkins model for 
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a gas of mass one .22/ Then” fore Sey 
CF 
o(E'—=E) = aT ERF(z) . (5.11) 


The use of the monatomic gas model has another advantage 


in that the scattering cross section defined by 


co 
0, (E) .| O(E=E')dE' , Corl») 
O 
is an analytic function given by 
Z 
o (E) =o (1 + l ) ERF(z) + eo (55153) 
¢ —7 Z, ; 
2y J Ty 
where 
yo = BA, (5.14) 


This analytic function for the scattering cross section is 
useful because it allows one to check on the accuracy of 
the numerical integration scheme employed. This aspect is 
more fully discussed in Sec. 5.2. 

The calculations performed utilizing the monatomic 
gas model for the scattering kernel are discussed in 
Chap. 7. Calculation of the error function necessary in 


the computations of this kernel is detailed in Appendix C. 


1/ M. J. Poole, M. S. Nelkin, and R. S. Stone, ‘'The Measure- 
ment and Theory of Reactor Spectra,'' Progress in 
Nuclear Energy, Series I, Vol. 2, pp. 9l- , Pergamon 
Press, New York, London, Paris, and Los Angeles, 1958. 


32/ E. P. Wigner and J. E. Wilkins, Jr., see Footnote 39. 
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5.2 The Nelkin Water Kernel 


The second mathematical model for the scattering kernel 
utilized in these calculations is that proposed by Nel kin22/ 
for protons bound in water. Although the computer cal- 
culation of the Nelkin water kernel takes a relatively 
long time, its use is justified in that infinite homogeneous 
media thermal spectra calculated with this model are in good 


34 / In addition, 


123/ 


agreement with those measured experimentally. 
diffusion coefficients calculated with the Nelkin mode 
agree well with measured sealers 

Physically, the mathematical approximation proposed by 
Nelkin for water treats the scattering molecule translation 
with the free gas model utilizing a translator mass (M) of 
18. The effects of the atomic motions within the molecule 
are approximated by simple harmonic oscillations, the 
phonon frequencies of which determine the relative strength 
of the interaction with a particular oscillator. The 
oscillators treated are as follows: A hindered rotational 
oscillator of mass (M) 2.32 and frequency (CO,) corresponding 
to an energy (w, = fiQ,) of 0.06 eV, where fi is Planck's 


Constant divided by 27; a vibrational oscillator of 


mass (M,) 5.84 and an energy (1) of 0.205 eV; and a 


53/ 
54/ 


M. S. Nelkin, see Footnote 41. 


J. R. Beyster, et al., "Measurement of Neutron Spectra 
in Water, Polyethylene, and Zirconium Hydride," Nucl. 
Selene. 9, 168 (1961). 
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second vibrational oscillator of mass (M5 = M2) 2.92 and 
energy (D5) of 0.481 eV. Different mathematical approxi- 
mations may be made to treat the translational, vibrational 
and rotational modes, dependent upon the incident neutron 
energy and the energy or momentum transfer. However, a 

single expression may be used to characterize the differential 
scattering kernel for scattering from an energy E' to an 
energy E (E' > E) through an angle with a cosine of nu. 
Parameters within the Nelkin water kernel then change 


dependent upon the specific approximation made, 


Timm 














°b E -< il Dye oe 
ou EF, U) = e » = 
m=0 
LL , (exe) nw ; 
Hm] Do ae ov [- ake ee 
£=0 N=-00 
+ X4 for 2 5.15 
Ti ars) =i 1 + nw ) : ( ‘ ) 


or) 


55/ G. P. Calame, "Diffusion Parameters of Water for Various 
Scattering Kernels,'' Nucl. Sci. Eng. 19, 189 (1964). 


Nn 


6/ E. Starr and J. Koppel, ‘'Determination of Diffusion 
Hardening in Water,'' Nucl. Sci. Eng. 14, 224 (1962). 
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Equation (5.15) gives the differential scattering kernel in 
barns/(eV-unit cosine) with Oo, as the bound atom cross 
section in barns. The I, (2) are the modified Bessel 
functions of the first kind and order n. As before, the 


quantity X is given by 
X = EtE'-2u VEE’ . (5.16) 


One may allow for all of the transitions for the 
translational and rotational modes, while utilizing a 
phonon expansion for the first (1 = .205 eV) vibrational 
mode and the elastic (unexcited) limit for the second 
(®, = .481 eV) vibrational mode by inserting the following 


values into Eq. (5.15): 








ie 1 , a 2 
— Mw M : 

w M_tanh(o, /2T) vee ve 
£ = T/M , 

Mwsinh(w/2T) 

(5.1 

tui 
Mi M 
mm= 0 , 
Lh=2, 
OS i Ce 
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and, 


where T is the moderator temperature in eV. The phonon 
expansion for the second vibrational mode is included if 
mm is set equal to 2 in the above set of parameters. 

To allow for all transitions of the translational 
and first vibrational modes, while treating the rotational 
mode with the short collision time approximation and the 
second vibrational mode with a phonon expansion, the 


following parameters are inserted into Eq. (5.15): 





1 al 2 
— = ——— + —— _, 
A wy De 
ay 

T r 1 1 
€-gte (mre t+ Zz]; 

M se er Ty 2 
Z = A 

Mosinh(w/2T) 

(5.18) 
Ll _é4t 1 
MOM” HM ° 
mm = 2 , 
Lh=0OQ, 
O = D4 
and 

M=M 
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The oxygen present in the water is accounted for by adding 
the monatomic gas model with mass 16, Eq. (5.1), to the 
Nelkin water kernel, Eq. (5.15). 

Angular integration of the differential scattering 
kernel, o(E'—>E,u), was performed numerically to obtain 
the scattering kernel in the form utilized in lattice cal- 
culations. For the non diagonal terms in the scattering 
matrix (E'+E), the Gauss quadrature formula was used for 


the angular integration: 


oC.) - | o(E'—~E,u)du 
-1 
kk 
av [OEE uy) + o(B'—=E, -14)), 


k=1 
(5.19) 


where w, is the weighting coefficient corresponding to My 


57/ For terms where 


a zero of the Legendre polynomials. 
the energy transfer is large (E << E'), an eight point 
(kk = 4) integration is used, while for small energy 
transfers where the fine structure of the differential 
kernel is a rapidly varying function of angle, a sixteen 
point (kk = 8) quadrature integration is performed. 
Calculation of the scattering kernel diagonal terms 
(E = E'), where there is no energy transfer, requires some 
a7! 4. N. Lowan, N. Davis, and A. Levinson, "Tables of the 
Zeros of the Legendre Polynomials of Order 1-16 and 


the Weight Coefficients for Gauss's Mechanical Quad- 
rature Formula,'' Bull. Am. Math. Soc. 48, 739 (1942). 
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modification of the above procedure. Inspection of Eq. 
(5.15) shows that for E' = E, X = 2E(1l-u), and there is 
a singularity in the differential scattering kernel at 
u=l. 

The use of the Gaussian quadrature angular integration 
scheme introduces numerical errors in this situation. 
Therefore the following calculational technique is utilized 
as suggested by Federighi and Coleen) Equation (5.15) 
is considered in the limit as X goes to zero (u goes to 1). 
The singularity is subtracted from the differential kernel 
and integrated analytically, the remainder being integrated 
with the sixteen point Gaussian quadrature rule. The 


analytically integrated singularity is then added to the 


result. This process is indicated below: 





1 
2 d 
o(E—~E) -| f(n)du + b if a (5.20) 
Te ; E(1l-u 


o 4 fre J 2E (1-1) 


where, 


°D 
f(u) = o(E—~E,y) - Gal) 


4 [TeX 


and X = 2E(l-uw). The diagonal terms are then approximated 





by 


. 38/ F. D. Federighi and D. T. Goldman, see Footnote 49. 








Sec. 5.2 2 


8 
(E—~E) e( +E ( ) ime 
O —— a = W - 

| ke | 2 fate 





C722) 

In the calculations of the Nelkin water kernel made for 
incorporation into the lattice studies Eqs. (5.17) were 
utilized for incident energies, E', less than 0.33 eV, 
while Eqs. (5.18) were used for E' greater than 0.33 eV. 
This is essentially the same approximation made by Honeck 
for the water kernel calculations incorporated in THERMOS .22/ 
The phonon expansion for the second vibrational mode along 
with Eqs. (5.17) is not used because only a small increase 
in accuracy at higher incident energies is accomplished 
with a large (by a factor of 3) increase in computer time 
through its utilization. 

Figure 5.1 shows the results of a calculation which 
Medea zed both Eqs. (5.17) and (5.18) in Eq. (5.15) for an 
incident energy of 0.366 eV. The results are expressed in 
barns/unit lethargy, Eo(E'—>E), for down scattering, 

E' > E. However, the basic difference between the two 
approximations is plainly evident in Fig. 5.1. By con- 
sidering all transitions for the hindered rotator, Eqs. 
(5.17), a scattering kernel is predicted with a complicated 


fine energy structure having oscillations at energy 


39/ H. C. Honeck, ''THERMOS - A Thermalization Transport 


Theory Code for Reactor Lattice Calculations," 
BNL 5826, Brookhaven National Laboratory (1961). 
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intervals (E'— E) equal to the frequency of the rotator, 
oa 0.06 eV. In going to the short collision time approxi- 
mation for the hindered rotator, Eqs. (5.18), this fine 
structure is considerably smoothed, however the cross 
section derived by the energy integration of the kernel, 
Eq. (5.12), is essentially the same for both models. The 
choice of 0.33 eV as the energy above which Eqs. (5.18) will 
be used and below which Eqs. (5.17) will be used is made 
with this fact in mind, but primarily on the basis of 
computer economy. In order to predict accurate cross 
sections utilizing all transitions for the hindered 
rotator, it is necessary to use a fine energy mesh in the 
calculations (AE < .06eV), but if one is interested in 
energies in the range up to 1.0 eV this means that an 
excessively large number of energy points must be used 
in the calculation. This problem is further discussed 
below and in Chap. /7. 

The first moment of the scattering kernel, O7 (== EB) 


is calculated for down scattering utilizing the same 


techniques as above, where, 
1 


CCE aE) -| wOo(E'——E,u)du . (5723) 
-1 
As before, a sixteen point Gaussian quadrature angular 
integration is made for the off diagonal terms when the 


energy transfer is small, while eight points is considered 
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sufficient for large energy transfers. The diagonal terms, 
E' = E, for the first moment are calculated using the 


function defined in Eq. (5.21): 
iL 


1 
O 
(E —>E) = £ Ale e eel di ; 


-1 





8 

£(n,)-u£ Gn) |e —— 
ow WwW = = 

L, ele eed 6 [TEE 


(5225) 

The upscattering terms for both the scattering kernel 
and its first moment are calculated using detailed balance, 
Eq. (2.5). It is to be noted that a slight, but 
acceptable inconsistency is introduced in this manner, 
since the use of detailed balance means that in the up- 
scattering portion of the scattering matrices, the 
approximation used, Eqs. (5.17) or Eqs. (5.18), actually 
depends on the neutron final energy rather than its initial 
energy. 

Energy integration was performed on the scattering kernel 
and its first moment to determine the average cosine of 
the scattering angle, u(E): 

E 


Cc 
if 0, (E—>E) dE’ 


= O 
u(E) = a ‘ (5.24) 


Cc 
ij o(E—>E) dE! 
O 








sec. 5.2 56 


A trapezoidal rule numerical integration scheme, Appendix C, 
is used for this calculation and the results are shown in 
Fig. 9.2. Figure 5.2 also shows the transport cross section 


of water as calculated using the expression 
o,(E) = 0, (E) = o(E PEt Cas (3325)) 


In the energy integration of the scattering kernel, 
Hq. (5.12), to obtain the scattering cross section, J, (E), 
certain numerical inaccuracies were encountered. Specifically, 
numerical integration of the Nelkin water utilizing the 
trapezoidal integration scheme discussed in Appendix C con- 
sistently yields scattering cross sections which were higher 
than those measured ee perimentaite 22) for energies greater 
than 0.1 eV. Essentially the same percentage error is seen 
when the numerically integrated monatomic gas kernel is com- 
pared with the analytic expression for the cross section 
given by Eq. (5.13). These numerical errors are attri- 
butable to the fact that at higher incident energies, the 
diagonal terms of the scattering matrix are large compared 
with the first off diagonal terms for the energy mesh used 
in the calculation as discussed in Chap. 7. Use of the 
trapezoidal integration rule in this region therefore yields 


scattering cross sections for both kernels which are higher 


60/ D. Hughes and B. Schwartz, "Neutron Cross Sections," 


BNL 325, 2nd ed., Brookhaven National Laboratory 
(1958). 
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than experimentally measured or analytically calculated 
values. Auxiliary calculations showed that good agreement 
between See) calculated with the Nelkin water kernal and 
measured experimentally or between oO, (E) calculated by the 
monatomic gas model and 0, (E) given by Eq. (5.13) could 

be obtained by either utilizing a much finer energy mesh 

or by using more sophisticated integration schemes. However, 
both of these approaches create large increases in the com- 
putational time since the kernel calculation time is pro- 


2 where n is the number of energy points con- 


portional ton 
sidered, and the more complicated integration schemes which 
involve the fitting of exponentials or higher order 
polynomials to the calculated kernel points require 
additional computer time due to the uneven energy spacing 
used in the calculations. It is therefore desirable to 
retain the simplicity of the trapezoidal rule and correct 
for its inability to handle accurately the large diagonal 
terms encountered at high energies. This is done by 
recognizing that the actual value of the diagonal terms, 
E' = E, has no effect on the spectrum calculated with 

Eq. (2.1) or Eq. (2.9). This is seen by substituting the 
definition of the scattering cross section, Eq. (2.2) 

into Eq. (2.1) or (2.9) and letting E go to E’. Further- 
more, integration of both sides of Eq. (2.1) or Eq. (2.9) 
over the energy variable E indicates the scattering cross 


section used in calculations must be the one derived from 
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the scattering kernel by the numerical integration scheme. 
Since some of the integral properties that are to be 
calculated are strongly dependent upon the use of an 
accurate scattering cross section, we propose to modify 
the calculated diagonal terms in such a way that the 
scattering cross section calculated with the trapezoidal 
numerical integration scheme will equal the measured values 
for the Nelkin water kernel and Eq. (5.13) for the monatomic 
gas model. Thus, the diagonal terms, o*(E->E), utilized 


in the calculations will be given by: 
o*(E—>E) = C(E)o(E—E) , (5.26) 


where o(E—~E) is calculated using Eq. (5.22). The correction 
factors, C(E), are easily found for the trapezoidal inte- 
gration rule from the difference between the measured or 
analytic scattering cross section and the scattering cross 
section using the diagonals, o(E—E). These correction 
factors are assumed to be temperature independent and those 
calculated for room temperature may therefore be applied to 
the diagonals, o(E-~E) calculated at any temperature. It 
is to be noted that this technique is not applied to the 
calculation of u(E), Eq. (5.24), since the most important 
effect is the relative magnitude of the correctly calculated 
values of CG (Eh FFE) and o(E'-—E). 

Appendix B gives a listing of the scattering matrix 


as calculated using the Nelkin water model. The physical 
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temperaturesof 20°C and 40°C are considered. 

The source of thermal neutrons, S(E), due to down 
scattering from energies above the thermal cutoff energy, Ea 
can be calculated in accordance with Eq. (2.3) by taking the 
short collision time approximation on all of the modes of 
oscillation in the Nelkin model. Assuming a 1/E flux above 


E. and neglecting upscattering, Eq. (2.3) yields an un- 


normalized source of: 


S(E) = ERF E/T, F (5.27) 


where ERF(z) is defined in Eq. (5.9). 


The effective temperature, Wee is given by: 


E E 2E 
il Se ; (5.28) 


: : &R 


(5.29) 





Because of the relative magnitudes of , and 5), T, may be 


approximated by: 


1 + 1 
eor/T_y 2) 





Ww a) 
er ae 
Te “Mt OU 7 My, 
(5.30) 


At a physical temperature of 20°C the effective temperature 


has a value of 0.1169 eV. 





He => 





CHAPTER 6 
THERMAL LATTICE SPECTRA AND INTEGRAL PROPERTIES 


Having selected appropriate mathematical models for the 
scattering kernel, Eq. (5.4) or Eqs. (5.15) and (5.19), and 
for the relative fuel absorption, f(E)/[1-f£(E)], Eqs. (3.8) 
and (4.21), one proceeds to look for solutions to Eq. (2.9), 


which is rearranged below. 


E 
Cc 
S(E) | =(E'—>E)o (E')dE' 
= —————————— re (6.1) 


SU (E) + SE) | 1 + __£(E) 
: - 1 = £¢E) 


where S(E) is given in its unnormalized form by Eq. (5.27) and, 


=(E'—E) = No(E'—E) , (6.2) 


N being the number density of the moderator. Solution of 

Eq. (6.1) for the spatially averaged, energy dependent thermal 
flux in the moderator is accomplished with Program COUTH 
(Cornell University Thermalization), a Fortran 63 code written 
for use with the Control Data Corporation 1604 digital com- 
puter. A Fortran statement of COUTH appears in Appendix E, 
while the numerical methods employed are discussed in 

Appendix C. This chapter describes the lattice variables and 
integral parameters calculated in COUTH. The results of a 

set of actual lattice calculations are described in Chap. /. 


Panis 
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Determination of the spatially averaged flux in the 
moderator through the solution of Eq. (6.1) leads to the 
calculation of several other useful quantities. The 
spatially averaged thermal flux in the fuel, oO (E), and in 
the cladding, 9° (E), may be determined from the moderator 
spectrum by equating absorption rates in the various 
regions. Utilizing the relationships in Eqs. (4.53) and 
(4.54), 





£(E) | ceoreaen =) 
ae a "| ea Sn a ’ 


v(E) = GH | 
V = 1-£(E) P(E) (1-P, (E)) 


(653)) 


f are the volumes of the moderator and fuel 


where V" and V 
respectively, 51 (E) is the fuel absorption cross section, 
Bee «/ (1-f) is the absorption rate in the fuel element, 
and eGo bee) / CL=P .)P is the ratio of absorptions in the 


fuel meat to absorptions in the fuel element. Similarly 


for the cladding: 


m. a 
Serny - ¥™ Za(®) | £(E) |-[2 I (E) (1 | =n) 





1-f£(E) eee Gy (Clio s) (UB) 


(6.4) 
where V° is the cladding volume and = (E) is the cladding 
absorption cross section. 

Average neutron velocities in each region of the cell 


may now be defined utilizing the relationship between 
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neutron density and flux, 
N(E) = 9(E)/v(E) , (6.5) 


where v(E) is the neutron velocity. The average velocity in 


the moderator, Va? is then expressed as 


Eo Eo 
if v(E)N"(E)dE if © (E)dE 
0 ee. =p Oe 
E 


vv = O 
m E 


Cc Cc 
if N (E)dE if o (E)/v(E)dE 
O O 


Similar expressions define the average velocities in the fuel, 


(6 .6) 


V,, and in the cladding, v,, 
Definitions of disadvantage factors follow directly by 
expressing the total spatially averaged neutron density, N, 


as shown below for the moderator: 


E 
Cc Cc 
N’ -f N'(E)dE -/ O (Ey Cb aka (6.7) 
O 


O 


Using analogous expressions for the total densities in the 
fuel, Ni, and the cladding, hee one can calculate neutron 


density disadvantage factors for the moderator, Bu 


n? 
om Wyn (6.8) 
and for the cladding ans 
so = NC/N (6.9) 
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The disadvantage factors are related to the familiar flux 


disadvantage factors through the relationships: 


Op TOOT = VL On/Ve (6.10) 
and 
oS = o/o* = V0°/¥, , (6.11) 


where m, c, and f again refer to moderator, cladding, and 


fuel, and 
E 
c 
-| o (E)dE . (6.12) 


oO 


“E 


In terms of quantities already defined the thermal 
utilization, f, and the more basic quantity, nf, may be 


determined, where 


thermal absorption rate in the fuel rod 
LO en oor (6.13) 
thermal absorption rate in the cell 


and, 


thermal fission neutron production rate in the cell 
Te 
thermal absorption rate in the cell 


(6.14) 


These definitions lead to the following expressions: 
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/ 52 (E) o: (E)dE 
f= ——9 eee 6.15) 
J 


om (E) E 4 Efi O(E)dE 


and, 
E 
Cc 
Vv vf SF (E) 01 (E)/[1+a(E) ]dE 
nf = “E P (6.16) 
vf 5 (E) E + Hs ) © (E)dE 


O 


where v is the number of fission neutrons produced per 
thermal fission, and a(E) is the energy dependent capture 
to fission ratio for the fuel. 

A thermal neutron lifetime, Ls may be defined as 
the ratio of the total number of thermal neutrons in the 


cell to the thermal absorption rate in the cell, 


Le See ’ (6.17) 


Cc 
al 5™(E) jie mee | OU(E)dE 


O 


where N is defined by Eq. (6.7) for each region. 
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Since one of the prime concerns in reactor calculations 
is the accurate determination of various reaction rates, it 
is desirable to define average thermal cross sections which 
are suitable for use in calculations utilizing the thermal 
flux or the 2200 meter/second flux. In terms of the 
spatially averaged flux or neutron density in a given region, 
a reaction rate, RR, per unit volume of the region is 


expressed by: 


Bo 


E 
@ 
RR | =(E)o(E)dE | =(E)N(E)v(E)dE . (6.18) 
O 


O 


As before, superscripts on the variables will denote the 
region considered (m = moderator, c = cladding, f = fuel), 
while subscripts will denote the interaction considered 

(a = absorption, s = scattering, t = total, tr = transport, 
f = fission). Thus, the thermal absorption rate in the 


moderator will be given by 


E Eo 


Cc 
RR, .| =, (E)O. (E)dE ~ if =, (E)N' (E)v(E)dE (6.19) 
O O 


L/ 


In the notation of Wescott,—’ reaction rates may now be 


formulated in terms of an effective cross section, =, and 


61/ C. H. Wescott, "The Specification of Neutron Flux and 


Effective Cross Sections,'' AECL 352 (1956). 





a 
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a mean cross section, ©, for use in a calculation with the 


spatially averaged 2200 meter per second flux, 9., and the 


0’ 
spatially averaged thermal flux, Ops respectively. A 


given reaction rate is now defined by 


KS = a : (6.20) 
where, 
Ee a 
oe Nv, ee if N(E)dE = v0 | O(E)/V(E)dE , 
O O 
(ons 
and, 
E. E. 
Or = Nv / N(E)v(E)dE = | O(E)dE . (6.427)) 
O O 


The velocity v, denotes 2200 m/sec, while v is the mean 
velocity defined by Eq. (6.6). This leads to the definition 


of the effective cross section as: 


E 


Cc 
f =(E)9(E)dE 
O 


oS (6.23) 


Cc 
vol N(E)dE 
O 


while the mean cross section is defined as 
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E 


Cc 
| =(E)o0(E)dE 


= 40 
S= 2 séF (6.24) 


Cc 
[ o(E)dE 
O 


Consideration of Eqs. (6.23) and (6.24) shows that By, = Sv. 
Applying this nomenclature to absorption rate per unit 


volume of the moderator as an example yields: 


RRO = 25 95 = 24 OF (6.25) 
where 
Eo 
[ 3, (E)9" (E)dE 
= EC (6.26) 
LG 
vo} 9 (E) /v(E) dE 
O 
E 
Cc 
| 3, (E)9 (E)dE 
ee = —____—— ; (6.27) 
C 
| 9 (E)dE 
O 
E 
Cc 
% = Ve if oO (E)/v(E)dE , (6.28) 


O 
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E 
Cc 
Or = if 9 (E)dE , (6.29) 
O 
and, 
am _am- 
24, * ov, (6.30) 


This approach is extended to define effective and mean 


H and sh sub - 


cross sections for the homogenized cell, > 
scripts denoting the interaction considered. In this 


manner 
Eo 
| ae | v™="(E)9"(E) + V°sC(E)SC(E) 4+ visk cEyo* cE) | 
ah O 
5 


E 3 
Cc 


Vv SEY | v™e™(E) + V°O°(E) + ViO*(E) | 


(6.31) 
and, 
Eo 
‘f dE ymsm (Eee) + V°S°(E)o°(E) + vist Erste | 
_ PF 


> = E 


Cc 
| aE] V"G(E) + VOO°(E) + vie CE) | 
O 


(6.32) 





- 
=> 


‘-7 
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The mean velocity of neutrons in the homogenized cell is de- 





fined as 
Bo 
| aE | v"a(E) + V°O°(E) + via ED 
— O 
V = ' 
h BE. 
vol 5 | V"o™(E) + V°OC(E) + viet | 
O 
(6.333) 
so that, 
sh sh— 


To complete the calculation of the thermal parameters, 
one would like to define a thermal diffusion area, L2, E Ot 
the homogenized cell in a meaningful way. In the one 
velocity, one region problem, iL may be defined unambiguously 


as 


Z 


L~ = D/=.,, ; (6.35) 


where D = Ly Sa. However, in the polyenergetic, homogenized 
cell problem, one must decide whether L2 in Eq. (6.35) should 
be calculated utilizing mean or effective cross sections, and 
whether D should be calculated by averaging the transport 


cross section or the transport mean free path. 
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Meeeieres” tas indicated the desirability of utilizing 
the mean homogenized cross section, in the calculation 
of L?. He also notes that several quchore ae argue in 
favor of averaging the transport mean free path, while 
S$}ustrand22/ has shown a physical basis for averaging the 
transport cross section in the determination of D. The paradox 
between these two approaches has recently been clarified by 
Berean ll! and on the basis of that discussion this work 
averages the reciprocal of a homogenized transport cross 
section, sf (E), over a homogenized spectrum in the deter- 
mination of the diffusion coefficient, D. The homogenized 


transport cross section is defined as 


aM (E)O™(E)V™ + 5° (E)O°(E)V® + 5f.(E)9"(E)V" 


h 
2 aie eae ET eT FF ’ 
o™(E)V™ + G(E)VS + O° (E)V 


(6 .36) 
while the homogenized flux, o'(E), is defined by 
yn. | OCEDV™ + BCE)VS + OF (EDV 
SE) = —~_____,——_.. (6.37) 
v@ 4 yo 4 y 
62/ 


—' J. A. Larrimore, "Temperature Coefficients of Reactivity 
in Homogenized Thermal Nuclear Reactors," Massachusetts 
Institute of Technology Thesis (1962). 


A. M. Weinberg and E. P. Wigner, The Physical Theory of 
Neutron Chain Reactors, p. 514, The University of 
Chicago Press, Chicago, 1958. 


64 / C. D. Petrie, M. L. Storm, and P. F. Zweifel, "Calculation 
of Thermal Group Constants for Mixtures Containing 
Hydrogen," Nucl. Sci. Eng. 2, 728-744 (1957). 
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The diffusion coefficient is then expressed as 


Eo 
l a 
3 O° (E) /52_ (E) dE 
. (6.38) 
c 
if O° (E)dE 
O 
The thermal diffusion area is calculated from 
L? = D/=; (6.39) 


where 5? is defined by Eq. (6.32). 

The calculation of all of the parameters and variables 
discussed in this chapter is carried out by program COUTH. 
The results derived from these calculations for an actual 


lattice system are discussed in Chap. 7. 


65/ K. S. Singwi and L. S. Kothari, ''Transport Cross Section 


of Thermal Neutrons in Solid Moderators,'' Second 
Geneva Conference Proceedings, 16, 325 (1958). 
66/ N. G. Sjtstrand, "Definition of the Diffusion Constant 
in One-Group Theory,'' J. Nucl. Energy, Part A, 
Reactor Science, 151 (1960). 
oily G. C. Pomraning, "On the Energy Averaging of the Diffusion 
Coefficient," Nucl. Sci. Eng. 19, 250 (1964). 








CHAPTER 7 
APPLICATION OF THE METHOD TO ACTUAL LATTICES 


7.1 Room Temperature Calculations 

In this chapter we discuss the numerical results obtained 
with program COUTH when the cell theory formulated in Secs. 
2.1 and 2.3 is applied to an actual lattice system. The aim 
of this work is to calculate at various temperatures, the 
Spatially averaged thermal spectra, and the integral parameters 
defined in Chap. 6 for all of the cores incorporated in the 
Cornell University Zero Power Reactor eves ee The 
Cornell ZPR cores are light water moderated and utilize high 
density, low enrichment uranium dioxide fuel, clad in aluminum. 
Several water to uranium oxide volume ratios (nominally 1:1, 
1.5:1, 2:1, 3:1 and 4:1) are provided with sets of grid plates 
which vary the spacing (lattice pitch) between the fuel 
elements arranged in a hexagonal array. More specifically, 
the fuel rods consist of 0.60 inch diameter uranium dioxide 
fuel of density 10.38 gem/cc, enriched to 2.096 atom per cent 


of uz? 


The fuel rods, formed by 48 one inch long fuel 
pellets, are clad in 0.028 inch thick aluminum to form a 
fuel element with an outside diameter of 0.666 inches. A 
list of the parameters of the various lattices is given in 


Appendix A. 


68/ D. D. Clark, ''The Cornell University Nuclear Reactor 


Laboratory," CURL-1, Cornell University (1961). 
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Initial calculations have been carried out utilizing 
as a scattering kernel both the Nelkin water kernel, 
Eq. (5.15), and the monatomic gas model, Eq. (5.11), for 
a gas of mass one. The two physical temperatures of 20°C 
and 40°C have been considered. The results of an elementary 
study utilizing the monatomic gas model have previously been 
peeorred 0! 

In the solution of the neutron balance equation, 
Eq. (6.1), to determine the spatially averaged moderator 
flux, » (E), the heterogeneous effects of the lattice are 
accounted for by the relative fuel absorption, f(E)/[1-£(E)], 
as given in terms of escape probabilities by Eq. (2.29). 
The moderator escape probabilities, P(E) and P(E); are 
given by Eqs. (3.4) and (3.6), respectively. This yields 
Eq. (3.8) for the relative fuel absorption. The fuel element 
escape probability to be used in Eq. (3.8) is calculated by 
means of Eq. (4.21), which in turn requires the calculations 
indicated by Eqs. (4.22), (4.19), and (4.8). Equation (4.8) 
calculates the fuel rod escape probability utilizing the 


parabolic flux approximation. 


69/ S. S. Berg, "Initial Experiments on the Cornell University 
Zero Power Reactor Cores,'' Cornell University Thesis 
(1964). 

/0/ K. B. Cady and C. R. Mac Vean, "A Simplified Formulation 

of Space-Energy Cell Theory," Trans. Am. Nucl. Soc. 6, 
234-235 (1963). 
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The source of neutrons due to down scattering from 
energies above ae that given by Eq. (5.27), normalized 
to one neutron per second per unit volume of the moderator. 


That is, 


ERF E/T, 


S(E) = ==. Gly 


¢ 
if ERF E/T, 
O 


The calculations were performed utilizing an energy 
mesh containing forty-five energy points ranging from 
0.001 eV to a thermal cutoff energy, Ea? of 1.609 eV. In 
the stipulation of a useful energy mesh a compromise must 
be reached between several conflicting requirements. In 
the first place, a reasonable range of energies must be 
covered. Secondly, the total number of energy mesh points, 
n, may not be too large, since the computer time necessary 
for the calculation of the scattering kernel increases as 
n squared. Furthermore, it is desirable to have a finely 
spaced energy mesh in the region above the thermal energy, 
iee= kIT, where the 1/E tail joins that part of the spectrum 
resembling a Maxwellian. Finally, the use of the Nelkin 
water kernel places certain restrictions on the size of 
the energy spacing, AE, between the energy points in the 
mesh. With the use of the approximation allowing for all 
transitions of the hindered rotator, Eqs. (5.17) (utilized 


for incident neutron energies less than 0.33 eV) it is 
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necessary to use a AE which is less than the frequency 

(®, = 0.06 eV) of the hindered rotator. This requirement 
insures that the fine structure of the water kernel is 
preserved in the calculations and that energy integration 
of the scattering kernel will yield a fairly accurate 
scattering cross section. In order to satisfy this require- 
ment, the energy spacings utilized in the calculations are 
less than one quarter of w, Or .O15 eV for all energies 
below 0.33 eV. In light of these restrictions, the energy 
mesh used in the calculations consists of 35 energy points 
below E = .33 eV with 34 energy intervals corresponding 

to equal intervals in velocity. Each of these low energy 
velocity intervals corresponds to 0.1 Vy or 220 m/sec. 

For energies above E = 0.33 eV the ten energy intervals 
correspond to velocity intervals which exponentially 
increase from 220 m/sec to 2200 m/sec as energy increases. 
A complete listing of the energy mesh utilized is given 

in Appendix B. 


235 


The U absorption and fission cross sections used 


iy 


are those compiled by Suich.—' The absorption cross 


238 


sections for U , the moderator and the cladding are all 


assumed to be 1/v and the 2200 m/s values of these cross 
sections are given in Appendix A. The scattering cross 


it / 


—' J. Suich, ''Temperature Coefficients in Heterogeneous 
Reactor Lattices,'' Massachusetts Institute of 
Technology Thesis (1963). 
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sections for all of the materials other than the moderator 
are calculated utilizing the analytic expression, Eq. (5.13), 
derived from the monatomic gas model. Scattering cross 
sections for the moderator are computed by integration of 
the scattering kernel in accordance with Eq. (5.12). The 
use of the monatomic gas model yields a water scattering 
cross section given by Eq. (5.13), while the use of the 
Nelkin water kernel results in the scattering cross section 


given by BNL 325 .22/ 


The free atom scattering cross sections 
used are given in Appendix A. 


The transport cross sections are calculated using 


Sey(E) = 5,(E) + 2,(E) (1 - HEY) . (7.2) 


The average cosine of the scattering angle, u(E), is cal- 
culated from the Nelkin water kernel, Eq. (5.24) and Fig. 
(5.2), for the water moderator. The water transport cross 
section derived from the Nelkin model is also used in the 
calculation of the moderator escape probabilities with use 
of the monatomic gas kernel. For all other materials u(E) 
is given by 

u(E) = pw = 2/3A , (7.3) 
A being the mass number of the scattering material. 


Calculation of the relative fuel absorption f(E)/ 


earth), Eq. (2-29), using the Pep, shown in Fig. (3.1), 


12! y Hughes and B. Schwartz, ''Neutron Cross Sections," 
BNL 325, 2nd Ed., Brookhaven National Laboratory (1958). 
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the Pacey) shown in Fig. (3.2), and the EAD shown in Fig. 
(4.2), yields the apparent lattice absorption cross section, 
5, (E) 1 + £(E)/[1-f£(E)]]|. This quantity is shown in 
Fig. 7.1 for the five Cornell ZPR cores. Comparison with 
the absorption cross section of pure water, =. (E); shows 
the non 1/v behavior of this apparent cross section as well 
as the effect of the 0.273 eV resonances in us? It is 
this non 1/v behavior which yields a temperature coefficient 
of the thermal utilization, f. 

The apparent lattice absorption cross section shown 
in Fig. (7.1) is substituted into the neutron balance 
equation, Eq. (6.1), along with either the monatomic gas 
or Nelkin water scattering cross section. The computer 
solution of Eq. (6.1) is discussed in Appendix C, the result 
of the solution being the spatially averaged moderator 
spectrum, © (E), evaluated at the discrete energy points 
in the energy mesh used. Figure 7.2 shows the results of 
this calculation utilizing the Nelkin water kernel for the 
five water to uranium oxide volume ratios of the Cornell 
ZPR. Since all of the results are normalized to one neutron 
per second slowing down per unit volume of the moderator, 
Fig. 7.2 shows the flux suppression as well as the spectral 
hardening experienced in the tighter (lower water to oxide 
volume ratios) cores. The difference between moderator 


Spectra predicted with the monatomic gas model and the 


Nelkin water kernel is shown in Fig. 7.3 where the 9 (E) 
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Sec. 7.14 
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Figure 7.2 Spatially Averaged Moderator Thermal Spectra 
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Figure 7.3 Comparison of the Moderator Spectra Predicted with the 
Nelkin Water and Monatomic Gas Models 
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predicted by these two models are plotted for the extreme 
lattices of the ZPR; that is the 1:1 and 4:1 water to oxide 
volume ratio cores. Figure 7.3 indicates that harder 
spectra are predicted with the use of the Nelkin water 
kernel than with the use of the monatomic gas kernel. 

Using the spatially averaged moderator flux, the 
spatially averaged fluxes in the fuel, o(E), and the 
cladding, 9°(E), may be calculated using Eqs. (6.3) and 
(6.4) respectively. Results typical of this calculation 
using the Nelkin water kernel are shown in Fig. 7.4 where 
O'(E), O° (E), and 9. (E) are plotted for the ZPR 3:1 core. 
Spectrum hardening as well as depression, giving rise to 
disadvantage factors, is experienced as one goes from the 
moderator to the cladding to the fuel. 

The average neutron velocities in the moderator, 
cladding, and fuel as defined in Eq. (6.6) are shown in 
Fig. 7.5 for all of the ZPR cores, these calculations 
being made with the Nelkin water kernel. Again the 
spectral hardening produced in the tighter cores is 
evident. Figure 7.5 also shows that for the tighter 
cores Vo is well represented by the relationship which 


3 


Foneck— has verified for the Brookhaven uranium metal 


cores, 


73/5, C. Honeck, "The Calculation of the Thermal Utiliza- 


tion and Disadvantage Factor in Uranium/Water 
Lattices,'' Nucl. Sci. Eng. 18, 49-68 (1964). 
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Figure 7.4 Spatially Averaged Cell Spectra 
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V, = 5(v, + ¥,) . (7.4) 
The error in the result predicted by Eq. (7.4) increases 
with the water to oxide volume ratio to a maximum of 1.8% 
for the 4:1 core. 

The average neutron velocities, Vp> for the homogen- 
ized cell, Eq. (6.33), are shown for calculations with 
both the monatomic gas model and the Nelkin water kernel 
in Fig. 7.6. Again the harger spectra predicted by the 
Nelkin water kernel is evident as the average neutron 
velocities for the homogenized cell calculated with the 
monatomic gas model are lower than those given by the 
Nelkin model by 5.4% for the 1:1 core and 2.4% for the 
4:1 core. 

The neutron density disadvantage factors for the 
moderator, hig and for the cladding, ahs as defined in 
Eqs. (6.8) and (6.9) respectively, are shown in Fig. 7.7 
as calculated with both the monatomic gas kernel and the 
Nelkin water model. In this case, the discrepancies 
between the results computed with the two models is 
entirely attributed to the softer spectra produced with 
the monatomic gas model since the transport cross section 
for water used in both calculations was that calculated 
with the Nelkin water kernel. The on predicted with the 


monatomic gas kernel are 1.4% higher for the 1:1 core and 
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0.8% higher for the 4:1 core than those calculated with the 
Nelkin model. Similarly, the differences in the 5 cal- 
culated by the two models is 0.8% for the 1:1 core and 

0.3% for the 4:1 core. The fact that the disadvantage 
factors monotonically increase as a function of water to 
oxide volume ratio is expected! +! Since a non-reflecting 
cylindrical cell boundary condition is incorporated in this 
formulation. 

The results of the calculation of the thermal utili- 
zation, f£, and nf defined by Eqs. (6.15) and (6.16) 
respectively for the Nelkin water model are shown in 
Figs. 7.8 and 7.9. The results predicted with the mon- 
atomic gas model closely approximate those computed with 
the Nelkin water kernel, the differences being less than 
0.2% for all cases of f, 7, and nf. This indicates that 
these integral properties are essentially insensitive to 
the scattering model used in the averaging process. 

Figure 7.10 shows the variation of thermal lifetimes, 
Eq. (6.17), as a function of water to oxide volume ratio 
calculated with the Nelkin water kernel and the monatomic 
gas model. Again relatively small variations, always less 
than 0.4%, between the Nelkin water and monatomic gas kernel 
calculations are present. ‘The lifetimes are seen to decrease 
in the tighter cores because of the increased apparent 


absorption, Fig. 7.1, present in these cores. 


74) H. C. Honeck, p. 56, see Footnote 73. 
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Figure 7.8 Thermal Utilization 
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Figure /.11 shows the effective absorption cross 
section, Sa Eq. (6.31), and the mean absorption cross 
section, : Eq. (6.32), for the homogenized cell as 
calculated using the Nelkin water and monatomic gas 
models. Variations encountered in the calculation of 
these results when the monatomic gas model was utilized 
were less than 0.4% for the effective absorption cross 
sections and 2.0% for the mean absorption cross sections. 

The thermal diffusion coefficient, D, and the 
thermal diffusion length, L, as calculated by averaging 
a homogenized mean free path in accordance with Eq. (6.38) 
are shown in Fig. 7.12 using both the Nelkin water and 
monatomic gas kernels. From Fig. 7.12 it is seen that 
the monatomic gas model predicts a D which is about 3.0% 
less, and an L which is about 4.0% less than those cal- 
culated with the Nelkin water kernel. The differences 
between these results are strictly due to the relative 
hardness of the spectra predicted by the two models, since 
the transport cross sections utilized were the same in 
both cases. It is still not clear how these calculated 
diffusion lengths would compare with those measured in a 
buckling experiment; however, it is noteworthy that the 
minimum in L for the 2:1 core is the same phenomenon 


fou, 


measured for the Brookhaven uranium metal lattices .— 


5/ R. L. Hellens and H. C. Honeck, ''Summary and Preliminary 
Analysis of the BNL Slightly Enriched Uranium, Water 
Moderated Lattice Measurements,'’ Light Water Lattices, 
p. 54, International Atomic Energy Agency, Vienna (1962). 
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This also indicates that the customary receipe for cal- 
culating a homogenized diffusion area, 12 = (1-f) es 16/ 
lacks validity in that it predicts a monotonically 


2? ° ° ° ° 
with increasing water to oxide volume 


increasing L 
ratios. 

A tabular listing of the parameters presented in this 
section appear in Appendix B as calculated with the Nelkin 


water kernel at 20°C. A discussion of the accuracy of 


these results is given in Chap. 8. 


7.2 Temperature Coefficients 


In order to estimate the effect of temperature on the 
various integral parameters defined in Chap. 6 and pre- 
sented for 20°C in Sec. 7.1, all of the calculations were 
remade at a physical temperature of 40°C. In these cal- 
culations any changes in the fuel and cladding densities 
and dimensions have been neglected. In addition, the 
expansion of grid plate supporting the fuel elements leading 
to a change in the equivalent cell radius has also been 
neglected. All of the results are expressed in terms of 
temperature coefficients defined for a given variable, x, 


as 

ie 1 x40) eo 

x AT x (20°) y) G ° (7.5) 
76/ 


—' A. M. Weinberg and E. P. Wigner, The Physical Theory of 
Neutron Chain Reactors, p. 724, The University of 


Chicago Press, Chicago, 1958. 
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The temperature coefficient of the average neutron 
velocity in the homogenized cell, Ve Eq. (6.33), is shown 
in Fig. 7.13 for both models of the scattering kermel. 

This indicates that although the Nelkin model predicts a 
harder homogenized spectrum at room temperature, the 
monatomic gas spectra hardens at a faster rate with 
increasing temperature. It is therefore possible to 
postulate that at some elevated temperature the integral 
parameters dependent primarily on the degree of spectral 
hardening will be the same for both models of the scattering 
kernel. The temperature coefficient of Vp predicted with 
the monatomic gas model is 23.8% higher for the 1:1 core 
and 7.8% higher for the 4:1 core than those predicted by 
the Nelkin water kernel. 

Figure 7.14 shows the temperature coefficients of 
the moderator and cladding disadvantage factors, oF and 
m 
(6.9) respectively. As expected, they are negative; the 


these parameters having been defined in Eqs. (6.8) and 


effects of harder spectra being seen in Fig. 7.7. The 
monatomic gas model predicts larger negative temperature 
coefficients because of the greater relative spectral 
hardening encountered with that model. The maximum differ- 
ence between the scattering models occurs in the 1:1 core 
where the temperature coefficient of oF predicted by the 
monatomic gas model is 30% higher than the prediction by 
the Nelkin water kernel. In the case of a the results 


differ by 20%. 
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Figure 7.13 Temperature Coefficient of the Average Neutron 
Velocity for the Homogenized Cell 
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The variation of the temperature coefficients of the 
enermal utilization, f, Eq. (6.15), and HE; Eq as(o.06) mas 
shown in Fig. 7.15, the difference between the two curves 
for each scattering model being equal to the temperature 
coefficient of 7». Again the results predicted are more 
positive with the use of the monatomic gas model because of 
the greater relative spectral hardening given by that model. 
It is to be noted that the temperature coefficient of nf 
is slightly negative for the 1:1 core. 

The final temperature coefficients worthy of note are 
those of the diffusion coefficient, D, Eq. (6.38), and the 
thermal diffusion length, L, Eq. (6.39). These are shown 
in Fig. 7.16; the difference between the monatomic gas and 
the Nelkin water models being quite evident. Again, the 
temperature coefficients predicted by the monatomic gas 
model are higher than those produced by the Nelkin water 
kernel. The deviation is the most pronounced for the 1:l 
core where the monatomic gas results are higher than the 
Nelkin water kernal results by the following percentages: 
temperature coefficient of D, 68%; and temperature co- 
Peticient of L, 34%. 

The differences between the monatomic gas and the 
Nelkin water kernel results are more pronounced for the 
D and L temperature coefficient calculations than for any 
of the other parameters previously discussed. This may 


be attributed to the fact that these calculations are 
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Figure 7.15 Temperature Coefficients of the Integral Prop- 
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Figure 7.16 Temperature Coefficients of the Thermal Diffusion 
Coefficient and Length for the Homogenized Cell 
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doubly sensitive to the degree of spectral hardening pre- 
dicted by the two models, since the homogenized cell trans- 
port cross section, se (E); Eq. (6.36), used in these cal- 
culations is derived by flux and volume weighting the 
transport cross sections of the various regions in the 
cell. The reciprocal of this homogenized transport cross 
section is then averaged over a homogenized spectrum, 
Eq. (6.37), making the resultant D highly sensitive to 
the degree of hardness in the spectra utilized in these 
averaging processes. 

An estimate of the accuracy of the temperature 
coefficient calculations is discussed in Chap. 8, while 
a tabular listing of all of the temperature coefficients 
calculated with the Nelkin water kernel is given in 


Appendix B. 
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CHAPTER 8 
ACCURACY OF THE RESULTS 


In an attempt to estimate the accuracy of the results 
predicted with program COUTH, as presented in Chap. 7 and 
listed in Appendix B, two lines of approach have been 
utilized. The first of these involved the calculation by 
COUTH of several integral properties of the Brookhaven 
uranium dioxide lattices, a comparison then being made between 
these results and those calculated by fenee “en the 
THERMOS transport theory code. The second of these was an 
investigation of the sensitivity of the Cornell ZPR results 
to changes in the inputs to the COUTH code such as the den- 
sities, cross sections, and dimensions of various regions 
within the lattice cell. 


78/ 


The Brookhaven cores——' studied are generically related 
to the Cornell ZPR cores in that they are light water mod- 
erated and utilize low enrichment uranium dioxide fuel. 

Five water to uranium oxide volume ratios are provided, 
these being 1.291:1, 1.600:1, 2.050:1, 2.803:1, and 4.000:1. 
The fuel rods are 0.444 inches in diameter and clad with 
0.028 inch thick stainless steel giving the fuel element 

an outside diameter of 0.500 inches. The uranium dioxide 


fuel is enriched to 3.0 atom per cent, its density being 


9.3 gm/ec. 


Lith) H. C. Honeck, "The Calculation of the Thermal Utilization 


and Disadvantage Factor in Uranium/Water Lattices," 
Nucl. Sci. Eng. 18, 49-68 (1964). 


-102- 


Ghap. 6 103 


The integral thermal parameters compared are the thermal 
utilization, f, the ratio of the number of fission neutrons 
produced to the number of thermal neutrons absorbed in the 
fuel, », and the moderator neutron density disadvantage 
factor, one The THERMOS calculation of these parameters 
utilizes a thirty point energy mesh ranging from 0.00025 to 
0.632 eV, while COUTH uses a 41 point energy mesh ranging 
from 0.001 to 0.645 eV, as discussed in Sec. 7.1. Both 
programs utilize the same cross section input data and both 
perform their calculations utilizing a Nelkin water model 
for the scattering kernel at 20°C. 

A comparison between the results predicted by THERMOS, 
using integral transport theory, and COUTH, using the 
simplified cell theory, are shown in Fig. 8.1 for the thermal 
utilization, f. In Fig. 8.1 it is seen that the f's cal- 
culated with the two codes compare very well, the maximum 
deviation in the results occurring in the 4.0:1 core where 
the THERMOS f£ is 0.39% higher than the £ calculated with 
COUTH. 

The results of the calculation of 7 by COUTH and 
THERMOS are seen in Fig. 8.2 to be even closer than those 
for f. In this case the maximum difference between the 
results predicted by the two techniques is found in the 
1:1 core, the THERMOS result here being 0.19% higher than 


the COUTH estimate. 


18/ H. C. Honeck, see Footnote 7/7. 
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Figure 8.1 Thermal Utilization 
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Figure 6.2 Integral Property, 7 
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In the case of the moderator disadvantage factor, anes 


the agreement between the two codes is not quite as close 
as in the calculations of f and 7. Figure 8.3 shows the 
results derived, indicating that THERMOS predicts higher 
disadvantage factors for the tighter cores than COUTH, while 
this role is reversed with the larger water to oxide volume 
ratios. The THERMOS oF for the 1.29:1 core is 02627 inalghes: 
than the COUTH number, while for the 4.0:1 core the COUTH 
oF is 0.97% higher than that given by THERMOS. oneet 
assigns a combined uncertainty to the THERMOS calculation 
of the disadvantage factor of 1.5%, meaning the COUTH cal- 
culation of this parameter has an uncertainty less than 
2.5%. The deviations between COUTH and THERMOS in this 
calculation are believed to be attributable to COUTH's use 
of simple diffusion theory in the moderator region, com- 
bined with an overestimation of the fuel element escape 
probability. Resolution of these differences requires 
further study and is discussed in Chap. 9. 

The parametric study made on the Cornell ZPR for 20°C 
using the Nelkin water kernel investigated the sensitivity 
of the calculated parameters to changes in the input data. 
Changes where made in the fuel density, the fuel enrichment, 
the fuel rod radius, the cladding radius, the cladding 
Scattering and absorption cross sections, the moderator 
density, and lattice spacing or pitch. In an attempt to 


dle) H. C. Honeck, p. 60, see Footnote 77. 
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determine the effect of using various numerical integration 
schemes, the results utilizing an integration method which 
fit exponentials, i.e., straight lines on log paper, were 
also studied. The effect of these changes were noted for 
several parameters, these being: The thermal utilization, 
f; the ratio of the number of fission neutrons produced to 
the number of neutrons absorbed thermally in the cell, nf; 
the moderator neutron density disadvantage factor, 53 the 
thermal diffusion length, L; the thermal lifetime Les the 
mean neutron velocity for the homogenized cell, Vp3 and the 


mean total cross section for the homogenized cell, xh The 


= 
results of the parameter variation study are shown in 
Table 8.1 where the percentage changes in the input para- 
meters required to produce a positive 0.1% change in the 
resultant integral properties are listed. The results 
presented in Table 8.1 represent the most sensitive cases, 
the error analysis having been carried out on the 1:1 and 
4:1 ZPR cores. Thus the percentage changes shown which 
involve changes in the fuel element parameters come from 
the 1:1 core analysis, since a relatively large portion of 
the unit cell is composed of fuel for this core. The l1:l 
core is also most sensitive to changes in lattice pitch. 
On the other hand results from the 4:1 core are given for 
the sensitivity of parameters to changes in the water 
density, water playing a relatively inore important role 


in thre core . 
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The use of the numerical integration scheme involving 
the fitting of exponentials rather than straight lines as 
in the trapezoidal integration rule produces the following 
maximum percentage changes in the integral parameters: 


£, 0.004%; nf, .009%; 5", 0.12%; L, 1.1%; &£., .069%; 


n ee 


0.12%; and =a 2.5%. The large percentage change in 


Vy te 
: may attribute to the fact that the numerical integration 
scheme using exponentials predicts scattering cross sections 
for water by integrating the scattering kernel which are as 
much as 3.0% higher than the experimentally measured 
values. 

Table 8.1 shows that the numerical results are 
essentially insensitive to the cladding cross section 
data, while they are most sensitive to the lattice pitch. 
The actual tolerances iver for the Cornell ZPR lattice 
pitches are less than 0.10% so that the output parameters 
will vary no more than about 0.24 because of variations 
in this dimension. The fuel element parameter, density, 
enrichment, fuel rod radius, and fuel element radius are 
known to within 0.5% and reference to Table 8.1 indicates 
that in general variations of the thermal integral prop- 
erties will be less than 0.5%. 

In light of the results of comparing the predictions 
of the THERMOS and COUTH codes and the results of the 


input parametric studies, the following probable uncer- 


tainties may be assigned to the 20°C calculations performed 








Chap. 3s 


Jeet 


by COUTH utilizing the Nelkin water kernel as a scattering 


model. 
ia 
De 


je oemercyeta: Gapes 40 yh. 

Disadvantage factors for the moderator and 
eladding,.2607. 

Thermal lifetimes; 3.0%. 

Mean velocities; 3.0%. 

Averaged cross sections, including the calculation 


of Deando) +s 542 


With regard to the room temperature calculations made 


utilizing the monatomic gas model, the deviations from the 


Nelkin water kernel results were quoted in Chap. 7. Using 


these percentages, one can then assign the following prob- 


able uncertainties to the 20°C monatomic gas calculations. 


die 


2 
3 
4. 
> 


ian manana ORS. 

Moderator and cladding disadvantage factors; 4.0%. 
Thermal lifetimes; 4.0%. 

Mean velocities; 8.0%. 

Averaged cross sections, including the diffusion 


coefficient and length; 10.0%. 


On the basis of these two sets of uncertainties, one 


must reconcile the accuracy obtained with the computer time 


required for a given calculation. A monatomic gas scattering 


kernel, in the form of a 45 by 45 matrix takes only two and 


one half minutes of computer calculation, while the Nelkin 
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water scattering kernel takes on the order of fifteen 
minutes. Thus if one is undertaking a detailed tempera- 
ture study wherein multiple scattering kernels must be 
calculated, use of the monatomic gas model is advocated 
for one primarily interested in f and 7f. However, if 
emphasis is to be placed on the other integral parameters, 
and in particular the averaged cross sections, then the 
more complicated Nelkin model should be used. 

The assignment of uncertainties to the temperature 
coefficients presented in Chap. 7 presents a slightly 
different problem. In this case one must estimate the 
uncertainties in a calculation at a given temperature and 
then estimate the uncertainty in the degree of spectral 
hardening predicted by a certain model. The fact that 
the THERMOS and COUTH results at room temperature are 
in agreement for all values of water to fuel volume ratios 
indicates that COUTH can adequately handle spectral harden- 
ing effects if the scattering kernel model can predict the 
appropriate temperature changes in the energy transfer 
process. We therefore, apriori, assign uncertainties to 
the temperature coefficient calculations as being twice 
those stated above for the room temperature calculations 
with both the Nelkin water model and the monatomic gas 
kernel. 


80/ ¢ S. Berg, "Initial Experiments on the Cornell University 


Zero Power Reactor Cores,'' Cornell University Thee? - 
ESiacn) = 








CHAPTER 9 
SUMMARY AND CONCLUSIONS 


A simplified polyenergetic cell theory has been formu- 
lated for use in the determination of the spatially averaged 
energy dependent thermal fluxes in the various regions of a 
lattice cell. The derived spectra are then utilized in the 
calculations of the thermal integral parameters and average 
cross sections required for reactor computations. 

The cell theory, as formulated, postulates an infinite 
moderator region with the absorption cross section of this 
region appropriately modified to account for the neutron 
leakage into and absorption by the fuel element in the 
actual lattice. The modifications to the moderator absorption 
cross section, have been formulated both in general terms 
an in terms of moderator and fuel element escape probabilities, 
the latter approach offering physical transparency and ease 
ee vcalcullation. 

Analytic expressions for the escape probabilities have 
been presented, integral transport theory having been applied 
in the fuel element region, while diffusion theory was utili- 
zed in the moderator region. Using these analytic expressions 
the theory was applied to actual lattices in the form of the 
Cornell University Zero Power Reactor cores to determine room 


temperature parameters and their temperature coefficients. 


ee 
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In an attempt to estimate the accuracy of these results, 
the method was also applied to the Brookhaven uranium 
dioxide cores and the results then compared with those 
predicted by multigroup transport methods. The remarkably 
good comparison between these two sets of calculations in 
conjunction with a parametric study which yielded a relative 
low sensitivity of the calculated parameters to changes in 
the input data has led to the assignment of rather small 
uncertainties in the results calculated with the simplified 
cell theory. 

Although the estimated uncertainties tend to instill 
confidence in the calculational results, one must insure 
that these events are not just fortuitous and that the 
simplified cell theory is in fact predicting behaviors 
which closely approximate the physical processes taking 
place. 

In this regard, there are three areas in which further 
work is necessary in the simplified cell theory as pre- 
sently formulated. The first of these concerns the validity 
of the assumptions made in the calculations of the modera- 
tor escape probabilities, or alternatively, in the calcu- 
lation of gradient of the moderator flux at the fuel 
element surface. In these calculations it was assumed that 
the shape of the flux in the moderator was given by mono- 
energetic diffusion theory applied at each energy point. 


Although multigroup transport calculations have given some 
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support to this assumption, its accuracy can certainly be 
questioned particularly in the tighter cores where fuel 
element spacing becomes comparable with a mean free path. 
Thus there is a need to accurately determine the flux 
gradient, or more precisely, the net current as a function 
of energy at the fuel element moderator interface. Any 
study such as this should include the effects of a non 
uniform spatial source of slowing down neutrons. The use 
of multigroup transport theory specificially oriented 
toward this problem is suggested for this future study. 
The second area of concern involves the calculation 
of the fuel rod escape probability. In the expression 
utilized in these calculations an isotropic flux was 
assumed incident upon the fuel rod. This in turn leads 
to escape probabilities which are too large; the flux in 
an actual lattice is more peaked toward the rod. The use 
of smaller escape probabilities for the fuel rod would 
lead to larger values of f and 7f, along with larger 
values for the calculated disadvantage factors. Thus 
there is need for at least a parametric study utilizing 
all of the analytic approximations for the fuel rod 
escape probability presented. One could then turn his 
attention to variations caused with changes in angular 
distribution of the neutrons incident on the fuel 


element. 
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The final area where further study is needed concerns 
the calculation of temperature coefficients. The immediate 
need is to compare temperature coefficients calculated with 
integral transport theory with those predicted by the 
simplified cell theory as was done with the room temperature 
case. However, if both calculations utilize the same 
scattering model, it is believed the comparison of the 
temperature coefficient computations will be comparable 
with those for room temperature, since the simplified cell 
theory has demonstrated its ability to handle spectral 
hardening in five different cores at 20°C. The accuracy 
of the temperature coefficient results thus reduces to a 
question of how well the scattering model predicts spectral 
hardening as a function of temperature. This question is 
perhaps best answered by the comparison of pure moderator 
spectra predicted by the model at various temperatures 


with those measured experimentally. 
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APPENDIX A 


LATTICE PARAMETERS FOR THE CORNELL UNIVERSITY ZERO POWER 
REACTOR 


This section details the input parameters which were 
utilized by program COUTH in the calculations described 
in Chap. 7 for the Cornell University Zero Power Reactor 
cores. Appendix D deals with the Fortran FORMATS used in 
COUTH for reading data into the program. 

Table A.1 lists the fuel element parameters for the 
Cornell ZPR while Table A.2 gives the dimensional param- 
eters of the five hexagonal lattices of the ZPR. Material 
concentrations and cross sections are given in Table A.3. 


The absorption cross sections, o listed are evaluated 


a’ 


at E, = .0253 eV while the scattering cross sections given 


are the free atom values, Oe. The number of neutrons 


produced per thermal fission for u?>, v, is 2.43 while 


the value of the y2? fission cross section at Eo is 


81/ 


580.9 barns. Appendix B gives a complete listing as 


a function of energy for the scattering and transport 


cross sections of water, the absorption, fission, and 


Zo5 


scattering cross sections of U , and the scattering 


238 


cross sections of U and oxygen. 


81/ H. C. Honeck, "The Calculation of Thermal Utilization 


and Disadvantage Factor in Uranium/Water Lattices," 
Neel sea. Eng. 18) 66 (1904). 
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TABLE A.1L 
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CORNELL UNIVERSITY ZERO POWER REACTOR FUEL ELEMENT 


PARAMETERS 2/ 
Fuel Material 
Fuel Rod Diameter, in: 
Fuel Density, em/cc: 
Fuel Enrichment, Weight Percent U 
Fuel Enrichment, Atom Percent ue. 
Oxygen to Uranium Number Ratio: 
Fuel Molecular Weight: 
Cladding Material: 
Cladding Density, em/cc: 
Cladding Inner Diameter, in: 


Cladding Thickness, in: 


Fuel Element Diameter, in: 


a/ 


KG 52 


Uranium Dioxide 


0 .600+0.003 
1073S T20e57. 
2207 O20 eee 

2 .096+0 .90 
20a 
270.26(physical) 
6061-76 aluminum 
270 

0.610+0.002 

ORO 7 S=0 008 


0.666+.006 


S. S. Berg, "Initial Experiments on the Cornell University 


Zero Power Reactor Cores,'' Cornell University Thesis 


(1964). 
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TABLE A.3 


CORNELL UNIVERSITY ZERO POWER REACTOR MATERIAL CON- 
CENTRATIONS AND CROSS SECTIONS 


Maeind Concentration 9, (E,)2! o A! 
(atoms /barn-cm) (ea ermtey 
y23> Roca” 679.0 10.0 
y238 2.2654x1072 2.71 8.3 
Oxygen in UO, 4 6510x1074 0.0 3.76 
Pain’ nun 6 .0205x1072 0.230 1.40 
Hydrogen in H50 Oe3525 20.4 
20°C S67 ele 
40°C 6 .6346x1072 
Oxygen in H50 Oo Bee 6) 
20°C 3.3373x107- 
40°C 3.3173x10~2 


a/ H. C. Honeck, "The Calculation of the Thermal Utilization 
and Disadvantage Factor in Uranium/Water Lattices," 
Nuc? Sei. Ene. US. c0, Goa. 
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APPENDIX B 


NUMERICAL RESULTS FOR THE CORNELL UNIVERSITY ZERO 
POWER REACTOR 


Chapter 7 discussed some of the results of the cal- 
culations made on the Cornell University Zero Power Reactor 
using the simplified cell theory as incorporated in COUTH. 
This section presents all of the results in a tabular form 
derived from the calculations utilizing the Nelkin water 
model for the scattering kernel. 

Table B.1 lists the energy, velocity and lethargy 
mesh used in the calculations by COUTH. The derivation 
of this mesh is discussed in Sec. 7.1, and all of the cal- 
culations are carried out using the lethargy mesh. 

Table B.2 gives the scattering kernel in matrix form 
as derived from the Nelkin water model, Eq. (5.15), for 
a physical moderator temperature of 20°C. The scattering 
matrix is expressed in barns/unit lethargy and it is 


related to the scattering matrix in barns/eV by, 


ae (B.1) 


o(u,—= 1; ) = E,O(E7—~ Ey 


The subscripted indices i and j correspond to the lethargy 

and energy points listed in Table B.1. The diagonal terms 

of the scattering matrix listed in Table B.2 are those cal- 
culated with Eq. (5.22) and they are not modified in 


accordance with Eq. (5.26). 
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TABLE B.1 


LZ 


ENERGY, VELOCITY, AND LETHARGY MESH UTILIZED BY COUTH 


_J. ENERGY ELTA EF VELOCITY DELTA V tETHARGY DELTA U 
4 4-699E 000 3.782E-001 17547.0 22n0 oO 0000 2679 
2) EE ce 6445- 5347 747.5 679 2418 
3 9.6A7F<901 1.873F-091 13599.5 esac .5097 ©2153 
= dine 7 
2 6-45nE-n01 9.770E-0n2 11108.7 875.8 9143 ~1642 
e — 9 = 4 
7 4.754E-901 5. 350E-092 9537.2 552.6 1+2194 ©1194 
Bol = = 
9 3.847E-901 3.051F-092 8545.6 348.7 44389 ~0833 
41 3.279E=n01 4.796E-092 7929.0 290.0 1.5940 60563 
12 S-099E=“nol 1.746E=-072 7700-0 220-0 1.6473 0580 
a3 9.904E=-n01 1.695E-092 7480.0 290-0 4.7083 °0597 
15 2.500E-901 1.594F-092 7040.0 2290.0 1.82465 .0635 
4 = ba 
sl v4 2.277E-901 1.493E-002 6600.0 220.0 1.9556 -0678 
e tad 44 
19 4-9R3E-001 1.391E-092 6160.0 290.0 2.0936 s072en 
2 2 -Of 
ragl 1-74 9E=-n01 1.290F=-0902 5720-0 220-0 2.2448 ~0784 





22. _4.594F-n01 4 240E-092 5509.0 22 Qi De, 8 Dik Bie 














23 1°-S57E-n01 1.189E-0902 5280.0 220-0 2.401419 ©0851 
24 4-S%RE=901 = a 
25 4°-204E-n01 1.088E-092 4840.0 220.0 2.5759 20939 
26 41+1446E=-4 7E- 2230.0 2.6690 °#£.0976 
27 1°0127E=-n01 9.866F-093 4400.0 220-0 2.7646 1026 
28 9:129E-n02 09.360F-003 418n.0 220-0 2,86914 ~10814 
09 8.194E=-992 8, 854E-093 3969.0 220.0 2.9773 21148 
30 7-344E-902 8.346F-993 3749.0 220-9 3.0946 1212 
34 6.476F=4)92 7.842E-093 3529.0 220-0 3.21298 01294 
32 5.609E-992 £47.336F-0n3 £#3300.0 ) 449 80 
33 4.95RE-902 6.830E=-0n3 3080.90 220.0 3.47909 14982 
35 3.643E=-992 5. 818E-093 2649.0 220-0 3.7882 -1749 
36 3.9A1E= - 96 
37-2. 5XnE-n02 4.807F-0n3 2200-0 220-0 4.1599 2107 
38 = = 626 5 
oe 1:649F=#902 3.795E-093 Seaee 0 eae 0 4.5994 02674 
40 1-24nE=n02 3,289E-0n3 1549.90 2290-9 4.8642 ~-3083_ 
41 9-107E-903 9.783E-0903 1320.0 270-0 5.1745 ° 3646 
42 6-304F-n03 ?.277E-0n3 1170.0 220-0 5.2001 ©4463 
43 4-N48E-993 1.771E-093 889.0 220-0 5.9854 5754 
44 2.277F=n038 4-265F-093 660.9 2290-0 6.568 -8109 
45 1°-047F-903 ~ONGE Onn 449.0 -0 72.5747 -0000 
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Appendix B 128 


Table B.3 gives the Nelkin water kernel scattering 
matrix for a moderator physical temperature of 40°C. Its 
format is the same as that described for Table B.2. 

Table B.4 lists as a function of energy the scattering 
cross section, Eq. (2.2), and the transport cross section, 
Eq. (5.25), for water. The energy dependent scattering 


cross sections for y23>, y238 


» and oxygen are also listed 
in Table B.4 as calculated using Eq. (5.13) with the param- 
eters listed in Table A.3. 

Table B.5 gives as a function of energy the absorption 
and fission cross SE er ey for ue??, as well as the 
uranium dioxide fuel fission and absorption cross sections. 
One plus the capture to fission ratio, [l+o(E)], for the 
fuel is also listed in Table B.5, along with the normalized 
slowing down source of neutrons, S(E), Eq. (7.1 ), used in 
the calculations. 

Table B.6 through B.10 list the energy dependent 
thermal spectra calculated by COUTH for a moderator 
physical temperature of 20°C. The moderator, Eq. (6.1), 
fuel, Eq. (6.3), and cladding, Eq. (6.4), spectra are 
given for the five Cornell ZPR cores. 


Table B.11 gives the integral properties defined in 


Chap. 6 for the five ZPR cores at a temperature of 20°C, 


82/ J. Suich, "Temperature Coefficients in Heterogeneous 


Reactor Lattices,'' Massachusetts Institute of 
Technology Thesis (1963). 
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Appendix B 136 
TABLE B.6 
SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 
ZPR 1:1 CORE 
I ENERGY MODERATOR CLADDING FUEL. 
SPECTRUM SPECTRUM SPECTRUM 

a 1.-6093E 00 5.4517E-91 5.45499E-01 92 387SE-01 

2 1s2o11e OD 6.158S3E-91 6.1323F-01 6.9196E-01 

3 9.6667E.01 8.0549E-01 7.99149E-O1 7.80649E-01 

4 7.7941E-04 1.005SE Go 9.98499E-01 9.7670E-O01 

5 &.4591E-0)4 1.2364E 90 1.22656 00 1°-1975E O00 

6 5.4731E-04 1.4605E 90 1.449626 00 1-4084E 00 

r| 4.7542E-04 1.6733E 00 1.65266 00 1.-6032E 00 

8 4,29 2E= 01 1.8733E 0 1.8446E 09 1.7B814E 00 
maeCtéi‘“‘C‘ SC ITER 2.0569E 0 2.q188F 00 1-9400E 00 
10 3.9119E-01 2.22751E 00 2.1/7295E 09 2-n78SE 00 
11 3.2786E-04 2.4269E 90 2.3654 00 2.2485E 00 
12 3.0990E-0141 2.2621E 90 2-4900F O00 2e3562E 00 
43 8 9244E-04 = O7HVE 00 = 2 5 98GE OO 2.4525E 00 
14 9,.7549E-01 2.8720E 90 2.78568E 00 2.6278E O00 
15 2.5995E-01 S.1Sloc 0 S.O4U0E 00 | 2.8708E 00 
16 9.49311E-01 3.3868E 90 5.2923E 00 $-1152E 00 
17 2.2/768E-01 3.6805E 00 $.5815-€ 00 3-3940E 00 
18 9.1275E-01 4.0585E 00 $.9501E O00 3-7445E 00 
19 1-983SE-01 4.5281E 90 4.4053F O00 4.1731E 00 
20 4-8442E-01 5.0873E 00 4.9465E— 00 4.6815E 00 
21 1.71n1E~-01 5.7628E 00 5.59826 00 (5.2910E 00 
22 1.5811E-01 6.6162E Q0 6.42U7E 00 6.9586E 00 
23 1.4571E-014 7.7219E G0 7.4837F 00 7.9480E 00 
24 1.5382E-01 9.1329E 90 8.a386E 00 8.305S5E 00 
25 1.2244E-01 1.0860E 91 1.9490EF 01 9.B2em = 00 
26 4-1156E.-01 1.2SS0e Of 1.2422- O1 1.1599E O1 
27 1.0119E-014 1.92S59E O21 1.4646E O01 1.3619E 01 
28 9.13525E-02 1.7990E ol 1.7259EF O1 1.5961E O01 
29 R.1965E-02 2.1205E 91 2.9247€ 01 4.8650E O14 
50 7.S110E=-02 2.4876E gi 2.3666F O01 2.1683E O01 
34 6.4762E-02 2.8784E o1 2.7259 O1 2.4807E 01 
32 5.6990E-02 3.2558E g1 3.96/76E 01 2-7716E O01 
3S 4.9583E-02 |§$.5759E 01  $-3489F 01 S-N001E 01 
34 4,.27/53F-02 $.81S2e O41 $.5455c¢ U1 3.1375E O01 
$5 3.6499E-02 3.9529E oi S.64U5E 01 3.4844E 04 
36 2. 0670E-02 3.9848E o1 $.6256E€ 01 3.1277E O1 
37 ?.929BE-02 o. Sone te. — $.49491E 01 —6-209651E 01 
38 9.0491E-02 3, 6699EF i 3.24456 01 £=2.6991E 01 
39 1.6191F-U2 S.5619E 01  2.8878F 01 2-34356E 01 
40 4.2996F-02 2.9266F€ gt 2.44266 01 1.9203E 01 
44 Cen Ter =U S 2. ee OM 1.94949F O1 1-4659E o1 
42 6.5244E-03 1.8522E€ g1 1.4278€ 01 °° °®&42«©£1-9146E O01 
43 4.0476F-03  i1.2954E o1 9.3837E 00 6.1077E 00 
44 9.27h8F-03 | 7.8840E g0 5.9262F 00 2-9522E 00 
45 1.0119E-03 _3.7721E 00 2.1714F 00 9.4648E-01 
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TABLE B./ 


ys yi 


SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 
ZPR 1.5:1 CORE 















































I ENERGY. MODERATOR 
SFECTRUM 
1 1.60983E000 Saeed 
2 ieee ort te U4) 6.2156E-01 
3 0.6667E-U1  8.1724E-91 
4 7./941E-014 1.-9227E GD 
5 &.4501F-01 1. 260Ce. ae 
6 B.4731E-01 1.494S5E 90 
a _.4.7/5426-01 ie anOne an 
8 4,2192E-U1 1.9S5c2E GD 
9 Be ese 2.1327E 00 
10 3 com ieee 1 2. ST90E Go 
tal 3. 27ROHE RY 2.-5420E Qo 
12 Z.0990E-61 2.6954E og 
13 2.9244E-U1 2. 8287E 00 
14 9.175496-04 S.U463E Ga 
15 . Baoepce<01  S.3362 Gm 
16 Aas te 4 3.6304E go 
5AM pe islet sl 5.97 73e e 
18 9.1275F-U1 4.4366E 90 
19 41 Fee Sei Saea7e iG 
20 4.8442F-04 SIT 752cF Tt 
et DS Ae MESA 6.6668E 90 
22 1.59811E-01 7.8577E go 
23 {eco ee 4 9.435 i 
24 1-0 SRZE HU 1.1493E o1 
25 4.2244F-01 1.4978E gi 
Ze 4;711506-01  ilaveie oe 
27 1.0:119F-01 2.0941E gi 
28 - SSeS =e 2. 2iemee o1 
29 R,1965F-602 5.f615E 1 
30 7 wo ee Ue 3.669254 01 
on 6.4762F-62 4.5973E gt 
$2 5, 6999F-U2 4,9486E Qg1 
33 4,.95R3F 2 Sse 2s el 
34 4.2753+~02 5.9801E 01 
35 SpCeee U2 §=6 -egeee ans 
36 3. U610E-02 6. 4S68E G1 
37 si (sti S29KF-U2 = 6. 3762E 91 
38 Oya tee 2 OG Ee it 
$9 4.6191F-U2_ So JO oe 01 
40 4.2396E-b2 4.9154E gi 
44 C.1072FE-U3 4,.Uu71¢2eF Ql 
42 4.38244E-03  3.1386E G1 
4.0476F-03 Beene Gal 
44 2.27/48F-U3 tf. oss jl 
45 4.0119F-03 — 6.5281E 0 


CLANPDING 
SPECI RUM 


5447/9 24. 
O77 64 5E = 0m) 
O.enGSo3F-Ui 
1.9155F Ud 
1. oaeSse G6 
1 e477 Se OH 
1.-49454F Ug 
1.-FOSBB8F 9 
Oo.nBO9FE Udo 
2.9592F O00 
2.46/3F U0 
2.GU75F Ug 
2aZvoter i100 
2.-9399F 9 
3.92S50F Ud 
J.oLiZe 09 
S.R8511F UO 
4.99576 C4) 
4.,R862NEF UN 
SS eics > ef) 
6.4577" 0 
7.5773F us 
¥.j828F UO 
Teer ail 
1.5491 F 
1-AS44E U1 
LOS o7 Fe 04 
2.49725 vl 
2.88Y5F Ui 
S.4509F G1 
.n2e2F bi 
42o0Ser U4 
Senos » td 
5. 490 Sr Se 
2.6618F Ui 
D9 AGJELE WI 
5.520¢7F ul 
2.14450F Ul 
4.58/46 Ut 
eB Pee = UL 
sonore Uy 
oP 512F U1 
~-4497F C4 
Pa cee ver al Se Vie 
Toe oe jt 


¥WzqQ ~ rN GW &® 


FUEL 
SPECTRUM 


5.-4101E-01 
6.9706E-01 
7.91 08E-O01 
9.9151E-01 
1.919SE 00 
1-43587E 09 
1.643S38E 9 
Le BSS 5 EG 
2.9055e 09 
2.15822 Oi 
2+ 3454E UD 
2.45 702 0 
2-5776E QQ 
2-7726E QQ 
3-74S6E GQ 
5$.-322S5E G0 
S.-44950E 09 
4.n72ece Wi 
4.4057E 1G 
5.29558E 00 
6.08455 UG 
Jo Se OE 
GE.5541E gg 
1.9°%37CE 07 
1-29541F G1 
1-5553E Ot 
1,.85S0E 91 
Zoe fea 
2-AS517E U1 
5.1 50CE sud 
RP teltieons | Cal 
4.4469E G1 
4.5509E U1 
9.87516 U1 
9.9524E G1 
4.910SE O1 
4.46866E Q1 
4. 3E 27 ed 
oS. 7/2 oUe mol 
5. 7499E 01 
os the Comp d 
1.555985 1 
9.3960 QU 
C.436SE Vy 
1.3728E QG 
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TABLE B.8 
SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 
ZPR 2:1 CORE 
i ENERGY MODERATOR ~ CEAnDING. FUEL 
SPECTRUM SPECTRUM SPECTRUM 

il 1.0093—E un  j5.4898E=01 .S.48876=501 # Saeed 

2 4.2511F UO 6.2448E-g1 6.908/7F-01 6.9946E-01 
Se 9.6667F-01. 8.232860).  teereaii 7-9606E-01 

4 72/7 9o 1. S0ase Git 1.n209F O00 9.9865E-01 

5 6.45n1F =14 1.2750 GD _ duped SiGEs OD 1-2298E 00 

6 5.4751 EWU1 1.97 Poe Ho  -2.4924E 00 1.-4534E 00 
as ate! SAE a lage mean 1.7148F 00 «26 663519 O10 

8 4,2102F-Ut 1.9644E go 1.92536 OO 1-859SE Q0 
_oe 2 OL 7d 2672 oo 26.12uU2F 00 2205.72 Es 00 
10 %,5149F-04 2.3689E 0 2.3003F JO 2.1975E 00 
ie. Z.27R86F-0i 2.60396 00 .2-5080F 00 £2. 390SE es 
12 Z,NGOGFAUL 2./7681E oo 2.6668E 00 2-5234E 00 
13 2.9244F.U1  2.9124E 90 2.7996E 00 2-6422E 00 
14 9.7549F-U4 Sei deve ie 5.n220F O00 2-8501E 00 
ee 29.5995F=00  Sp4587E a0 3.324B8EF 00 3.14398E 00 
16 9.493416-U1 3./806E 00 SASYIE 00 3.4441E 00 
177 ea SL Melle al 4.17172 oO. Fee do  3.-8096E 00 
18 oe Joes U1 4.6967E go 4.5206F 00 4.2910E 00 
See = (CU CBRE RIG = Ge OBSOE QO 5.1862F 00 4.9128E 00 
20 4.8442F-U4 6.2615E 00 6,9253F 00 5.700/7E 00 
21 4,71ni16-01 i 7,3919E QO 7.n999E 00 6.7102E 00 
22 eeyo eel 8". 3925 Esp &.5271F 00 8.9462E 00 
23 4. 4op 1b 204 1.G0909E 01  $$=|.‘1.9458F 01 9.8308E 00 
24 4.,S53R2E-U1 1.55 79e. 01 1.2958F 01 1.91756 0m 
2 ec | ee ei 1.61516 01 1.513S5E 01 
26 1 aS Ge U4 2.1 7eee ot 2e-n068F 01 1.-87S9E O1 
ae Aes Ea. 2.6233E 91 Poe h= Tel 2.3030E 01 
28 een eos Slip, $.252cem0n S.nS/8EF O1 2.8124E O1 
29 R.1965F-02 3.9520E 901 3.4932F 01  3.40919E 01 
30 7 Oe 2 4./797E gi 4.4400F 01 4.7 S6ue. 0s 
coy 6.4762F-02 j##$5.6769E 01 > Oo Segiee 0 ft 4.7649E 01 
32 5.6990E-U9 6.5794E gi 6.9195F O01 5.4377E Q1 
Sep  4,95a3E-092 7.97 he Gd 6.7090F O01 6.9102E 01 
34 4,2753F-02 8.LU958E ai /.2402F O01 6.4123E o1 
35 3.6499E-u2 8159316. 01 7.5681F U1 O.6199E 01 
36 3.067 0EF-02 8. SS aterm y. 7.A3/B8E O1 6.5887E 01 
i _PLezeeeele #8 8=5§=8. Sev teegt 7.42>30F U1 6.3005ERon 
| 9.0401E-U2 8.4918E 01 6.9214E U1 5.7580E 01 
Com .«i“‘éiétyt1«OAOUESUR 0 6=©=©6O7 BASE OL 6.45/74F U1 4.9971E 01 
40 1. 2090c=U? 6.8997E ol D1 eee OT 4Fo7LUc yw Ot 
41 9. U0 ea S 5. / 2g 22) ie 4.9814F 01 $.qn762E O01 
42 &.3244E-U3 —4.4241E 01 2.996 ier UT 2.n913SE 01 
43 4,U476E-U3 3.N975E gi 1 AS lees UT 1.292749E O12 
44 D oui S 1.8772E gi 1.900740 1 5.6906E 00 
45 1.0119F-U3 8.8666E 00 3,9564F U0 1.7245E 00 
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TABLE B.9 
SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 
ZER 3:1 CORE 
ir ENERGY MODERATOR © CLADDING FUEL 
See Usy SPECTRUM SPECTRUM 

1  eteaS0.9 SE 0g 5.5925E-901  $(5.4979E-01 i 5.4298E-01 

2 1.2311E O00 6.2790E-01 6.29299F-01 6.4154E-01 

3 9.6647E-01 8.2930E-a1 8.1943F-01 8.q045E-01 

6 7.7941F=01 1 .USS6E (90 1.92756 O00 1.9051E 00 

z, 6.4501E-04 1.2855E 0 1.2694E 00 1-2399E 00 

6 5.4731E-01 1.5295E 90 1.5062E 00 1.4668E 00 

Z 4,.7542E-01 $1.7671E 00 1.7333 00 1-6814E 00 

8 4.2192E-U1 1.996Se"00 1.9491E 00 1-8824E 00 

9 _ 388070 Bet 2. 248BE 0  .2eds0GEn00 2-9664E 00 
10 3.9119E-04 2.4208E 00 2.3379F 00 2-2334E 00 
ial 3.2786E=-01 2.6693E 00 2.56516 00 2.4384E 00 
12 3.0990E-01 2.8464E 0 2.7251E& 00 2.-5767E 00 
13 2.9244E.01 -3.0952E 90 _2.8671E 00 2-7059E 00 
14 2.7549E=01 3.26055 "aU 3$.10/77E 00 2.-9308E 00 
15 2.5995E-01 3.6043E 00. 3.4384E 00 3.9471E 00 
16 2.4311E-01 3.9724E 90 3.7952E 00 3.5910E 00 
17 2.2768E-04 4.4338E 900 4.24U04F 00 4.9184E 00 
18 OP. le7 seed, 5.07126 oo 4.8493E 00 4.5969E 00 
19 4.-9833E=-01— _5.9S326E 00 5.6662E 00 5.-3675E 00 
20 1.8442E-04 7.0673E 00 6.7408E 00 6.3797E 00 
24 4.7191E=-01  8.5809E 00 8.1684E 00 7-7201E 00 
22 1.9811E-01 1.0640E O11 1.9106F 01 9.5357E OO 
23 4.4571E=014 1.5451E ol (1.2740E 01 1.-4999E 04 
24 4.9382E=-04 1.72202 ¢1 1.6262— O01 1.S5281E 01 
25 4.2244E-04 2.2119 G1 2.0796E 01 1.9489E 014 
26 4.1156E-04  2,.8257E gl 2.6453F 01 —264701E O14 
2? 4.0119E=-04 —3.5816E ot $.33542E 01 3-1004E 01 
28 9.1325E-09 4.4981E o1 4.1626E 01 $.a540E 01 
29 8.1965E-02 5.5868E o1 SMISWCE 01 4.7266E 01 
30 7.5110E-02 6.8387E o1 6.2343 01 5.7119E 01 
St 6.4762E-02 8.20950€ 01 i i=($7.4117€ 01 6.7450E 01 
32 5.6920E-uU2 9.5986E o1 8.5787EF 01 7-7509E 01 
33 4.9583E-02  =1.0909E 02 9.6262EF 01 8.6235E 01 
34 4.2753E-02 —41.29036E 92 1.0458E& 02 9.2546E 01 
35 3.6429E-02 1.2889E 02 1.0966E 02 (9.5918E 01 
36 3.0610E-02 1.936426 Ge 1.1087E 02 9.5639E O41 
37 2.5298F-02 1.3434E 92 (1.0771E 02 9-1404E 01 
38 a 0401E-02 1.2997E 02 1.9010E O02 8.3273E o1 
39 .6191E-02 1.2954E 02 8.8597E 01 7-4902E 01 
40 .2396F-02  $$1.0644E o2 7.3949E 01 5-8138E 01 
41 x 1072E-03 8.8606E o1 5.7/754E 01 4.3530E 01 
42 &.3244E-U3 6.8496E 91 =4.4127E 01 2-9224E O1 
43 4.0476E-03 jj 4./7951E 01 2.5877E 01 1-6845E 01 
mt 2.2768E-03 2.9002E 01 1.35866 01 7.6744E 00 
45 1.0119E-03 1.3635E o1 5.2270E 00 2-278S5E 00 
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SPATIALLY AVERAGED THERMAL SPECTRA FOR THE CORNELL 
ZPR 4:1 CORE 























MODERAILOR 


SPE cnn y 


2 oS ee = 01 
6 Ree oue =o ft 


8.5241E-o1 


1-0430E 90 


1-2920E 00 


1.5386E 00 
1./795E o0 
2.0129E 0 
2.2350E 00 
2.4489E 0 
2.7049E 00 


2.8902E 00 


$.0592E 90 


3. S50 FE wo 
3$.6990E 00 
4.1049E o0 


4,.62965E 00 


5.J605E 00 


6.37 2a e 


7./365E 90 


_9.9931E 00 


Teel oven 


1 1-9 685E 01 


2.0456E gi 
2.6706E got 
3.49641E 1 
4.49483E oi 


55647 56) gil 


7.U749E qi 


8.7178E gi 
1.0516E 92 
1. 236°%E a2 
1. $1222 le 
1.5658E 2 
1.68449E 92 


— £.7557E 92 


1./7687E 92 
1.7157E g2 
1.9947E ge 
1.49104E 02 


i VS Sea ie 


9,0891E gl 
6.S611E 01 
J, C4SUE OL 


fis ENERGY 
‘i 1. CO SSE200 
2 4.2311£& UN 
3S = 9, OHK7E-04 
4 he 7941F- Ut 
2 6, 4 Sine = al 
6) 5. 4731F-u4 
7 — 4.7542E-91 
8 A eerie ud 
Lae, ae Sealy i 
10 Aco ti 9e Un 
41 %.27R6F-04 
12 %,0995F-01 
Dloeet(‘<‘<‘é éw« C«C PASE RO 
14 9./549E-U4 
mee ONDE -U 1. 
6 pe tlle sh 
17 9 27TASE=}01 
18 aoe 7 Ole tee 
a eae a 
20 ,8442E-01 
21 = - 71n1F-04 
22 Ae See Det 
23 fre O7 TES 
24 4,935R2F-U1 
2 1.2244E-01 
26 4.1156F-01 
eet—“<i«é‘wC«CL GA VE - 91 
28 9. 1S 5E so 
29 2.1965E-02 
30 7.5110E-02 
Sa 6.49762E-02 — 
52 5.6990E-u2 
33 4,95R3F.02 
34° 4.2753E-02 
35 %,60499E-02 | 
36 %,.0619F-02 
37 APS eam 
SS ——-9.04901E-02 
39 i. Sete =e 
40 ,2306F-U2 
44 i. 1N72E-03 
42 A. S2a4E203 
43 4,U476F-03 
44 2. 27A4BF-U3 
45 4.0119F-03 


ie one se an 4 


CLANUING 


SPECTRUM 


oy syipleey = = ji) al 
OF 2 000F Um 
Sie eA 
12H Suse. Ua 
57356 70 
1.51206 UO0 
1.7409F 00 
1.9569E U9 


2.16027F UD 


2 s55S5F UG 
2.5854F UN 
27 4559E 0 
2.90U0F UD 
Satizsoe Ut 
$.505949F Un 
S.ROODF Ut 
4.39576 UD 
5.n914F ud 
6.0455F UD 
7.3254F vO 
9 .V6ceewus 
12ers 
i. 47 2 Sie 
1.91441F O1 
2.486076 yi 
5.00 /7 Feee4 
4.qn918F ul 
el elles =i 
6.4047F Ui 
7.8221F U1 
9.,3594F U1 
1.08366 Ue 
1.?P192F 02 
e. 32/7 Sees 
1.3926F U2 
IWa4eSeaue 
1.36566 G2 
1. 9625E 502 
1.1117F 02 
vee el= (Cal 
7 4495 01 
D1 oOo G4 
5.44352F UL 
V.65e25F Ul 
CO. POC oe Ut 


roi 


SPECTRUM 


5.433S5E-01 
6.1235E-01 
6.n22uE-01 
eel (ii 
1 2S E D0 
1.4725E O00 
1.6889E 00 
1.-8913E 00 
2-978S5E Gu 
2.-248S5E QOD 
2-45/76E 90 
2-6007E Qy 
2.7369E QO 
2.97S57Ee OU 
$e 3105E Ou 
3$.6364E O00 
4.1556E OU 
4.8264E 00 
Oe FeO Ce eon 
G6.9351E OG 
8.55498E Qg 
4.n805E 01 
ies a70e 04 
1.79860E Q1 
2.-3299E O11 
2.-9955E 072 
S$.8049E Qt 
4.7761E O1 
5.R995E o1 
7.156/7E O1 
8.49497E O1 
9.7907E 01 
1.9922E 02 
1.1746E 02 
1. 973822 362 
1.2136E 02 
Tet o7 te coe 
1 qo eewn 2 
eeyedil iene: «(iil 
7.29485E G1 
5.3889E 01 
S$.5862E U1 
2-f459E O21 
9-92S57E6 00 
2.71835E 00 
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TABLE B.11 


CORNELL UNIVERSITY ZERO POWER REACTOR INTEGRAL PROPERTIES 


Integral 

Properties Water/Oxide Volume Ratio 

ee eo oom aga i a eo 
f 92185 .88785 .85480 .79331 .73687 
n .72300 1.72867 1.73165 1.73461 1.73605 
nf .58835 1.53480 1.48021 1.37607 1.27925 
oF 25149 1.30644 1.351200) 1e4717 sie ay 
57 .14087 1.14917. 1.15325 1.15677 1.15807 
D, cm. 3699 .3110 2770 2402 2206 
L, cm. 2.0086 1.9356 1.9253 1.9659 2.0298 
fs usec. 28.69 aye yi Pale \s 53.43 64.37 
vj» m/sec. 3667. 3318. 3134. 2945. 2849 . 
Vv.» m/sec. 3802. 3478. 3310. 3146. 3067. 
Ve, m/sec. 3971. 3626. 3446. 3265. Sul) - 
7, m/sec. 3802. 3426. 32038 3012. 2902. 
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while Table B.12 gives their temperature coefficients as 
defined by Eq. (7.5). 

The effective cross sections, Eqs. (6.23) and (6.31), 
are given at 20°C for the five cores in Table B.13, the 
temperature coefficients being given in Table B.14. 

Table B.15 lists the 20°C values of the mean cross 
sections, Eqs. (6.24) and (6.32), while Table B.16 lists 


the temperature coefficients of these cross sections. 
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TABLE B.12 


CORNELL UNIVERSITY ZERO POWER REACTOR INTEGRAL PROPERTY 
TEMPERATURE COEFFICIENTS 


Temperature 
Coefficients Water/Oxide Volume Ratio 
(107*°c7H) eee ee Don Bail 4:1 

£/£ +.266 +.546 +.813 +1.37 41.95 
n/n -.310  <.386 =.41 204 7cmen 0 
nf /n£ =.044 “4.1608 #2400 reo cmemiece 
on /on -2.76 =3.71 =<Qudeme=-5 S5ouemecneG 
5° / 8S -1.53 =l978 2 oe? 0c eet 
D/D +3.23 +4.40 “—F5 my Cmemer a o7 meron 
iL 45.51 46.89 +7.82 +9.10 +9.90 
b/2, +.048 =.411 S07 meeC eee 
Vv. Vin +8.32 410.1 $411.2 412.6 413.3 
v./¥. +7.76 49034) Evie mmEerIs ames es 
V¢/¥ 46.92 +8.55 49.58 +4+10.7  2+11.5 
v,/Tp 47.76 49.63 4 alOMQuer I > some sionO 
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TABLE B.13 
CORNELL UNIVERSITY ZERO POWER REACTOR EFFECTIVE CROSS SECTIONS 


Effective 
Cross Sections Water/Oxide Volume Ratio 
fone) ind ayo 2:1 oye Ac] 

a .0219 .0219 .0219 .0219 .0219 
2 .0139 .0139 .0139 .0139 .0139 
= 3641 3680 3700 3719  .3728 
oF .1967 .0195 9254 8295 7806 
ae 70 .1346 .1283 71220 .1190 
6685 6112 581i 5510 5364 
oe .9376 8927 2.8688 8445 ee: 
a8 niS7 2 1451 .1389 .1328 .1298 
ae 0182 9660 9385 9110 .8976 
ae 2186 .0414 .9473 .8513 .8025 
= .1609 1485 1421 1358 laze 
e 0327 9792 9511 9230  .9092 
Bie 2852 2618 2637 2655 =. 2664 
> 1584 1292 1095 0851 .0706 
30 4138 2.6804 8600 0800 3.2103 
. 9152 2.0965 2.2227 3807. 2.4755 
Pe 5722 2.8097 2.9694 1651 3.2809 
a 0816 0667 


- L036 


0482 


aa A TT I a LI a6 < a 


$372 
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CORNELL UNIVERSITY ZERO POWER REACTOR EFFECTIVE CROSS SECTION 
TEMPERATURE COEFFICIENTS 


Temperature 
Coefficients 


mice °c) 
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Water/Oxide Volume Ratio 
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Table B.14 (continued) 


Temperature 

Coefficients Water/Oxide Volume Ratio 
(107*°¢74) ie. 1.5:1 2:1 il 4:1 

o/s. +1.66 +2.19 +256 +3.05 Ae 37 


Sn /Se ane 55 +1.20 42.28 nee 
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TABLE B.15 
CORNELL UNIVERSITY ZERO POWER REACTOR MEAN CROSS SECTIONS 


Mean 
Cross Sections Water/Oxide Volume Ratio 

fen) ence ot 3:1 i 
a .0131 §©.0145 .0154 .0164 #.0169 
i .0080 .0088  .0092 .0097 +#.0099 
sf .2017 .2233. «2363S. 2506 ~~. 2582 
a 2.5177 226648 2.75600 se276G0bezeoL® 
se :0851 .0852  .0852 0 eSamuer oc 
a 3704  .3708 237100 meaydmmaria 
- 1.7624 1.9178 2.0141 2.1248 2.1867 
ae 0910 .0918 .0923 .0929 .0931 
Se 5641 .5860 .5992 .6138 .6215 
Se 2.5309 2.6794 2.7713 2.8768 2.9359 
a 0931 .0939 .0944 .0950 .0953 
=e 5721 .5940 .6073 .6219 .6296 
x 1430 .1588  .1684 .1789 .1844 
- .0917.  .0830 .0747 .0622 .0535 
a 1.3970 1.7214 1.9524 2 ocgmmoN sen 


). 1084 °° 1.3463 1.517450 soil orice 


Mi 


1g 


1.4687 1.8044" 2.0271 2. 3120me ao 


MI 
rp aD cto 


hUDo9 .0524 -0455 PU so2 T0Ze2 
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TABLE B.16 


CORNELL UNIVERSITY ZERO POWER REACTOR MEAN CROSS SECTION 
TEMPERATURE COEFFICIENTS 


Temperature 
Coefficients Water/Oxide Volume Ratio 
(107*°c7t) len Si 1 ya No. fl 
ola hee 0 Mey, Sia ALS 0 Gea 
=e By) re Boni -9.78 =oes -11.6 
B/S -8.13 -9.78  -10.9  -12.2  -12.9 
Se /he -6.16 EU /esiL i ~7.75 -8.51 eee 
ae pe +.16 +.14 HEL omelet ce) 
fe. a Vi Pers PaaG +.09 +.08 
ee 2 -6 .96 -8.36 “98 Sli). 2 -10.9 
Bee. - .60 -.76 = 2 12.08 =i oe 
o/s. Bean neuen Woe hk OD 
2a am -6.19 a7 bs =F) 7S -8.55 -8.98 
be / > - 54 275 -,90 -1.05 -1.15 
a /5t -2.79 -3.61 =< isms semesm 
e/a -8.46  -10.2 = Degeneres 
sh jsh 7.69  -8.97 =9083 ue medlo=cnemmeaninel 
a a a e e e r e 
sh sh Eno 7.19 -8.01 -8.93 -9-40 
S S oa e a a e e e 
sh 75h E736 Sli ee 2S S105" se ee 


[es a ae 
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—s Table B.16 (continued 


Temperature 

Coefficients 
(107*°c7) 1:1 

zh sh 

2 Oa -6.01 


sh ;sh 
el m— -7./76 


Water/Oxide Volume Ratio 


Leo Zon Beal ak 
-7.28 -8.08 ~8.97 -9 .44 
-8.87 -9 44 -9 .66 -9.70 





APPENDIX C 


NUMERICAL METHODS IN COUTH 


This appendix describes the computational methods 
employed by the COUTH code in solving Eq. (6.1) for the 
spatially averaged moderator flux 9 (E), as well as the 
subsequent calculation of all of the parameters defined 
in Chap. 6. A Fortran listing of COUTH appears in 
Appendix E utilizing the Nelkin water model for the 
scattering kernel. The input to COUTH is discussed in 
Appendix D. 

COUTH, as presently instrumented begins its cal- 
culations with the computation of the energy, velocity, 
and lethargy mesh as described in Sec. 7.1 and listed 
in Appendix B. Inspection of these meshes shows that 
the lethargy mesh has the finest relative spacing in 
the region of maximum interest, .025 eV < E < .6 eV, 
while both the velocity and energy mesh tend to have the 
finest spacing at low energy end of the mesh. For this 
reason all calculations in COUTH are done in lethargy, 


u, Space utilizing the relationships 


p(u) = Eg(E) , (ORL) 

o(u'—eu) = Eo(E'—~E) , (CED 
and 

aby) Be (5) =, (C.3) 


-150- 


a 








Appendix C Alps eik 


The heterogeneous effects of the lattice, Sec. 2.3, are 
incorporated into COUTH in terms of the moderator and fuel 
element escape probabilities. The escape probability of 
source neutrons, P(E); is calculated in the subprogram FUNCTION 
PSLOF utilizing diffusion theory as expressed by Eq. (3.4). 
P(E); the moderator escape probability of neutrons incident 
from the fuel element is calculated utilizing the reciprocity 
relationship, Eq. (3.6), by the subprogram FUNCTION PMODF. 

The linear extrapolation length into a black hole, dy: 


necessary for these calculations is approximated yee 


026229 


1M) — —_ 
d,(E)2,,.(E) = r= + 0.7104 , (C.4) 


7m 
I4+2 482, (EB) re 


where Ey is radius of the black hole. This expression equals 


m™m = ' 84/ m = 
4/3 at 2 ri, = 0, Zaretsky's value of 0.889 at 2 ry = 1, 
amd 0.7104 for TAs te approaching infinity. Leslie's 


SPECTROX mathode2! corresponds to the special case of A = 2/3. 


The fuel rod escape probability, P(E), is calculated 
from the parabolic flux approximation, Eq. (4.8). This com- 
putation is carried out by the subprogram FUNCTION PRODF. 

The subprogram SUBROUTINE PESCF calculates the fuel element 
escape probability, PAs utilizing the thin region transport 
theory embodied in Eq. (4.21). 

In the explicit calculation of escape probabilities, the 

modified Bessel functions of the first and second kind and of 


order 0, 1, and 2 are required. Calculation of the functions 





a 
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I, (2), I, (2), I,(2), K, (2); Kj (2), and K, (Z) is accomplished 
by the subprogram SUBROUTINE BESS. 
The modified Bessel functions of the first kind, I, (2), 


86/ 


are computed in the following manner. Noting that for 


any given z, 1 () approaches zero as n gets very large, we 


40 


set I.(z) equal to zero and I-1(2) equal to 1x10 “~ where 
for any z, n is determined as follows. For z < 45, n = 
[2z+10], where [x] is the largest integer contained in x. 
For z > 45, n is taken equal to 100. Lower order I's are 


then calculated utilizing the recurrence relationship®// 


T7(z) = 1 (z) + Lai (2) - (C.5) 


Recognizing that the I's calculated in this manner bear the 
correct relationship to each other, but lack normalization, 


the results are normalized through the expression 


oe) n 
e” = » I, (2) ~ I, (2) + 2 » Th) . (G26) 
=-00 m=1 


This same technique is utilized in the calculation of the 
- s utilized in the computation of the Nelkin water kernel, 
Hee (5.15). 


The zero order modified Bessel function of the second 


kind, K (2); is computed using the following powenaeicas 


83/ A. M. Weinberg and E. P. Wigner, The Physical Theory of 


Neutron Chain Reactors, p. 765, The University of Chicago 
Press, Chicago, 1958. 


BE yy Imee Zaretsky, PICG 5, 529°( L955)" 
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Hor 2 > 2: 


2 2 DD ee 7 
(2) a -* 11 wheat ale anal 
2! (8x) 3! (8x) 


3 


™ 
OB 








(C.7) 
the number of terms used in this expansion being set equal to 


mez. For z < 5, 


50 
K,(2) av iy 6 pom 62) : (Cae) 
m=0 
where 
5, = -(7 + log Z/ 26 (C.9) 
¥ = Euler's Constant — 90.57/21. (c.10) 
5, = 2, (C.11) 
and , 
6, = eo. m>2 (C.12) 


The first and second order modified Bessel functions of 


the second kind are calculated using the Wronskian relation- 


ship, 82/ 


ke 


K (2) Ty (2) + Koy (21, (2) == - (cmley, 


85/ D. C. Leslie, ''The SPECTROX Method for Thermal Spectra in 


Lattice Cells,'' AEEW-M211, United Kingdom Atomic Energy 
Authority (1962). 

86/ F. D. Federighi and D. T. Goldman, "KERNEL and PAM-Programs 
for Use in the Calculation of the Thermal Scattering 
Matrix for Chemically Bound Systems,'' KAPL-2225, Knolls 
Atomic Power Laboratory (1962). 
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The calculation of the slowing down source, S(E), as 
given in Eq. (7.1), and of the monatomic gas model for the 


scattering kernel, Eq. (5.9), requires the evaluation of 


the error function for various arguments. The error 


function, ERF(z), is defined by 


Z 
2 
enr(a) = 72. | é dee (G.14)) 
O 


In the COUTH code, ERF(z) is calculated by FUNCTION ERF 


in the following manner! For |z|>4, 


ERF(z) = 121 , 





(C.15) 
miele for 2 < |z| < 4, 
a9? 
ERF(z) ~ L2t eee 1-24 413-135 
l2| fr az2 2x4 = 3 
l 32507 
ean (C.16) 
eae, ) 
For |z| < 2, 
= En. GzP) 4 
-1 2(n-1l 0 
wD ZL -1 » VA 7, 
ERF(z) 2) =r [BNL O75) + BOTY G1) 707 : 
n=1 
(Cr. 


87/ H. B. Dwight, T 


ables of Integrals and Other Mathematical 
Data, p. 189, The Macmillan Company, New York, 1961. 
88/ M. Goldstein and R. M. Thaler, ''Recurrence Techniques for 
Calculation of Bessel Functions,'' Mathematical Tables 
and Other Aids to Computation, 13, 102 (1959). 
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Numerical integration in the lethargy space used in 
COUTH is carried out by the subprogram FUNCTION TRAP. The 
trapezoidal rule is utilized to simplify the calculations 
with the uneven lethargy spacing encountered. Integration 
of an arbitrary function of lethargy, f(u), which has been 
evaluated at the discrete points f(u,), i =P; 2)3e. ne 


is approximated by 


u n-L 
n 


/ eGo A plea +454) Au, . (Ca?) 
uy i=] 


The lethargy intervals, Au,,.are defined by 


Au; = U;4,-U;, L = 1,233.54. 0- 0 C6713) 


Calculation of the scattering kernel using the monatomic 
gas model, Eq. (5.4), is performed by the subprogram SUB- 
ROUTINE FKERN. Only the down scattering in energy (or up- 
scattering in lethargy) portion of the kernel is calculated 
in FKERN, the upscattering portion being calculated by 
detailed balance, Eq. (2.5), in the main program of COUTH. 
Calculation of the monatomic gas kernel in the form of a 
45 by 45 matrix as needed with the lethargy mesh used takes 
about 25 minutes of computer time. If the Nelkin water 
kernel is used, then the scattering kernel, Eq. (5.15), is 
calculated with the subprogram SUBROUTINE WATER. Again, 
only the downscattering portion of the matrix is computed, 


89/ M. Goldstein and R. M. Thaler, see Footnote 88. 
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the upscattering portion being calculated by detailed 
balance in the main body of COUTH. The actual methods 
used in the calculation of the Nelkin water kernel are 
discussed in Sec. 5.2. Since the calculation of 45 by 45 
matrix representing the Nelkin water kernel takes approxi- 
mately 15 minutes at a given temperature, the down- 
scattering portion of the matrix for a physical tempera- 
ture of 20°C has been reproduced on cards so that numerous 
lattice studies may be made at this temperature without 
necessitating the recomputation of the water kernel. Con- 
trol cards read into COUTH (see Appendix D), determine 
whether the downscattering portion of scattering kernel 
shall be calculated using the Nelkin model, or whether the 
Nelkin results for 20°C shall be read directly in. 

For calculations at temperatures above 20°C, the 
average cosine of the scattering angle for water, u(E), 


“ for every centigrade degree increase 


is decreased by 1x10 
in temperature above 20°C at all energies. The slowing 

down source of neutrons, S(E), is calculated according to 
Eq. (7.1), while the moderator scattering cross section, 
Ze(E), is calculated using Eq. (2.2) and as discussed in 
Sec. 7.1. The neutron balance equation, Eq. (6.1), 
rewritten below in terms of the lethargy mesh is then solved 


91/ 


using an iterative process .— 


20/ H. B. Dwight, p. 136, see Footnote 38/7. 


94/ M. J. Poole, M. S. Nelkin, and R. S. Stone, ''The Measure- 
ment and Theory of Reactor Spectra,'' Progress in Nuclear 


Energy, Series I, Vol. 2, p. 99, Pergamon Press, New York, 
London, Paris, and Los Angeles, 1958. 
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u 
Cc 
o— —T1) 
S(u; ) +f =(u,—>u,;) o Cu, )du, 


O 


ie Ghen )) 
S. (u,; ) Fe 2, (u,;) 1 + -Ft,) | 


The cutoff lethargy, Uo: is 7.372 corresponding to a low 
energy of .00101 eV with the thermal cutoff energy, Ew of 
1.609 eV being the zero lethargy point. The relative fuel 
absorption f/[1-f£] is calculated from the escape probabilities 
in accordance with Eq. (2.29). 


An initial guess is made that the moderator spectrum 


is a constant in lethargy: 
go (u;) = 1, AD én oot ahs 


(62.20) 
This guess is then normalized to satisfy neutron conservation 


as expressed in Eq. (2.4). 


94 (u,) 
$7 (u,) = = 1=1,2.34ne 


: m £(u;) —m 
2, (0) | Eo eee a; 97 (a; du; 


O 


(G2) 


The normalized flux, @7(u,) is then substituted into the right 


hand side of Eq. (C.19) and the next generation approximation 


ae 





Oo 





ae 


Oe 


_ 


Oo 
- a 
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to the moderator flux, 95(u,), is calcuiated. This process 


is repeated utilizing the general expressions 


u 


S(u,) | : Z(u;—= uso, (u, )d(u, ) 
Pye (Ug) = a, ii ln, 
st(u,) + 32 (u,) E + rte 
(C.22) 
and, 
Ou, ) ape ee ee 
[ : s™(u,) ; + ore | op (u, )du, (C.23) 


O 


The equation is considered solved when the successive 
iterations satisfy the following convergence criterion at 
every other lethargy point: 

am 
P47 Cu, ) 
th E pte ae Pe ee em ee! 50) 
1, Cu; ) 
(C.24) 
The convergence criterion, ¢, utilized in the COUTH cal- 
culations is ion In general solution of the integral 


equation, Eq. (C.19), has taken about 13 iterations with use 





ss 
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of the monatomic gas model, while the use of the Nelkin water 
kernel with its more complicated fine structure usually 
requires about 27 iterations for solution. 

The fuel and cladding spectra, the mean velocities, 
the disadvantage factors, the lifetime, the thermal utilization, 
nf, L, D, and the averaged cross sections are calculated in a 
straightforward manner as defined in Chap. 6. 

An entire lattice calculation, including the computation 
of the escape probabilities, the spatially averaged thermal 
spectra, and all of the parameters defined in Chap. 6 takes 
about 40 seconds of computer time with use of the Nelkin water 
kernel, and about 30 seconds with the monatomic gas model. 
Calculational times for the kernel are not included in these 
figures. Compiling time for the Fortran source program of 


COUTH runs about two and one half minutes. 








APPENDIX D 


INPUT-OUTPUT FORMAT FOR COUTH 


A Fortran statement of the COUTH code for use with the 
Nelkin water kernel as a scattering model is listed in 
Appendix E. This section describes the input required 
and the output received for a typical lattice calculation. 
The input data is listed below, the first column being the 
item or data card number, the second column is the Fortran 
FORMAT as it appears in COUTH, the third column is the 
Fortran list of the input variable, and the fourth column 


is a description of each variable. 


1. (14) IND Identification 
number. If IND =l1, 
the previously cal- 
culated down 
scattering portion 
of the Nelkin water 
kernel at 20°C will 
be used in the cal- 
culation. If IND > 2, 
the Nelkin water 
kernel will be cal- 
culated at the tempera- 
ture specified in 
Item 2. 


-160- 
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Zz. 


C2ERHOn 6515, E12.4) 


(13F6.1) 


(13F6.1) 


(7E11.3) 


(20F4.3) 


T 


RHO 


EPSI 


XA (T) 


XF(L) 


SK(J, I) 


AMUBAR(T) 


iRonk 


The moderator tempera- 
ture in degrees abso- 
lute. 

The moderator density 
at T in em/cc. 

The number of lethargy 
intervals. N = 44 gives 
mesh in Table B.1l. N 
must be < 54. 
Convergence criterion, 


€ Eq (G6. 24)— 


yu??? absorption cross 


section in barns. 


235 


U fission cross 


section in barns. 


Downscattering portion 
of Nelkin model 
scattering matrix for 


iA TOF 


Average cosine of the 
scattering angle in 


water, u(E, ). 
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i. COE W7eAcr ) 
o.. (8F8.4,E16.4) 


oe 


(10A8 ) 


CF(L) 


DEN 


ENR 


RO 


RCI 


RF 


SSC 


SAC 


FNC 


TITLE 


162 


Correction factors to 
the SK(J,1) diagonals 


see Sec. 7.1. 


Fuel rod density in 
em/cc. 

Fuel enrichment in 
weight per cent. 

Fuel element diameter 
in inches. 

Fuel cladding inner 
diameter in inches. 
Fuel rod diameter in 
inches. 

Cladding scatterer 
mass in amu. 

Cladding free atom 
cross section in barns. 
Cladding 2200 m/sec 
absorption cross 
section in barns. 
Cladding number density 


in atoms/barn-cm. 


80 character title. 
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me (2F10.5) WO The nominal water to 
oxide volume ratio 
for the lattice under 
consideration. 
AP The lattice pitch or 
rod to rod spacing (in.) 


in a hexagonal lattice. 


COUTH is instrumented to reread a series of items 10 until 
a negative value of Wd is encountered. If this occurs, 
items 8 and 9 will be reread with the subsequent read in of a 
series of items 10. Calculation will stop when the computer 
encounters a negative value of Wd, item 10, followed directly 
by a negative value for DEN, item 8. This allows one to cal- 
culate the scattering kernel at a given temperature, then con- 
sider a fuel element of certain specifications and calculate 
parameters for all of the lattices utilizing this fuel element. 
New fuel element specifications are then read in and the 
lattices utilizing this fuel element are analyzed. There is 
no limit on the number of fuel elements or the number of 
lattices for each fuel element type which may be considered. 
The output for a typical 1 lattice calculation using a 
45 point energy mesh is outlined below. The first column 
represents the page number after the execution of the source 


program, the second column gives the Fortran list of the output 





Appendix D 


164 


variable, and the third column gives a description of the 


output variable. 


“SJ ON U1 BS Co 


E(I) 
DE (1) 
V(I) 
DV(I) 
U(L) 
DU(T) 


BLANK 


eKCI;1) 


E(T) 
SX(I) 
SB(L) 
SS25(I) 
SS28(T) 
SS16(I) 
AMUBAR(I) 


Index identifying lethargy mesh points. 
Energy mesh point, E.. 

Energy interval, AE; = E..,-E, ; 
Velocity mesh point, Vs 
Velocity interval, Av,. 
Lethargy mesh point, u;. 


Lethargy interval, Au, - 


Scattering matrix, o(U(J)—~U(1)), 
J running horizontally and I running 


vertically. 


Energy mesh point, E,. 
Moderator scattering cross section in barns. 


Moderator transport cross section in barns. 
Ze) 


U scattering cross section in barns. 
PBN: ; : : 

U scattering cross section in barns. 

ot scattering cross section in barns. 


Average cosine of the scattering angle, 


u(E.), for water. 


|< 
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10. 


E(T) 
XA(T) 
SF(L) 
SIGMAF(I) 


FSIGMAF(T1) 


AA(T) 


SN(T) 


E(37) 


PR 
FI 
igi 
FO 
FIN 
SIGS 
SIGT 


PU 
PUNIF 


165 


Energy mesh point, E, . 
eae absorption cross section in barns. 
ue? fission cross section in barns. 


e e ° -1 
Fuel absorption cross section in cm “, 


£ 

=-(E;)- 

Fuel absorption cross section in ona 
f 

(Ey). 

1 + a(E;), a(E, ) = fuel capture to fission 

ratio: 

Slowing down source of neutrons, S(E), 


EG ema) = tie 


E(37) = .0253 eV. This page shows E = kT 
values for the variables used in the 
cladding calculations. 

CAB (CE IS 


Eq. (492208 


Fare Eqs. 


Jap 


Cladding thickness in cm. 


Solid angle factor, Eq. (4.26). 


(4.48). 
1 


Fraction of neutrons scattered in, Eq. 
Cladding scattering cross section in cm. 
Cladding total cross section in ga. 
e, Higa eco) a 
T, Eqs. (4u25)m 
ge cecrone 


Eq. (4.44). 


ae 


Pmif? 
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ED. 


Ee 


i. 


14. 


RSUM 
TAU 
RHOSUM 
CTE 


BLANK 


BLANK 


Wh 

Th 

GITLE 
E(T) 
FLUX (I) 
FC(L) 
FS(I) 
AB(T) 


wo 

T 
TITLE 
E(I) 
PS (TI) 
PM(T) 
PR(T) 


166 


ts dt Dae Isp (4) BIS 
%, Equa, (4a40)e 
Pan ++ Pee Eq. (4647) 


Cladding absorption rate to fuel element 


absorption rate ratio, Eq. (4.55). 


Input item 10. 
Input item 2. 
Input item 9. 
Bie iee ya. E;- 
@ (E;); Eq. (6.1). 
®,(E;), Eq. (6.4). 

o (E;), Eq. (6.3). 
5, (E) p14 3 -3 |: 
Input item 10. 


Input item 2. 


Input item 9. 


Energy, E;. 
P(E; ); Eq. (3.4). 
P(E.) hae eomeye 


PAE; ) > Eq. ((4Gy ee eene ye. 
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HD. 


PE(I) 
TU(I) 
RFA(T) 
RTE(T) 
FIO 


wo 
it 
TITLE 


ETA 
ETAF 
DIS 
DISC 
DIFFC 
DIFFL 
AL 
ALI 
VM 
VF 
VC 


FN25 
FN28 
FN16 
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P(E; ); ECiemce Ze 
f(E,), Eq.s(3n00e 

£/(1-£), Eq. (3.8). 
Rod abs. rate/element abs. rate, Eq. (4.53). 


J, (E,), Bq.) (4e220r 


Input item lO. 
Input item 2. 
Input item 9. 
Ey Eq oris). 


ae 

NE, Eq. 6Gr16)2 
m 

me HG (Oro. 
Cc 

o> Eq. (6.9). 


Ds Eq (ons co 
i, Eqv Cons 9)e 
LS Equrtomly) : 


‘a 
Ls L/nt. 
Var Eq: (6.6) 


Veo Eq. (6.6). 
Vo» Eq. (6 <6) 
Vi: Ean. (6..36ne 
y23° 


ASHE: 


number density. 
U number density. 


Oxygen in U0, number density. 
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luo 


FNC 
FNM 


EPSI 
RHO 
ENR 
AP 
R1 
RO 
VOLM 
VOLC 
VOLF 


TITLE 
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Cladding number density, input item 8. 
Moderator number density. 

Input item 2. 

Number of iterations required to solve 
Xa] Per Ie 

Input item 2. 

Input item 2. 

Fuel enrichment in atom per cent. 
Input item lO. 

Equivalent cell radius, cm. 

Fuel element radius, cm. 

Moderator area, om” 

Cladding area, ome. 

Fuel rod area, onan 

Actual water to oxide volume ratio. 


Input item 10. 


Page 16 lists the effective cross sections 
for each region of the cell, Eq. (6.23), 
and for the homogenized cell, Eq. (6.31); 
and the mean cross sections for each 
region, Eq. (6.24), and for the homogenized 


cell, Eq. (6.32). 
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Pages 12 through 16 are repeated for each lattice con- 
sidered, while consideration of a fuel element with new 
specifications will cause pages 10 and 11 to be printed 
with the subsequent printing of pages 12 through 16 for 


each lattice considered. 








APPENDIX E 
FORTRAN LISTING OF COUTH 


This appendix lists the Fortran statement of the COUTH 
source program as instrumented for calculation with the 
Nelkin water kernel as the scattering model. SUBROUTINE 
FKERN, which calculates the monatomic gas model for the 
scattering kernel, is also listed. For calculations with 
the monatomic gas model, SUBROUTINE WATER, FUNCTION ONE and 
FUNCTION TWO are replaced by SUBROUTINE FKERN. In addition 
the main program of COUTH must be modified to read in the 
values of the transport cross section of water as cal- 
culated with the Nelkin model for use in the calculation 


of the moderator escape probabilities. 


-170- 
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PROGRAM COUTH 

DIMENSION U(655)5F(55)5PS(55) sPM(55).,PFE (55) »RFA(55)2TU(55) »ABR(55) > 
1SN(55 5 °SK (55,55) +SX 655) eo ALPHA (55) «PHIN(55) spU(55)s AA(55) 5 FEA CGS 55 
EGAMMA (55) »shEL TASS) «SALPHA(5S5) «> SPHIN(¢5S5 js FEUX (55) 422055 ee 
3XA (55 )4%F (55 )9FS(55)55ST (55) 65 AMURAR (55) sKN(5)2IM(5)sCF(55) 2 
4S51GMAF(55),FSIGMAF 155 34NF (55) 5s 0V (55 ).V (55) »S525(55)+35528(55)s 
53S$16(595)45N(55)-.0N(55); TITLE (10) sPR(55 ) 5 FYO(55)+RTE(55) + FC(55)s 
6CTE(55)°.FB(55)-.TRB(5535 

ORDER OF DATA INPUT IS AS FOLLOWS. 


IND =----—--—--—---—-— mee ee ee 14 

WW ol Sdieahe. 2 N EPS! > SS ane P(F 10°64 154 2a elo 
2 ad na al elearning LSGr Seay 

i CSCS PPP StS aS ee ee D2 ec 1 

Sa Se ae ee TE oe? 

AM jE AR meet mee ere ee ee eee 20(F 403) 

Ce 9S == = = SS SS ae Phe 7 oe) 

DENsENR>ROs enc 8(F8.4)s4XsbI12.4 
ee Se eee 190¢(A8) 

HOPE HH SSS SSS Ses oe a nn ae ae PAVE eh te.) 


IND = ia NELKIN WATER KERNEL FOR 20 DEGe C READ INe 
PND = 2+ NELKIN WATER KERNEL CAMCUEATER Alen. 
T = MODERATOR TEMPERATURE IN DEGREES KELVINe 
RHO 1S THE MODERATOR DENSITY IN GM/CCe 
N = NUMBER OF LETHARGY INTERVALS. 
EPSI = CONVERGENCE CRITERION FOR THE ITERATION. 
XA IS HONECKS U235 ABSORPTION CROSS SECTION IN BARNSo 
XF IS HONECKS V235 FISSION CROSS SECTION IN BARNSo 
Sees. Gti SCART eER ING MAR a. 
Sk 1S SiGMatusJ) TO Utijj)FUNIT LETHARGY OR B(1)*Sieiiace Cae 
FCT))/SUNIT EV. SIGMA IS IW BARNS. 
AMUBAR IS THE AVERAGE COSINE OF THE SCATTERING ANGLE. 
= 15 THE CORRECTION FACTOR TO THE Shep ince neo. 
NEN 1S THE FUEL DENSITY IN GM/CCe 
ENR iS THe FUEL ENRICHMENT EXPRESSED AS WEIGHT PER CENTo 
RO IS THE #UEL ELEMENT DIAMETER IN INCHES. 
Rey %S fM4e FUEL CLACDING INNER DIAMETER IN INCHES. 
RF 15 THE FUITL ROD DIAMETER IN INCHESo 
AMc iS (tHE CLADDING SCATTERING MASS IN AMUe 
SSc IS THE CLADDING FREE ATOM SCATTERING CROSS SECTION IN BARNS. 
SAC 1S THE CLADDING 2200 M/SEC ARSORPTION CROSS SECTION IN BARNS. 
FNC 1S THE CLADDING CONCENTRATION IN ATOMS/RBARN-CM. 
WO TS THE NOMINAL WATER TOG URANIUM OXIDE VOLUME RATIO. 
AE 15 FHE LATT’CE FITCH IN Ne tee. 
mo PORMAT (113 

199 FORMAT i60HISCATTERING KERNEL» SIGMACUCJY) TO U1))/ZUNIT LETHARGY: 
fot =] Pras ots Us Pe.4n Vo ee. 

BOOP ORMAT €144Pxs10CEL1 63) 

203 FORMAT (3741NwATER TO URANIUM Oxfoe VOLUME RATIOsE Sc ls vpn oleereneeeere 
14H ABSOLUTE TEMPERATURE = F80294127H THERMAL UTILIZATION (F) = FS. 
P5e/9 TH ETA = FB659799H ETA*XF = FB55s/4133H MODERATOR DISADVANTAGE F 
3ACTOR = F8.5./7+32H CLADDING DISADVANTAGE FACTOR = F8.59/330H DIFFU 
4S1ON COFFFICIENT (CM) = F8.49s7s25H DIFFUSION LENGTH (CM) = FBRe4 s/s 
5254 LIFETIMES (MICROSFCONDS) :4210Xe11H THERMAL = F8.253/7,10X+,8H MEA 
6N = £8,507,;27H AVERAGE VELOCITIES (M/SFC)+7310X3:13H MODERATC. = ee 


- 





= 
se 
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To2s/s10Xs8H FUFL = FB02s/510X%e12H CLADDING = F8.224110Ks20H HOMOGE 
BNIZED CELL = F862) 

204 FORMAT (/77530H NUMBER OF ENERGY INTERVALS = 14./7524H NUMBER OF ITE 
LRATIONS = 1437%925H CONVERGENCE CRITERION = F1l1039s/73s329H MODERATOR D 
CENSITY (GM/CC) = F8.594s24H FUEL DENSITY (GM/CC) = F805s4334H FUEL 
3 ENRICHMENT (ATOM PERCENT) = F8e59/%9s26H LATTICE PITCH (INCHES) = F 
48.59//+11H RANDIT (CM)9/4510X%019H FQUTVALENT celle = F8.554510xKs16H F 
SUFL FLEMENT = F8.59/510X%012H FUFL ROD = FB.5s/44914H AREAS (CM**2), 
67+10X.13H MOpnFRATOR = €8.5./7;5;10x%;1l2u CLAD ING .— F8,5./3:10X.12H FUE 
7L ROD = FB.59//:46H ACTUAL WATER TO URANIUM Ox1IpEe VOLUME RATIO = F 
BBeSeZ/) 

206 FORMAT (4Xs31HEFFECTIVE CROSS SECTIONS (1/4CM)s29Xs26HMEAN CROSS SE 
ICTIONS €140M)5/54X%931 (1H) 929X%s2601H-)949268X%s LOHABSORPTION,42X) of 
222(612% 9 GHMONHNERATOR 4 F10,.5929X% 5492612 sG4HEUEL OF 15.5329X% )3492(12X+ 8H 
BCL ADDN INGsFl1.5329XK)9%92(8Xs 1LOHSCATTERING: 42X% ) 9 4752(12X%s9HMODERATOR, 
F10.5929XK ) 9/2 2(12X% 9 4HE UPR F195 651 29X 1 492(12X%,8HCLADDINGsF! less 29u) 
59732 (8Xs9HTRANSPORT 543K )9/492(412Xs9HMONERATOR.F1065:29X ) 9492612 94H 
6FUFL 6 F15.5529X 9479201 2X% s BHCLANDING 9 F111 65929X)9452(8X se SHTOTAL 147X) + 
7792(12X s9HMOHERATOR 4 F10,.5929X%) 94s 2012X% se GHEUFL 0 F1505329X 949261 2X98 
BHCLADDINGsF11.65929X%)9/22(8X%s VHF ISSIONs45X%) 9/4920612X%s4HFUEL + F15.25929 
9X)) 

207 FORMAT (7/752(8Xs31HHOMOGENIZED CELL CROSS SECTIONSs21X)s/492(12X510 
LHARSORPT 1ON9F9.5929X% ) a 49 2(12X%s LOHSCATTERINGS F965 329X ) 9 4582(12XKs9HTR 
CANSPORT 3F1055129XK) s/s 2612X% sSHTOT AL 9 F14,65329X%) ae%s2(1leXsTHFISSIONePF} 
32065329X)) 

208 FORMAT (7534H NUMBER DENSITIES (NUMRER/RBARN-CM)s/510X%38H U235 = Fil 
12649/910XK98H U238 = F12.49475310K010H OXYGEN = £12.42/4510K+112H CLADD 
PING = £120649/+810X0113H MODERATOR = Fil2e4) 

334 FORMAT (11(F704)) 


815 FORMAT (6H 193% e+6HENERGY 3 7Xs THDELTA Fs6Xs8HVFLOCITY:s2X%s 7HDELTA 
1 V»3X sBHLETHARGYs8H DELTA Us///5 (2X281492X%92(F1126352X%)s2(FB0192xK) 52 
ClF Pe432X)) 


1000 FORMAT (72¢¢1X%s1497(F1504))) 

1001 FORMAT (13(F60e01)) 

1604 FORMAT (37H1WATER TO URANIUM OxIprF VOLUME RATIO.,FS5.1+58H TO 12099 
124H ABSOLUTE TEMPERATURE = F60l) 

1905 FORMAT (2(F10¢5)) 

1006 FORMAT (2¢(F1026)115s2XKsE100P7) 

1007 FORMAT (8(F804)94X%9FE1 204) 

10908 FORMAT (10(A8)) 

1009 FORMAT (/7+10(A8)) 

PIOO FORMAT (5H1 I 94Xs6HENERGYs 7X9 1SHU235 ABSORPTIONs 2X%s12HU235 FISS] 
1ON.1X%+1S5HPUFL ABSORPTION. 1Xs12HFUFL FISSION. 5xXs9H1 + ALPHA.4xX%s12HS 
PLOWING DOWN. /:23X%s4¢(13HCROSS SFCTION.,2X%)2919X%s6HSOURCE 2/9 26Xs2(7H(B 
3ARNS )+8X%)92(6HC17CM) 39X%) 179 3903H~ )) 

1101 FORMAT (5H1 1 +4Xs6HENERGYs 7X3 14HH2O SCATTERINGs 2X3 13HH2O0 TRANSPO 
Pewee xs PeHUeLS5 SCATTERs 3Xs1 2hUaee SCAT bemeoxs | Tow 1G SCATTERs7xXs6HMyYU 
2 RARs/4123X%95(1BHCROSS SFCTIONs 2X) 9491 26X% 9506 7H (RARNS ) + 8X )9/7339(3H- 
3)) 

1102 FORMAT (BHO 1598X+6HENERGYs 9X s9SHMONERATOR, 6X+B8HCLADDING+s 7X+4HFUEL » 1 
13X+16HABSORPTION CROSS6/5:27X% 33 (B8HSPFCTRUMs7XK)s2X%s18HSECTION (1/7CM 
PH20) 945332 (3H- )144 910 1486X%94(F12 045 232X)9F1004) ) 

1103 FORMAT (BHO 1,8x,6HENERGY >. 20. 2HPS.8x »2ZHPM, 8x, 2ZHPR, 8x, 2HPE.8x,2HTU 
fe 7X9 SHREAL7X se BHRTF 6X s5HY OUTs/.:38(3BH- 93743 114,.6XsF leet ee eee 
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ele 
ZG0'0 


4444 
9198 
9199 
9298 


90 


23 


Fol 


92 


93 


Zee) ) ) 


FORMAT (14) 

FORMAT (7(E116%3)) 

FORMAT (20(F4e23)) 

mommMAT (5Xe4HI = +s 100 1 AR SX je7 wat I ) 
FORMAT (1492XsS(E1103)) 

FORMAT (S5Xs4HJ = 65§(12Ps39X%)e493H I> 
READ 1212+IND 

READ 1006sTs»sRHOsNsEPS!1 

NM=N+ 1 

poe 1010] «(XACT)» l=leNM) 

READ 1001s (XFCI)eIl=ieNM) 

READ 40004 ((SK(Je1)sJ=191)et2=15NM) 
READ 44445 (AMURAR(1I)oI=1«NM) 
READ 3344(CF(I)el=19NM) 

T=8 0616 7E-5¥T 

OM1=0e?FP05 

OM2=00481 

OMR=0 006 

FM=18 e 

EMR=2 232 

EMV =3e/llorm1loe/1lBe~1lo/2e32) 

mae >= 1] 020 

FA2Z8=8.3 

mo Oo- 3076 

AM?P52=2350117/1 200896 
AM28=2380125/1 200896 

AM16=162e0/1 «00896 

Bee 360 1 25/7235'8117 

ANU=2243 

START CALCULATION OF ENERGY. VELOCITY,LETHARGY MESH. 
ENDS WITH STATEMENT 94.6 
meoe | SLOGF(10< ) 

UC1)=0609 

byt) )=001*2200. 

mel )=4406 

wee) =V (1 )4+DV(1) 

DV (NM )=0.00 

DU(NM) =0-0 

DE (NM)=0.20 

BO 90 1=2.34 

Beal y=DVC1) 

VOTF1)=VCIT)40V01) 

BO St 1=35.N 

DV(C1)=DV(1)*EXPF (B*FLOATF CI-34)) 
VOT+F1)=VC1T)40DV01) 

DO 711 T=1.>6NM 

Biya C1) 

Siely=OV (1 ) 

DO 92 I=1+«NM 

V(T)=D(NM4F1—-1) 

mney C1 ye * 2) HSOP5 P9777 (2200s = 42) 
DOr set =) 4N 

DV(1I)=DDC(NM~I ) 


173 


CALCULATION 
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94 


224 
ae 2 


USS, 


ey) 


445 


ND hee 
io 


447 


448 


DO 94 Tf=25NM 

U( I )=LOGF(E())/E(1)) 
DUC I-1)=U(1)-U(I-1) 
DE (I-1)s6&(I-1)-F C1) 
Saint 19 

PRINT 815° (IlsECI1) sDECT)sVC1) oMV6CT) sUCT) DUCT)» T=14NM) 
PRINT 19 

DO 712 T=1,.5 
KJ=100%#FLOATF(1) 
IMC I) =KJ 

KN( I) =KU-9 

IF(IND~2) 22252235223 
CALL WATER(E.sTsaNe SK) 
DO 224 I=1+5NmM 


AMUBAR (CT )=AMUBAR CI )-6100*T7-— 2 2526 )/8s61637 
Bot 10O J=2,NM 
JIK=J=—1 


BO 110 T=1,JK 

gio (Ce (] )7E CI) YA¥2)*EXPR( CXS) )7 FE ) 
SK(J.1)=SK(1,J)*BB 

DO 951 L=1.4 

NI=KN(L) 

P= IMCL) 

PRINT 199%9TsNIoMI 

PRINT 91984 CIsI=NIeMI ) 

PRINT 200e( Ts (SKC Ue] )sJeNIoM] bs l= eNM) 

NI=41 

MI=s45 

PRINT 1995TsNI MI 

PRINT 92985 (I eI1=NIeMI ) 

PRINT 9199¢4(14(SK(CIe1)eIJ=NI oMI ) eo T219NM) 

DO 445 J=1.NM 

meets tT )=SK (Ia CrCl ) 

DO 119 J=1.NM 

WO 116 T=1.NM 

exnt(!] )=SK (eI ) 

SX(J)=TRAP(SXKA.DU.N) 

SX 1S THE SCATTERING CROSS SECTION OF HPO IN RARNS. 
DO 447 JT=1.NM 

SS25¢C1)=FA2Z5*SCROSS(E C1) + Ts AM25) 
SS?8(1)2FA2PB*SCROSS(E (1) sTsAM28) 

SS16( 1 ))=FA16*SCROSS(ECI)+TsAM16) 

$S$25328+16 ARE THE SCATTFRING CROSS SFCTIONS OF VU235.U238,016 
IN BARNSe 
ST(1)=0656*SQRTF (.0253/7E (1) )+SX(1)* 01 60 *AMUBAR(TI)) 
ST IS THE TRANSPORT CROSS SECTION OF H2O IN BARNS.) 
DO 448 lT=1.NM 

DO 448 J=1+NM 

SK (Je 1) =SK (Jo 1) *RHO* 060247718292 

= EN eye 

PRINT 10000 (yor (T)9SX(])9ST( 1) 9SS25 (1 )*+SS28(1)+SS16(1)sAMUBAR( 1), 


1l=1«NM) 


TE =T/EM4+OM1 7 (20 *EMy )+OM2 JEMy+OMR* (1 oO Z(EXPF COMRST )—1 0 )+05) /EMR 
TE 1S THE FFFECTIVE TEMPFRATURE « 
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Pans 


47 


ao | 


141 


DO 8 JT=1-e+Nm 

ST¢(1I)=STCI)*®RHO#.60247/18-02 

SX(1)=SX(I )*RHO*#.60247/18.02 

SK ,sST»+AND SX ARE NOw IN (1/7CM).o 
AlL=SQRTF(E(I)/TE) 

SN( J )=E(1)*ERF (Al ) 

BE TA=TRAP(SN,DUwN) 

DO 9 IT=1.6NM 

SN(I)=SNC1I)/BETA 

SN IS THE NORMALIZED SOURCE 

T=T/8c°6167E-5 

READ 1007 «eNENeEFNReROSRCI »+RF es AMC+1SSC+SACeFNC 
Ma ePeEN) 41.41,.715 

READ 1008eTITLE 

AMC=AMC/ 1200896 

RO=RO*2 054/20 

RF=RF#2054/2 6 

RCI=HERCI*¥2054/2 6 

ENR=1000*ENR*¥RM/ (ENR¥RM+(1000-ENR) ) 
FN2PS=ENR¥* 060274 7#DEN¥*1 eo OE-2/7P70026 
FN28=FN25*(1000-ENR)I/ZENR 

FN16=2001* e-60247*DEN/270026 

FNM=se 6024 7*RHO/18°02 

DO 713 T=14NM 

SIGMAF (CT )=FN25*KA CI )+FN28*2071*SQRTF (20253/F(1)) 
FSIGMAF (1 )=FN2S5*xXF (1) 

PR(1)=PRODF (FN25.,FN28,.SIGMAF (1 )sFSIGMAF (I )+«DFNsREF»«RCI«RO, 


1$$25(1)+2eSS28(1)+SS16(1)) 


meee l )=SIGMAF CI )/FSIGMAF (I ) 
PRINT 11200 
PRINT LOO0s (yeF (IL )sXACT)+eXF CT) *SIGMAF CI y»sFSIGMAF CI) +s AAT) +SN(T)-> 


f= 1 ao NM) 


CALL PESCF (EseFNC+SSCo SACs AMC sRCI eROseNMeTePReFIO.PE SCTE ) 
DO 47 T=1+NM 

ee ot )=F JOC 1 )*(te-PRCT) SC PROI (1 Pe) 

RTF IS THE RATIO OF THE NUMBER OF NEUTRONS ABSORBED IN THE FUEL 
ROD TO THE NUMBER ABSORBED IN THE FUEL ELEMENT .e 
READ 10055Ww0O,AP 

TF (WO) 71437145141 

PRINT 19 

PRINT 1004,WO+T 

PRINT 10096TITLE 

R1y=SQRTF (CCC AP*¥ 2054) ¥**2)¥16732057 (2 0*¥ 32014159) ) 
Y=R1i/RO 

VOLM=3014159%( (R1*#*2)—(ROF*SD) ) 

VOLFH=a3Be14159* (RF¥%2) 
VOLCH3014159% ( (ROX*#P)—(RCIF*F2) ) 

VOGEL E=ae141594(R1*42) 

R=VOLM/VOLF 

RC=VOLM/SVOLC 

Pee yl = 1 «NM 

PS(IT)=PSLOF (ECT) sROsYsRHOsSTC(I)) 
PM(T)=PMODF (ECT) »®RO+YeRHOsSTC(!I)) 

TUCT) aC re~PECI)IFPSCI))ZC1e-~PM( I) *PE(T)) 
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a 


le. 


14 


16 
as 


27 


2B 


Eo 


2 2 


35 


Si 
60 


36 
Sie 
70 


v1 


we 


TU IS THE RATIO OF THE NUMBER OF NEUTRONS ABSORRFD IN THE FUEL 
FLEMENT TO THE NUMBER ABSORBED IN THE CFLLo 

aoe = | 1761 5—hU Giey 

RFA 1S THE RATIO OF THE NUMRER OF NEUTRONS ABSORBED IN THE FUEL 
ELEMENT TO THE NUMBER ABSORBED IN THE MODERATOR. 

FEACT Y=RFACI )»*® o656*RHO*¥,60247* SQRTF ( .0253B/4E(1))418e02 

FEA IS THE FUEL ELEMENT ABSORPTION CROSS SECTION IN 1/7CMo 

AB CT) =(C1e¢+RFACT))*0.656*RHO*, 6024 7#SQRTF (.02534E(1))718002 

AB 1S THE TOTAL ABSORPTION CROSS SECTION IN THE CELL PEFR CM H20. 
BEGIN ITERATIVE SOLUTION OF INTEGRAL FQUATION. SOLUTION ENDS ON 36 
DO tie T=1.«NmM 

PHIN( I )=12e000 

INITIAL GUESS AT SPECTRUM IS A CONSTANT 


<=] 
Woet3 [.[=1+6NM 
REePAACI)=PHINCI)#ABCT ) 


B= TRAP (ALPHA,DUsN) 

DO 23 IT=1+NM 

PHIN(I)=PHINCI)/B 

DO 15 IT=1+NM 

DO 16 J=lsNM 

GAMMA (J)=PHIN(J)*SK (4,1) 
DELTACI)=TRAP(GAMMAsDUsN) 

DO 27 fT=I16NM 

SPHIN( IT )= (SNC T)4DFL TACT ))ISZ(CARCI)+SX(1)) 
DO 28 IT=15NM 
SALPHA(1)=SPHIN(I)*AB(1 ) 

C=TRAP (SALPHAsDUsN) 

DO 29 JT=1+NM 

SPHIN(1I)=SPHINCI)ZC 

DO 32 IT=1lsNMee 

RA=SPHIN(1)/PHINCI) 
ABSDIF=ABSF (RA]=}1 0 ) 

IF (ABSDIF-EPSI!) 32:32:38 

een PNUE 

GO TO 36 

=k } 

IF (K~-100) 39+39436 

DO 60 J=14NM 

PHIN(J)=SPHIN(J) 

GO TO 14 

ie, ay, T=15NM 

FLUX CI)=SPHINCIDZECI) 

Oo FO [=ile«eNM 

ZZ(1)=FEAC(I )¥RTECI)*SPHINCI ) 
C1=eTRAP (2Z.+DUsN) 

DO 71 IT=1>+NM 

ZZ(1)=AB(I)*SPHINCT) 

CP=aTRAP (2Z2Zs0UIN) 

R=Cl7 C2 

DO 72 IT=1sNM 
FS(1)=SPHIN( 1) *R*¥FEA CT) *RTE CLT )ZSIGMAF (1 ) 
FC(L)=FEACLT)*CTF(CI)*SPHINCT) *RO*SQRTEF(F (1) )7(SQRTE (20253) *FNC*SAC) 
ZZ(1)=FEACT )¥RTE CE) *SPHINGI )ZAACT ) 





Aas 
- ~o Fg) Arr ees 
: 
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gs: 


74 


a7 


{oe 


ee: 


80 


81 


82 


a3 


— 


C3=TRAP(ZZ.DUsN) 

E TAF =ANU*¥C3/C2 

ETA=ETAF/F 

BG Fst SSI, 

ZZ(1)=SPHINC(I)/VC61) 

C4=TRAP (2ZZ+DUsN) 

DO 741 =16NM 

ZZ(1)=FSCI)/VCT) 

C5=TRAP(ZZ.DUsN) 

DO 75 I=1+NM 

MPU) yer CClL)/V (61> 

C6=TRAP (ZZsDUsN) 

Des=Cc4/CS 

DISC=C6/C5 

C7=TRAP(SPHIN+sDUsN) 

CB=TRAP(FSsDUIN) 

CO=TRAP(FC.DUsN) 

VM=C7/C4 

VF=C8/C5 

VC=C9/C6 

DO 77 Y=1+NmM 

ZZ( 1 )=SPHIN(I )*e656*RHO* .60247*SQRTF (.02537F(1))718202 
C1O0=TRAP(ZZsDU.N) 

DO 78 I=l1+NM 

ZZ(1)=FS(1)*SIGMAF (I) 

C11=TRAP(ZZsDU.wN) 

Boe 79 1=1>.N™ 

ZZC1)=FCCI )*SACKFNC*SQRTF (.02P53/E(1)) 
C12=TRAP(ZZ+DUIN) 

Zio. BO =] «NM 

ZZ(1)=SPHIN(I)*¥STC(I) 

Gi 3=TRAP(ZZ+DU.N) 

DO 81 JT=1+NM 

ZZ(1)=SPHINC(I)¥SXCT) 

C14=TRAP(ZZsDUeN) 

DO 82 T=146«NM 
ZZ(1)=FS(1)*(SS25 (1) *FN25+4+SS28 (1 »)*FN284+S816(1)*FN16) 
C15=TRAP (ZZsDUeN) 

DO 83 IT=1+NM 

ZZ(1)=FS(01)* (S825 0] ) FEN2S* (1 0-207 (36 *AM2S ) 9+59S28 (01 ) * EF N2B*¥(1.-2e7 
(36*AM28) )4+SS1601 )*FFNI6B*¥® (1 o-Po%(R2e*AM1I6))) 
C16=TRAP(2ZZ+sDUsN) 

T=T*8e6167E—-5 

DO B86 T=1eNM 

FR(1)=(SPHINCI )*VOLM4+FS(1)*VOLF4+FC(1)*VOLC)/VOLCELL 
TRACT) =CSTCT)*®SPHINGT ) *VOLM4 (FS 01)*S1IGMAF(1)4+Z2Z¢01))*VOLFF(F CCI) ¥¢ ( 


LSAC*ENCESGRTE ( .0253/7F 61)))4+CFNC*SSCFSCROSS (FCI eo Te AMC )* (107204 (30* 


ZAME ) ey) 17 CVOLCELESFB (1) 


S622 Gl =ee<')/ TRACI) 


CPO=TRAP (ZZ.DUsN) 
C21=TRAP (FBsDUsN) 
DIFF GC=C207 ( set Cr!} ) 
DO 84 T=1¢eNM 


R4 ZZ¢1T)H2FCC(I )Y*ENC#SSC#FSCROSS(FE CT)» Ts AMC) 
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C17=TRAP(ZZsDUeN) 

Va Aslven lee SS 
CIB=C1IT7*(Clomee/(30%AMC) ) 

DO 85S [=1+4«NM 

YATE Mle Sl ees eM iy et) 

C19=TRAP (ZZ+DU.N) 

Sr NM= (117 C7 

SBAF=C11/C8 

SBAC=C12/C9 

SBSM=C14/C7 

Seon -C1S7Ca 

Sooc=-C17/C9 

SeeM=C137C7 

SBTRF=(C11+C16)/CB 
SBTRC=(C12+C18)/C9 
SBTM=SBAM+SBSM 

SBTF=SBAF+SBSF 

SBTC=SBAC+SBSC 

Se-F =C19/7CB 

VO=220060 

SHAM=SBAM*¥VM/VO 

SHAF =SBAF*XVF/VO 

SHAC=SBAC#¥YVC/VO 

SHSM=SBSM*¥yVM VO 
SHSF=SBSFFRVFE /VO 

Siro > oBSCeVC/VO 
SHTRM=SBTRM*VM/VO 

SHTRF =SBTRFEF*VF sVO 
SereCc=SBERCAYCSVO 
SHTM=SBTM*¥VMSVO 
SHTF=SBTFRVFSVO 

SnCc=ss)rCeVCVO 
SHFF=SBFFRXVF/VO 

DZ=C7* VOLM+C8B*VOLF+COXVOLC 
SBA=(C1O0*# VOLM+C11*VOLF+C12*®VOLC)/NZ 
SBS=(C14*#VOLM+C1S5*VOLF+C1 7*VOLC)/DZ 
SBTR=(C13*VOLM4+(C114+C16)*VOLF+ (C124+C18)*VOLC)/DZ 
SBT=SBA+SBS 

SBF=C19*VOLF/DZ 

HV=DZ/ (C4*¥VOLM+CS5¥FVOLF+CE6*VOLC ) 
SHA=SBA*HV/VO 

SHS=SBS#HV/VO 

SHTR=SBTRFHV/VO 

SHT=SBT¥HV/VO 

SHF =SBF*#HV/VO 

DIFFL=SQRTF (DIFFC/SBA ) 
AL=((C4*¥VOLM+CS5*¥VOLF+C6*VOLC )*¥1 00E+4)/(C2*VOLM ) 
AL 1=AL/ETAF 

DO 76 T=leNM 

reChi=-e Cl )/ECT ) 

PS Cry =rPSC¢1) EC! ) 

PRINT 1102 ¢(ToFCLT)sFLUX(IT)seFCCIL)sFSCTtsAR(CLT)eT=al1eNMy 
PRINT 1004,.,WO+T 

PRINT 1009,.,TITLE 


178 
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Appendix E 179 


41 


PRINT 11030 ¢€1oF(P) sFPS Cy VoPM C1 Pe je CG Ty re es 


1FYOCT)sf=1eNM) 


PRINT 203sWOsTsFeo ETASETAFSDISs DISC »DIFFCODIFFEs AL s AL1s+VMsVFeVCsHV 
PRINT 208s+FN25+FN28e0FNI 65FNC oF NM 

PRINT 204sNeKsEPSIT +> RHO.DEN,ENR,AP,R] »RORE»VOLM,yOLC,VOLF,R 

PRINT 1009,TITLE 

PRINT 1004.WO¢eT 

PRINT 206>SHAMs SRBAMs SHAF + SBAF s SHACs SRAC*+SHSMsSRSM+s SHSF » SRSF4sSHSC » 


1SBSC+SHTRM,.SRTRM. SHTRE» SRBTRE sSHTREsSRTROs SHTMs SBTM. SHTFs SBTFsSHTCs 
PSBTC+SHFF » SBFF 


PRINT 2079SHAsSRAsSHS+SRSe SHTReSRTRs SHT se SRT so SHE » SRF 
PRINT 1009,TITLE 

PRINT 19 

Geom TO: 301 

CONT INUE 

END 


FUNCTION PROF (FN25«¢FN28s+SIGMAAsSIGMAF eDENeRE «ROT sRO.SS255SS28 » 


S16 ? 


THIS CALCULATES THE ESCAPE PROBABILITY FROM THE FUEL ROD USING 
THE PARABOLIC FLUX APPROXIMATION, SFE KeBe CANDY THESIS+ PG 118.6 
FN16=2001* e60247*DEN/S270 026 
SS=SS25*FN25+SSPRXFN2R+SS16%FN16 

SA=SIGMAA 

ST=SS+SA 

SF =S1GMAF 

C=SS/ST 

A=ST*RF 

B=1e-C 

IF (A-05) lels2 

CALL BESS(AsFIOsFIl1lsFI2sFKOeFK1 os FRO ) 
FILO=A*(FI1*FKO4+10-A¥#(FI1*FK1+FITOXFKO) ) 
FLL=(FII*FK14+16~2e*F10)/30 

FMOO=1e-2o0%*F 11 

PR=1l0-2o*XA*B*¥( }o-FMOO)/(10e-C#F MOO) 

SO. TO =3 

CALL BESS(A+sFIOseFI1sFI2eFKOsFKI se FKZ) 

FIO=A*¥(FIIL*XFKO+1 o-AX(FII*FKL+F LOXFKO) ) 

FILI=(FIIL#FK1+1 em~2e%*F10)/30 

FEO=(FI2*FKO—1 ot20%F10+20%F11)/30 

FPL=( (CF le*FK1/A)4+F11-F20)750 

Fee= (CCF IARFK AS (CAF*2) 1401 oA (3B e* (AXED) ) )—-DOo*F21 SAT 

FMOO=1 e-2o*F 1 1 

FMOP=1e-Go*¥F11+80%xF 21 
FMP2PH=10e~GeX*F114+240%F 2P14+2P40X%F CP—-Pe/(Ak*2D) 
FMO=Po0*¥A*(1e¢—-F MOO ) 

FMe=z=AX* (1e-FMO2 ) 
G=10-C¥(4e*®FMO0+4G o X¥FMP2—-G0 *FMO2) +(C¥¥2)¥ (Ge XEFMOOKFM22P—-306%* 
(FMO2**?) ) 

PR=1.-B* (FMO-C¥ (FMO¥ (4 6 FFM 22-30 *FMOD) +60 *FMe* (FMQOO-FMO2)))/G 


SP ROE =(la sr ?/7RC!I )+RF*PR/RC!I 


END 


Gy Cy CO} 


— 
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SURROUT INE PESCF (Es FNCs+SSC»SACs AMC oRE eROSNM,TsPReFYO,PE.CTE) 
DIMENSION £€(55)«PR(55 ) +FIJ0(55) »PE(55) +CTE(55) 
THIS CALCULATES THE ESCAPE PROBABILTY FROM THE FUEL ELEMENT GIVEN 


THE ESCAPE PROBABILITY FROM THE FUEL ROD. 
IS USED AS DEVELOPED IN KeBo GADY THEs hs. 
19 FORMAT (1H1) 
SO FORMAT (10H ENERGY = Fl264s/e31H FUFL ROD 
1597s319H ROD CURRENT OUT = F8.59/75335H FUEL 
2TY = FBeD59/%s27TH CLADDING THICKNESS 
3 FRACTION SCATTERED IN = F8e5s/4920H SIGMA 
SsiGMA TOTAL = FS8.59/15H G = FEsSs7 sone 
5 PUNITF = FB.5s/+14H RIN + ROUT = 
6+ RHOOUT = F8R0594s7H CTE = FBS) 


Bee 2> 14159 

GAM=e5772157 

Pre 2616 7E=5 

Zee JRO 

TT=RO~-RF 

p= 1o— (2ZExee) 
mO-2oeAS TNE (ZE)/P I 

FIl=ZE 

Pe- OTe o®ZEaSORIF (Fis) 7P I 
Beet ee Ze Mo—( ZEX eT Go ee 
OUP Ix (ZEeHe ) 746 
FLAM=LOGF (46/SQRTF(FK) ) 


THIN REGION THEORY 
APPEND oo > 


ESCAPE PROBABILITY = F8e 
ELEMENT ESCAPE PROBABIL 
FB659/95H FO =F80D94%925H 
SCATTERING = F8e5e/315H 
F8,59/4%,6H PU = F865e/7e9H 


FB e5e4e7H TAU = F8B6594918H RHOIN 


D2=(2 0/30 )* (100 7O*F K-So*% (FL AM=6 2S JX (ERKEX2)7 1 Gir — ao © (REAM ee 


1126) )¥(FK¥¥3)/640) 


FLis(F2-4o*D1/P1)/F 1 
FL2=( Be *F3/3e-2e¥FK#F1—-4o*D2)/F 1 


Cl=(PIXFO/20+ZE*SQRTF (FK) )/(20*ZE ) 


S2-De/(ZE**2 ) 
mini 2—=4>o~rCis7Pl—-ZeE 


Bie e=PeXFR+Ao* (ZE*% 2 )/3e-40k Ze ae? 


PmeeN=AS] Nr (eat ZE/ (1 otZE) YVR 
pow: O [=i sNM 
Sies=SsSC=FPNC*#SCROSS(ETI )«e Ts AME) 


SIGT=SIGSt+tFNC*®SAC*SQRTF (.00253/E(1)) 


SIGA=SIGT-SIGS 


Sahl sol Gi*RO*F 1 *F Li eS" ( (SIGIR jeer 2 jie 2 


m= StoTe (7 


PU=oS*xX* (GAM+LOGF (X 7-1 05)+10—-(X¥**D) Get (X**3) 4B oe — (X**4) S360 ot 


(X**5)/2880 6 


BUNT eS=ew7 (i e—-( SIGS (ieee) 47s nce 


Mx se O55) Sis alvse 


G= (Le ~F1)7~2o*¥ (1 o~F2)¥*¥SIGT#ROFG 0 *¥ (1 a -F 3) ¥((SIGT*RO) ¥*2) 73.6 


COTO 33 

G=-0027?1 P-. OP 66*LOGF (X )/S (ZEXRS 2) 
CONTINUE 
eoetlaolGS*( 1] o0—-G=-S)*PUNTF7s1 Ca 


TAU=1 o~SIGT#ROXFLL1+65% ((SIGT#RO )*#*2) *FLLE 


RHeswM=oiGst*( 1 >5>—-TAU)*PUNTEZS Tem 
RIN=FIN*®RSUM 

POUR oUM—R ITN 

RHOIN=F IN*RHOSUM 
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RHOOUT =RHOSUM-—RHOIN 

FJO(1)=PRCI)¥(S+tRIN)/(1 e~RHOIN*PR(1)) 

PE (1)=G+ROUT+FJ0(1)* (TAU+RHOOUT ) 

CTE )=HESIGA* ((16-G-S)FEYOCT 7% (1 TAU Ole SG e oes. 


te eye) 


a loo 10440510 


40 PRINT 19 


AG, 


86 


65 
S7 


88 


PRINT SOsE (C1) +PRC1L) eFYO(1) sPEC1)<TTsFOsFINS SIGS sol Gives oe sau hee 


1RSUM, TAUsRHOSUM.sCTE (TI) 


— 


ent) 19 
CONTINUE 
i176°0616/7E-5 
= INO) 


FUNCTION ERF(Z) 

ERF(Z) CALCULATES THE ERROR FUNCTION, 
eA Se OZ) 

IF (AZ-4e) 85386386 

FRF=AZ/Z 

RETURN 

IF (AZ~2oe ) 87388: 88 

i= le 

T2=T1*(Z*¥*¥2)/30 
T3=T2% (ZEK2 KB 0/ (2050 ) 
T4=T3% (ZERO )K5 0/21 0 

TS=T4* (Z¥*¥2)K7 0/360 
T6=TS* (ZEEE YEG 0/55 0 
T7=T6*(Z¥*2)*¥11 0/78 0 
TR=T7¥(Z*¥¥2)%¥1 3041050 

T9O=TBE(Z¥EX2 )¥K150/1360 

Mr O= 19% (Z**2 )*¥1 72/171 0 

T11L=T1O0% (Z¥¥2)¥#¥190/2100 

Pve=T1ii* (Z2¥*2)*218/ 2536 

Mpsstie* (24S )*¥23073006 

faa 1 ae (ZkEP ) KOS oS 351 6 
TIS5=7T14* (ZX*#2 )¥P7 o/S4068 

T16=T15¥* (Z¥*XD)¥290/465 

TY7=1 LOX (ZH*D 1*¥31 0/528 oc 
T1R=T17¥ (Z¥*2 )¥33 6/595 0 

T19=T18% (Z¥¥2)*¥35 0/6666 
T20=H=T1IO*¥(ZE¥D)H3B7 0/741 © 

mop 1720*(2*s>)*195.57G20-> 

ERF=26* 956418968 2¥ (T1-T24+T3-T44T5—-T64T 7-T84+7T9-T104T11-T12+T13 
=H) 15-TIGtTi7-~TISt119-lee77a2 
RETURN 

ZA=2Pot* (Z¥*2) 

Slaia 

Sea te7 7h 

S3=3B0/(ZA*XK2) 

S4=15¢/(ZA¥*¥3) 

S5=52 05/7 (ZA¥KG ) 
FRE=(10-25641896* (S1-S24+837 S44+85)/(AZ*FEXPF (Z*¥*¥2)))*®(AZ/Z) 
END 





4 


Ul 


52 
110 


si 
ee 


--- 
=a ps 
i> 


eS 


Lie. 


Lag 
ES 


ier 2 


140 


SURROUTINE BESS(ZsC+DeEs&sGoH) 
BESS CALCULATES MODIFIED BESSEL FUNCTIONS OF THE FIRST AND 


SEGOND KIND s 


C=10(2Z)sDET1 (Z)s0F aI 212) 6FH50.02Z 1s Gah eZ aes 2 Zz 
DIMENSION 8(100).BK(3) 


Wet z-456) Lele2 
N=XFIXF (C2 6%Z4+100) 
Gor To 3 

NM=100 

SN 1 )=1eOE—40 
DO 8 I=N4100 
Bt(I)=000 

SUM=2 ¢OE-~40 
No=N — ] 

DO 4 L=e2.NJ 
== |o 


BONE )=PLOATE (NL ) *20"%8 (NE ie Ze Es 2 


SUM=SUM+20%*B (NL) 
SUM=SUM-B(1) 
FACTOR=(EXPF(Z))7SUM 
POMS L=al«sN 
BC(L)I=S=FACTORZB(L) 


Mec (lL) =) 3 OE—307) 65655 


B(L)y=OeO 

CONT I NUE 

C=B(1) 

D=B(2) 

F=B(3) 

Meer —Sa) 81+82.82 
ieree2—-Ge) 11O¢elii«111 
<KL=10 

G@ 7O 140 

Bee Z—-f/o) VIAsli3sel113 
mE=12 

GO TO 140 

Heer eed) Li4e11S5e115 
mi 14 

Sono | |;46 

IF (Z-9e) 116091179117 
re 16 

Go TO 140 


IF (Z-100) 11891199119 


Se-16 

GO TO 140 
Rik=10 

GO io 140 
w=1 00 
SUM2=10¢0 

DO 83 IT=1+Ku 


We-W*((Pe*®¥FLOATF(I)~10)*#*2)/ (0 °(Z*¥B. )*FLOATEF (I) ) 


SUM2=SUMet+W 
SUM2=SUM2-W72 o 


RIK C1 )= (SUMS*SQRTF (3214159276547 (2e*Z)))/EXPF(Z) 


SOO. 2 


182 
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rom a7 


OY Oe 


81 


~ 


SUM1=0 00 

DO 7 123399992 

SU 1=SUM)I+B.0 1) /PLOA EEG) 
SUM1=4e*SUM1 

BR (1) =SUM1~—(.577215665+LOGF (Z72- ))*B(1)) 
BK (2)=2¢(10e/2Z)-RK(1)48B8(2))/BC1 ) 
F=BK (1) 

G=BK (2) 
BR(S)=((1le/Z)-BK(2)*B(3))078 (2) 
H=BK(3) 

END 


FUNCTION PMODF (EsRO+YsRHOsSIGMATR) 

THIS CALCULATES THE ESCAPE PROBABILITY OF NEUTRONS INCIDENT ON 
THE MODERATOR FROM THE FUEL. IT 1S CALCULATED FROM PSLOF By A 
RECIPROCITY RELATIONSHIPs SEE KeBe CADY THESIS* PG 886 
SIGMATR IS THE TRANSPORT CROSS SECTION OF THE MODERATOR 
FLAMBDA IS THE LINEAR EXTRAPOLATION LENGTH 

SIGMAMA 1S THE ABSORBTION CROSS SECTION IN THE MODERATOR 
R1,=RO*Y 

FLAMBDA= (6234 (1 0t2c048*SIGMATRERO ) )+07104 

SI GMAMA= e656 *RHO* 0.6024 7*SQRTF ( 20253/7E )418202 
D=l1o0/(3e*SIGMATR) 

DELTA=FLAMBDA/SIGMATR 

FK=SQRTF(SIGMAMA/D ) 

RK1=FK¥R1 

RKO=FK*RO 

CALL BESS(RKOsFIOOsFI10sFI20sFKOOsFK10sFK20) 

CALL BESS(RK1sFIOLsFI11 °F 121 sF KO] sFK11sFK21 ) 
B=(FK1L1*FLOO+F111*FKOO)/S(FKIO*FI1117FK11*F 110) 
C=B+FK*DELTA 

KEaPe/(C#HROK( (VY #¥2)—1 0) ¥FK) 

PMODF=1 0-26. *SIGMAMA¥RO* ( (Y¥¥2)—-10) *¥XE 

END 


FUNCTION PSLOF (E»ROsYsRHO>sSIGMATR) 

THIS CALCULATES THE ESCAPE PROBABILITY OF NEUTRONS APPEARING 
AT ENERGY (F&F) DUE TO COLLISIONS AT OTHER ENERGIESe DIFFUSION 
THEORY IS USED WITH A SPATIALLY UNIFORM SOURCE AND AN ENERGY 
DEPENDENT LINEAR EXTRAPOLATION LENGTHe 

FLAMBDA= (2662374 (1 ¢6+2048*®SIGMATR¥RO ) )+27104 

SI GMAMA= 2 656*RHO¥ e 6024 7* SQRTF ( 00 253/7E )/418002 

R1,=RO*Y 

D=10/4(3e*XSIGMATR) 

DEL TA=FLAMBDA/SIGMATR 

FK=SQRTF (SIGMAMA/SD) 

RK1=FK#¥R1 

RK O=FK*¥RO 

CALL BESS(RKOsFIOOsFIL1O0sFI2OsFKOOsFKI1OseFK2O ) 

CALL BESS(RK1sFI0O1sF111°9F121 »F KO] sFK11sFK21 ) 
B=(FK}Y1*F1LOO+F111*F KOO)/(FK1O*#FI14-FK11*F110) 

C=B+FK*DELTA 


tt 





tt 
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YO” 


LS 


PSLOF=2e/(C#RO#X( (Y#HD )~] 0 )*FK) 
END 


FUNCTION TRAP(Ss.DUeN) 

DIMENSION $(56)+DU(56) 

USES THE TRAPEZOIDAL RULE TO INTEGRATE THE N+1 DIMENSIONAL 
VECTOR Se OU IS THE LETHARGY INCREMENT © 

C=0e0 

DO 2@ IT=1.4N 

C=C+eS¥ (S01 )4S5S 0141) )* DUCT) 

TRAP=C 

END 


FUNCTION SCROSS(E+TsA) 

CALCULATES SCATTERING CROSS SECTION/FREE ATOM CROSS SECTION FROM 
THE MONATOMIC GAS MODEL. EF AND T ARE IN EVa A IS IN NEUTRON 
MASSESo 

Y=EXAST 

IF (Y¥~-700e) 14192 

X=SORTF CY) 
SCROSS2=(10+e¢5/Y) *F RF (X)4+10/4(X*EXPF (CY) *SQRTF (3214159) ) 

RE TURN 

SCROSS=(1e+05/Y) 

END 


SUBROUTINE WATER(CE+T«eNeSK) 

WATER CALCULATES THE NELKIN WATER KERNELe6 FUNCTIONS ONE AND TwO 
ARE NECESSARY o 

DIMENSION F (55 )+6SK (55.55) sWGT(12)+SKA(55) oKN(5).1M(5)+.AMU(12) 
OM1=0e205 

OM?P=0 e481 

OMR=0 606 

EM=186 

EMR=2232 

EMY =Bo/(lemilce/1Ge-1e/2e32) 

FORMAT (1H1) 


NM=N+ 1 

KN(1)=1 

KN(2)=11 
KN(3)=21 
KN(4)=31 
IM(1)=10 
IM(2)=20 
IM(3)=30 
IM(4)=40 


AMU(C1 )=0e9602899 
AMU (2 )=0° 7966665 
AMU(3)=025255324 
AMU (4 )=021834346 
AMU (5 )= 098940093 
AMU (6) =294457502 
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¥50 


80 
Cis) 
301 
300 


S02 
304 


Sls 


AMU (7 ) =e 86563120 
AMU (8 )=0 75540441 
AMU (9 )=0261787624 
AMU(10)=0e45801678 
AMU(11)=e28160355 
AMU(12)=0e09501251 
WGT C1 )=021012285 
WGT (2 )=e2223810 
WG6T(3)=03137066 
WGT (4 )=e3626838 
WGT (5 )=202715246 
WGT (6 )=006225352 
WGT (7)=2e09515851 
WG6T(8)=012462897 
WGT (9 )=014959599 
WGT(10)=016915652 
WGT(11)=0e18260342 
WGT€12)=218945061 
DO 102 I=1+eNm 

DO LOY J=1,] 
BEE=E(1)+E (J) 


C=SQRTF CECI )*¥E CJ) ) 


D=Oo 

H=06 

NJ=NMtJ 
NA=eS*FLOATF (NJ) 


IF CI-NA) 150150280 


KM=5 
ReHle 

Go TO 306 
KM=2] 

KL=4 


fete ¢J)—~0.33) 30043004301 
Meet )>o 34) 302330243038 


Mes | 

GO TO 304 

MJ=3 

OM=OMR 

AM=EMR 

LJ=3 

CM=1¢e¢/EM 

GX=T/EM 
B=EXPF (OMS (2 e*¥T) ) 


185 


At=((B*¥*2)4+1.)/ 00 (B¥*2 )-1 6) *EMREOMR ) +] oS (EMY*OM, 1420/7 (EMY*OMO ) 


Gey LO 305 
OM=0M 1 

AM=EMYV 

MIJ=e3 

LJ=] 
CM=10/EM+106/EMR 


GX=T/EM+OMR¥(1 oA CEXPF COMRST )-1 0) +05)/EMR 


B=EXPF (OM/(2c%#T) ) 


AT=10/(EMV*¥OM1 ) +207 (EMV#OM2 ) 
IF (I~J) 901.902.4901 





- 
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902 


901 


903 


ey A 


B02 


70 


10] 


1004 
BCS 


102 


| Shey 


oi ey 


eal 
200 


220 
2276 


= 


ALPHA=E (J)/SQRTF(30614159*GxX ) 

BE TA=E (J)*SQRTF(160/(3¢14159%T ) ) 

GO TO 903 

ALPHA=0.0 

BETA=O06 

ZOX=20* BS (AM¥OME (B¥¥ D1 ©) ) 

DO 151 KSKMeKL 

X=BEE+2 o¥C#¥AMU (K ) 

XX=SQRTF (X ) 

Z=X*¥ZO0X 

F=ONFE (CE (JY) sECI)eXeGXeOMsCMeMJ,LI,A1.B+Z,EMV) 
F=(Pot™Ble5%4o0 )*(F-ALPHA/SXX )+(4024/46)* 
(TWOCE CY) SE CI) eXeT )-BETAZXX ) 
X=BEE-2 eo *C¥AMU (K) 

Mx = SQRTF (CX) 

Z=X*ZOX 

G=ONF (CECI) es ECT 9 eo Xs GK se OMe CMeMIseLdIe Al eRe ZeEMV ) 
G= (2 .*81 05/40 )*¥ (G7~ALPHAS XX ) + (4024/44 )* 
(TWOCE (IY) eE CIT) oXeT) mBETA/XX ) 
D=D+WGT(K)*#(F+G) 
H=H+WGT (K)# AMU(K)* (GF ) 

IF CI-J) 70.802270 
FA=((20*81.5/40)* ALPHA+ (4024/46 )*BETA) *20/SQRTF(E(J)) 
D=D+FA 

H=H+tFA/36 

SK (J-1)=D 

SKA(J)=SK(deI) 

PRINT 1004,1 

PORMAT (4HI = 13) 

FORMAT (19H SCATTERING KERNEL /(7(13+F13¢e4)/)) 
PRINT 10054 (4U+SKA(J)eJ21¢!I) 

CONT INUE 

DO 110 J=2eNM 

JK=J~-1 

Be 110 la=leJK 

pipet (1 )7E (J) #2) * EXPE ((E CIE tho) 
SK (J.1)=SK(1I.,.J)*BB 

DO 951 L=144 

NI=KN(L) 

MI=IM(L) 

PRINT 199¢T 

FORMAT (25HISCATTERING KERNEL>+ KT = F7e4) 
PRINT 91994 CIe«I=NI «MI ) 

FORMAT (7Xe10(1299X%)) 

PRINT POOs (Te (SKCIYeT Pe IJeNIoM] ) eo fay eNM) 
FORMAT (1492X%010(E11063)) 

Naa 

MI=45 

PRINT 199¢T 

PRINT 91994¢(Te¢I=SNIoMI ) 

PRINT pPeOe (14 (SKt( Ie] )eIJ=NIoMI ) eo T=4 ¢NM) 
FORMAT (1442Kse5(E1163)) 

CONTINUE 

END 
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FUNCTION ONE (EPs Ee Xe GK sOMeCMeMIJeL Je Al seBeZeEMV) 


é THIS CALCULATES (SIGM~A(UP TO UsMU)y/(SIGMA BOUND ATOM74)) FOR 
Se NELKIN PROTONe EPsE+Xe AND T ARE IN EVe 

DIMENSION BI (200) +AA(3)*+CC(3) 

OM1=02205 

OM2=0.481 

AA(1)=1.6 


AAC2)=X7Z(EMV*OM1 ) 
AA(3)=(AA(2 FFE) /2 0 
CC(1)=1 o 
CC(2)=26*®XK/S (EMV¥OME ) 
G6E(3)=(CC(2)¥e2) 7/26 
PROD1=0e0 
SUM=0 6 
D=E*SQRTF (EZ (EP®X¥GK¥3014159)) 
N=10e¢5*7+100 
Se aCN)) =O'> 
SUM=1 00E~20 
BI (N-1)=SUM 
NN=N=1 
DO 11 J=2@eNN 
T=N-J 
BIC] )=2@ce*#FLOATF (I )*BICI4+1)/Z24+BI (6 1+2) 

11 SUM=SUM+BI (1) 
SUM=2 0*® SUM-BI (1 ) 
C=EXPF (Z)7SUM 
DO 12 J=1>5NN 

12 BI(J)=C#XBI (J) 
DO 1050 M=1 «MJ 
FM=M=1 
PROD=0 0 
DO 17 L=leLJ 
tal 1 
SUM=0 6 
DO 13 N=leNN 
EN=N=— 1 
BB=E-EP+X*CM+tENXOM+EL*XOM] +FM*OM2 
AZ= (BB¥*2)/(40*X*GX) 
ie CAZ=7058) 112091120113 

112 P=( (B**EN)*BI O(N) )ZEXPFCAZ) 
GO TO 114 

113 P=0.0 

114 SUM=SUM+P 
Ww=P/SUM 
[IF (W-1c00E-4) 14314413 

re CONTINUE 
PRINT SOO—N 

S00 FORMAT (SH N = 13) 

14 DO 15 N=2eNN 
EN=N-~1 
BB=E~-EP+X*CM-ENXOM+EL*OM1 +F M#OM2 
AZ=(BB¥*2)7(4e*X*GX ) 
Gears) 19591159116 

115 P=(BIEN)/(B¥*EN) )/EXPF (AZ) 
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116 
ete 


16 
17 


oS 0 


6 


1.8 


CO (HOw le 7 

P=0e0 

SUM=SUM+P 

W=P/SUM 

Pe Weite Oe 4) 1 65 Gat 
GENT I NUE 

PRINT SOON 

PROD=AA (L)*SUM+PROD 
CONT INUE 
PROD1=CC(M)*PROD+PROD}1 
CONT I NUE 

ONE =D*PROD1 /EXPF (X#AT ) 
END 


FUNCTION TWOC(EPsE«XoT) 
THIS CALCULATES (SIGMA(UP TO UsMU)7(SIGMA BOUND ATOM74) 4 
FOR A FREE GAS OF MASS 16+ FPsE€sX+sAND T ARE IN EVe 
B=E-EP+X/166 
D=E*SQRTF(E*160/(X¥T¥3014159%EP) ) 
TWO=D/EXPF ( (B¥¥2)*160/(Ge*X¥*T) ) 
END 

END COUTH 


SUBROUTINE FKERN(ONM.SFAs Ts As Es SK) 

SUBROUTINE FKERN CALCULATES THE SCATTERING MATRIX USING THE 
MONATOMIC GAS MODEL. 

DIMENSTON F(55) eSK(55955)9W(55 )9+RXKX(55)+EX(55) 
ETA=(At10)/(2e*SQRTFCA)) 

Rero=(A—te)/(2o*SOQRTE CA) ) 

DO 43 I[=1+NM 

Meet )=SQRTF CECI )ZT } 

RX(IT)J=RIO*FWCI ) 

ao? P=ETAW (I) 

DO 10 JT=ie«eNM 

Porvoo” J= 15 | 

SK (Se Tee S*SFAX(E TARE) RCFE CT) SECI) FC ERE (CEX CT) emRx (OY) +E RE OCEX CI 4+ 


1RX(4S) + CEXPE C CWI) #8 y— (WOT) %% 2B) ) eC ERE (EX CJ)-RX 01) ERE ORX (1 4 
PEX(I)))) 


END 
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APPENDIX F 
NOTATION 


Fuel rod radius in mean free path, a = Ser 


Fuel element area, Eq. (2.12) 

Ratio of scatterer mass to neutron mass, 
EG eee eok ai) 

Parameter for substitution into the Nelkin 
Water Kernel, Eqs. (5.17) and (5.18) 

Mass number, Eq. (7.3) 

Ratio of the scattering to the total cross 
section, p. 2/7 

Correction factors to the diagonal terms in 
the Nelkin Water Matrix, p. 59 

Eg on (4e750)) 

Bae (eo) 

Linear extrapolation length into the fuel 
element, Eq. (2.18) 

Linear extrapolation length into a black hole, 
Bae ( Cc 


Diffusion coefficient for the homogenized cell, 


Eq. (6.38) 
Bae (43.0) 
Eq. (4.32) 
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DCE) 
D¢ (E) 


120 


Diffusion coefficient for the moderator, Eq. (2.13) 
Diffusion coefficient for the fuel, p. 2/7 

Neutron energy in eV 

Thermal cutoff energy 


Thermal utilization, Eqs. (6.13) and (6.15) 


Eqe (5-21) 

Eqs. (2.7)>) (2226). andar.) 
Eq. (4.26) 

Eq. (4.27) 

Eq. (4.28) 

Eq. (4.29) 

Sec. 4.3, Eq. (4.45) 

Eqmeet 2.457) 


Planck's Constant 

Modified Bessel function of the first kind 

Current into the fuel element, Eq. (2.24) 

Net current entering the cladding from the 
fuel rod, Eq. (4.22) 

Net current entering the fuel element, 
Ha. (2212) 

Modified Bessel function of the second kind 

Thermal lifetime, Eq. (6.17) 

Diffusion length for the homogenized cell, 
Be oso 

Eq. (4.33) 
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—2 
Ly 


xy 
Pa 
M 
M(E) 


Eu 


Edin (42549) 
Ea. (4.36) 
Eq. (4.39) 
Did 


Maxwellian distribution, Eq. (2.6) 

Cylindrical moments, Sec. 4.1 

Deg bd 

Dec, 

Spatially averaged source, Eq. (2.11) 

Neutron source, Eq. (2.11) 

Spatially averaged energy dependent neutron 
density in region i, Eq. (6.5) 

Total neutron density in region i, Eq. (6.7) 

Fuel element escape probability, Chap. 4 

Moderator escape probability, Sec. 3.2 

Fuel rod escape probability, Sec. 4.1 

Escape probability for the fuel rod and void 
region, Sec. 4.2 

Moderator escape probability, Sec. 3.l 

Eq. (4.43) 

EGE) 

Fuel element radius 

Equivalent cell radius 

Inner radius of the cladding 


Fuel rod radius 





Appendix F 


in 

out 
1 

RR. 


S(E) 
S(E) 


S(t, 5) 


<0 -s 


12 


Sec. 4.3, Eq. (4.49) 

Seem 4 (SRR EG = (4-50) 

Reaction rate of type j in region i, Eqs. (6.18) 
and (6.19) 

Neutron slowing down source, Eq. (2.3) 

Spatially averaged slowing down source of 
neutrons, Eqs. (2.10), (5.27), and (7.1) 

Spatially dependent slowing down source of 
neutrons, Eq. (2.10) 

EGuat4 24) 

Moderator temperature in electron volts 

Sec. 4.300Eq.. Geos 

Lethargy 

Neutron velocity 

Average neutron velocity in region i, Eq. (6.6) 

The volume of region i 

Weighting coefficients for the Gaussian 
quadrature integration formula, p. 50 

Eq. (4.42) 

Ihe) 5 wiles a4), 

Povo. y = r,/r 

Ears) 

Eqs. (5.17) and@@auley 

eh eee, 


O 
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Greek Letters 
a(E) 
B(E) 


u(E) 


Mk 


Me. 


Fuel capture to fission ratio 

Eq. (4.2) 

Euler's Constant 

Moderator neutron density disadvantage factor, 
Eq. (6.8) 

Cladding neutron density disadvantage factor, 
Eq. (6.9) 

Moderator flux disadvantage factor, Eq. (6.10) 

Cladding flux disadvantage factor, Eq. (6.11) 

Convergence criterion, Eq. (C.24) 

Eq. (4.20) 

Number of fission neutrons produced per thermal 


absorption in the fuel, Eq. (6.16) 


Eda .7) 
Hii (an2o)) 
p. 15 

De 28 
Bacar 
Eaie (4.315) 


Cosine of the scattering angle 

Average cosine of the scattering angle, Eq. (5.24) 

A zero of the Legendre polynomials, p. 50 

Number of fission neutrons produced per thermal 
fission 

HiGueeC4s 23) 

Eqs. (5.17) and (5.18) 
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194 


Sra 59 

Eq. (5.8) 

Sec. 4.35) Ede (4-545 

See, 4.3, Equa 

Bound atom microscopic cross section 

Free atom microscopic cross section 
Microscopic cross section expressed in barns 


for interaction type j in material i 


o(E'—=E,.w) Differential scattering kernel, Eqs. (5.1) 


o(E'—=E) 
o* (E—=E) 


0, (E'—=E) 


zB; (E) 


si 

J 
3(E'—=E) 
48 


op (E) 


ang (5415) 
Scattering kernel, Eqs. (5.4) and (5.19) 
Modified diagonals for the Nelkin water 
kernel, Eq. (5.26) 
First moment of the scattering kernel, 
Eaei(.:2 5) 
Macroscopic cross section expressed in an 
for interaction type j in material i 
Macroscopic effective cross section for inter- 
action type j in material i, Eqs. (6.23) and 
(6231) 
Macroscopic mean cross section for interaction 
type j in material i, Eqs. (6.24) and (6.32) 
Macroscopic scattering kernel, Eq. (2.1) 
Sec. 4.3, Eqv-(4eq0) 
Energy dependent flux in an infinite homogeneous 


moderator, Eq. (2.1) 
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01 (E) 


Subscripts and 


a 


a> Fh Fh 


tr 


IES 


Spatially averaged energy dependent flux in 
region i 

Spatially averaged 2200 m/sec flux in region i 

Spatially averaged thermal flux in region i 


th 


iteration for the moderator flux at 
lethargy u,, Eq. (C.22) 

Oy. (a, ) normalized, Eq. (C.23) 

Eq. (4.5) 

Eqs. (5.17) and (5.18) 

Eqs (2,20) 

p. 47 

p. 47 


Dey 


Superscripts 
absorption 
cladding 

fuel 

fission 
homogenized cell 
moderator 
scattering 

total 


transport 





WORKS CITED 


BOOKS 


. M. Case, F. de Hoffmann, and G. Placzek, Introduction to 
the Theory of Neutron Diffusion, Vol. I, United States 
Government Printing Office, 1953. (See Footnotes 28 

and 35.) 

. B. Dwight, Tables of Integrals and Other Mathematical Data, 
The Macmillan Company, New York, 1961. (See Footnotes 87 
and 90.) 

. Glasstone and M. C. Edlund, The Elements of Nuclear Reactor 
Theory, D. Van Nostrand Company, Inc., Princeton, 
Toronto, London, and New York, 1952. (See Footnotes 20 
and 25.) 

. M. Weinberg and E. P. Wigner, The Physical Theory of Neutron 
Chain Reactors, The University of Chicago Press, Chicago, 


1958. (See Footnotes 1, 23, 26, 30, 63, 76, and 83.) 


JOURNALS 


. Amouyal, P. Benoist, and J. Horowitz, ‘Nouvelle Methode de 
Determination du Facteur d'Utilization Thermique d'une 
Cellule," J. Nucl. Energy 6, 79 (1957). (See Footnotes 
LO@antcde 2 7-9) 

. R. Beyster, et al., "Measurement of Neutron Spectra in 
Water, Polyehtylene, and Zirconium Hydride," Nucl. Sci. 


Eng. 9, 168 (1961). (See Footnote 54.) 


~196- 





Works Cited 197 


K. B. Cady and M. Clark, Jr., “Neutron Transport in Gyiendmieat 
Rods,'' Nucl. Sci. Eng. 18, 491-507 (1964). (See Footnotes 
2a. anders 2.) 

K. B. Cady and C. R. Mac Vean, "A Simplified Formulation of 
Space-Energy Cell Theory,'' Trans. Am. Nucl. Soc. 6, 234- 
235 (1963). (See Footnote 70.) 

G. P. Calame, "Diffusion Parameters of Water for Various 
Scattering Kernels," Nucl. Sci. Eng. 19, 189 (1964). 

(See Footnote 55.) 

L. G. de Sobrino and M. Clark, Jr., "A Study of the Wilkins 
Equation,'' Nucl. Sci. Eng. 10, 388 (1961). "Comparison 
of the Wilkins Equation with Experiments on Water Systems," 
Nucl. Sci. Eng. 10, 377 (1961). (See Footnote 47.) 

M. Goldstein and R. M. Thaler, "Recurrence Techniques for 
Calculation of Bessel Functions,'' Mathematical Tables 
and Other Aids to Computation, 13, 102 (1959). (See 
Footnotes 88 and 89.) 


A. Hassitt, ''Methods of Calculation for Heterogeneous Reactors," 


Progress in Nuclear Energy, Series I, Vol. 2, Pergamon 
Press, New York, London, Paris, and Los Angeles, 1958. 
(See Footnote 4.) 

R. L. Hellens and H. C. Honeck, "Summary and Preliminary 
Analysis of the BNL Slight Enriched Uranium, Water 
Moderated Lattice Measurements," Light Water Lattices, 
International Atomic Energy Agency, Vienna, 1962. 


(See Footnote 75.) 





—S>S ti‘ re . i 


=P 
- = 
ee 
— | 8 = 
_ 
== 





Works Cited 198 


H. C. Honeck, "The Calculation of the Thermal Utilization and 
Disadvantage Factor in Uranium/Water Lattices," Nucl. Sci. 
Eng. 18, 49-68 (1964). (See Footnotes 7, 15, 38, 73, 74, 
77, 78, 79, and 81; Table A.3; and Figs. 6.136 -2 sone 
Soe) 

A. N. Lowan, N. Davis, and A. Levinson, "Tables of the Zeros 
of the Legendre Polynomials of Order 1-16 and the Weight 
Coefficients for Gauss's Mechanical Quadrature Formula," 
Bull. Am. Math. Soc. 48, 739 (1942). (See Footnote 57.) 

M. S. Nelkin, ''Scattering of Slow Neutrons by Water,'' Phys. Rev. 
119, 741-746 (1961). (See Footnotes 13, 41, 46, 48, and 
oo) 

C. D. Petrie, M. L. Storm, and P. F. Zweifel, "Calculation of 
Thermal Group Constants for Mixtures Containing Hydrogen," 
Nucl. Sci. Eng. 2, 728-744 (1957). (See Footnote 64.) 

G. C. Pomraning, "On the Energy Averaging of the Diffusion 
Coefficient,'' Nucl. Sci. Eng. 19, 250 (1964). (See 
Footnote 67.) 

G. C. Pomraning, ''On the Validity of the Constant Source 
Assumption for the Cell Problem," Nucl. Sci. Eng. 18, 
400-403 (1964). (See Footnote 21.) 

M. J. Poole, M. S. Nelkin, and R. S. Stone, ''The Measurement 
and Theory of Reactor Spectra," Progress in Nuclear 
Energy, series I, Vol. 2, Pergamon Press), New YORK: 
London, Paris, and Los Angeles, 1958. (See Footnotes 


Joeem 51, and 91. ) 





= -&- 


——— — om 





Works Cited 199 


K. S. Singwi and L. S. Kothari, “Transperes@reccescce tenet 
Thermal Neutrons in Solid Moderators,'' Second Geneva 
Conference Proceedings 16, 325 (1958). (See Footnote 
Ge) 

N. G. Sjostrand, "Definition of the Diffusion Constant in 
One-Group Theory," J. Nucl.Energy, Part A, Reactor 
Science 12, 151 (1964). (See Footnote 66.) 

E. Starr and J. Koppel, "Determination of Diffusion Hardening 
in Water," Nucl. Sci. Eng. 14, 244 (1962). (See 
Footnote 56.) 

G. W. Stuart and R. W. Woodruff, "Method of Successive 
Generations,'' Nucl. Sci. Eng. 10, 388 (1961). (See 
Footnote 33.) 

R. S. Varga, "Numerical Methods for Solving Multi-Dimensional 
Multigroup Diffusion Equations," Proceedings of Symposia 
in Applied Mathematics, Vol. XI, Nuclear Reactor Theory, 
American Mathematical Society, Providence, 1961. (See 
Footnote 16.) 

D. F. Zaretsky, "Effective Boundary Conditions for Grey 
Bodies,'' First Geneva Conference Proceedings 5, 525- 


530 (1955). (See Footnote 84.) 


REPORTS 


D. D. Clark, ''The Cornell University Nuclear Reactor Laboratory," 
CURL-1, Cornell University (1961). (See Footnote 14 and 
68 and Table A.2.) 





a < s = —- oS =) 








Works Cited 200 


F. D. Federighi and D. T. Goldman, ‘'KERNEL and PAM-Programs 
for Use in the Calculation of the Thermal Scattering 
Matrix for Chemically Bound Systems,'' KAPL-2225, Knolls 
Atomic Power Laboratory (1962). (See Footnotes 49, 58, 
and 86.) 

D. Hughes and B. Schwartz, ''Neutron Cross Sections," BNL-325, 
2nd Ed., Brookhaven National Laboratory (1958). (See 
Footnotes 60 and 72.) 

H. C. Honeck, ''THERMOS-A Thermalization Transport Theory Code 
for Reactor Lattice Calculations," BNL-5826, Brookhaven 
National Laboratory (1961). (See Footnotes 6, 17, and 
59.) 

D. C. Leslie, "'The Calculation of Thermal Spectra,'' Proceedings 
of the Brookhaven Conference on Neutron Thermalization, 
BNL-719, Vol. II, pp. 592-609, Brookhaven National 
Laboratory (1962). (See Footnote 9.) 

D. C. Leslie, "The SPECTROX Method for Thermal Spectra in 
Lattice Cells,'' AEEW-M211, United Kingdom Atomic Energy 
Authority (1962). (See Footnotes 8, 19, 22, and 85.) 

Neutron Spectra in Lattices and Infinite Media, Proceedings 
of the Brookhaven Conference on Neutron Thermalization, 
BNL-719, Vol. II, Brookhaven National Laboratory (1962). 
(See Footnotes 3 and 5.) 

A. Radkowsky, ''Temperature Dependence of the Thermal Transport 
Mean Free Path,'' Physics Quarterly Report (April, May, 
June, 1950), ANL-4476, p. 89, Argonne National Laboratory 
(July, 1950). (See Footnote 42.) 








Works Cited 201 


The Scattering Law, Proceedings of the Brookhaven Conference 
on Neutron Thermalization, BNL-/719, Vol. I. Brookhaven 
National Laboratory (1962). (See Footnote 44.) 

E. P. Wigner and J. E. Wilkins, Jr., "Effect of the Temperature 
of the Moderator on the Velocity Distribution of Neutrons 
with Numerical Calculations for H as the Moderator," 
AECD-2275 (1944). (See Footnotes 11, 39, 45, and 52.) 

J. E. Wilkins, Jr., "Effect of the Temperature of the Moderator 
on the Velocity Distribution of Neutrons for a Heavy 
Moderator, CP-2481 (1944). (See Footnotes 12, 40, 45, and 
a0.) 

C. H. Wescott, "The Specification of Neutron Flux and Effective 


Cross Sections,'' AECL-352 (1956). (See Footnote 61.) 


THESES 


S. S. Berg, ‘Initial Experiments on the Cornell University Zero 
Power Reactor Cores," Cornell University Thesis (1964). 
(See Footnotes 69 and 80 and Table A.1.) 

K. B. Cady, "Neutron Transport in Cylindrical Rods," 
Massachusetts Institute of Technology Thesis (1962). 

(See Footnotes 29, 31, 34, 36, and 37.) 

J. A. Larrimore, ''Temperature Coefficients of Reactivity in 
Homogenized Thermal Nuclear Reactors,'' Massachusetts 
Institute of Technology Thesis (1962). (See Footnote 62.) 

J. Suich, ''Temperature Coefficients in Heterogeneous Reactor 
Lattices,'' Massachusetts Institute of Technology Thesis 


(1963). (See Footnotes 71 and 82.) 





a 
— 
— 
7 7 
—_ 
- _— 
———>- al —- 



































“AN 





